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Abstract 

We use a prescription to gauge the SU{2) Skyrme model with a U{1) field, 
characterised by a conserved Baryonic current. This model reverts to the usual 
Skyrme model in the limit of the gauge coupling constant vanishing. We show 
that there exist axially symmetric static solutions with zero magnetic charge, 
which can be electrically either charged or uncharged. The energies of the 
(uncharged) gauged Skyrmions are less than the energy of the (usual) ungauged 
Skyrmion. For physical values of the parameters the impact of the U{1) field 
is very small, so that it can be treated as a perturbation to the (ungauged) 
spherically symmetric Hedgehog. This allows the perturbative calculation of 
the magnetic moment. 
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1 Introduction 



For a long time now, much attention has been paid to the Skyrme model in 3 
dimensions. It is believed to be an effective theory for nucleons in the large limit 
of QCD at low energies. The classical properties as well as the quantum properties 
of the model are in relatively good agreement with the observed properties of small 
nuclei [|, I, |]. 

Gauged Skyrme models have been used in the past. The U{1) gauged model 0, ^ 
was used to study the decay of nucleons in the vicinity of a monopole 0, while the 
SU{2)l gauged model was used to study the decay of nucleons when the Skyrme 
model is coupled to the weak interactions 0. The Skyrme model has also been used 
to compute the quantum properties of the Skyrmion where the gauge degrees of 
freedom are quantised to compute the low energy eigenstates of a Skyrmion. These 
states were identified as the proton, the neutron and the delta. 

The aim of this work is to show that the Skyrme model can be coupled to a self 
contained electromagnetic field and that this U{1) gauged model has stable classical 
solutions. In addition to these solitons with vanishing magnetic and electric flux, we 
show that this system supports solutions with nonvanishing electric flux which are 
analogous to the dyon solutions of the Georgi-Glashow model, just as the uncharged 
solitons are the analogues of the monoploles [§] of that model. The electrically charged 
lumps have larger energy, or mass, than the uncharged soliton, just like the Julia- Zee 
dyon [P] has larger energy, or is heavier, than the (electrically uncharged) monopole. 
We shall refer to these lumps as charged U{1) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Maxwell gauged Skyrme 
model, the present work is also an example of a soliton in a (i-dimensional SO{N) 
gauged S"^ model with N < d for the case d = 3, N = 2 , extending the results of 
Ref. which were restricted to the N = d cases. (The work of Ref. |TU| consists 



of establishing topological lower bounds for the generic case, encompassing earlier 
examples in two |TT] and three [|T^, 13] dimensions respectively.) The gauging pre- 



scription used here by us coincides precisely with that used in Ref. and permits 
the establishing of a topological lower bound which did not feature in Ref. |^ and 
which is carried out here to establish the stability of the soliton. Such lower bounds 
are absent in the other prescription of gauging the Skyrme model as in Refs. 0]. 
(Notice that we name the sigma models after the manifold in which the fields take 
their values rather than using the name of the symmetry group for the model. Thus 
what is sometime called the 0{d+ 1) model in the literature will be refered to as the 
S"^ model.) 

The U{1) gauged SU{2) Skyrme model is described by the Lagrangian p 

C = ^Tr{D,UD,U^) - ^Tr{[iD,U)U\ {D.UW^f - 

where the U{1) gauge covariant derivative is 

D^U = d^U + ieA^[Q,U], (1) 
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where T^l,v = d^Aj, — d^A^ and where the charge matrix of the quarks is expressed 
as Q = + Ca). This differs from the covariant derivative of Ref. only in the 
unimportant matter of the sign of i in (|1]), which we have chosen for consistency of 
the convention used in Ref. [|10|] . 

In what follows it will be more convenient |jlO| to parametrise the Skyrme field as 
an valued field 0" = (0", (j)^), a = 1, 2, A = 3, 4 subject to the constraint = 1. 
The two fields U and (j) are related to each other via the following expression 

U = (ly^r", U-^ = U^ = (2) 

where r" = {ia°',ia^,l) and = {—ia",—ia^,l), in terms of the Pauli matrices 

The gauge covariant derivative now can be re-expressed as 

= d,r + D,^^ = d,^^. (3) 

where = eAf, and F^,y = eT^u- 

The Lagrangian for the U{1) gauged Skyrme model can then be written as 

C = -XoFl^ + \i\D^r? - \2\D^^rD.^(t>'? (4) 

where the square brackets on the indices imply (total) antisymmetrisation and where 
Aq ^ = 4e^, Ai = and A2 ^ = 8a^. The late Greek indices /z label the Minkowskian 
coordinates, while the early Greek indices a = 1,2 and the upper case Latin indices 
A = 3, 4 label the fields 0" = (0", 0^). 

The static Hamiltonian pertaining to the Lagrangian (^) is 

n = x^Ff^ + A1IA0I' + A2|/^[.0"/^,]0T 

(5) 

+z2Ao|9,Aop + A2{Ai|0-|2 + 16A2[|012|9,0^P + ||9.(|0^p)|2]}, 

where the indices i = a, 3 label the space- like coordinates. 

To find the static solutions of the model, one would usualy solve the Euler La- 
grange equations which minimise the Hamiltonian @, but because of the electric 
potential Aq, one must solve the Euler Lagrange equations derived from the La- 
grangian @). We then look for static solutions, but, as for the Julia-Zee dyon 0, we 
have to impose the proper asymptotic behaviour for the electric potential to obtain 
static solutions which are electrically charged (in the classical sense, i.e. solutions 
where the flux of the electric field is non zero). 

When the full equations of motion are written down, one finds as expected that 
there are static solutions for which Aq = 0, i.e. solutions for which the electric field is 
identically zero. For these solutions in the temporal gauge, the equations of motion 
reduce to the equations obtained by minimising the Hamiltonian (^. We study the 
solutions of unit Baryon charge of the U{1) gauged Skyrme model with and whithout 
an electric field, for various values of the U{1) coupling constant (or equivalently the 
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Skyrme coupling). For physical values of these parameters in the model, we find that 
the energy (mass) of the gauged Skyrmion does not differ significantly from that of 
the ungauged charged-1 Skyrmion, namely the familiar hedgehog This implies 
that for these values of the physical parameters, the U{1) gauged Skyrmion can be 
regarded as a perturbation of the (ungauged) hedgehog, enabling the computation 
of the magnetic moments of the gauged Skyrmion (i.e. the Neutron) and the shift 
of the energy of the gauged Skyrmion away from the energy of the hedgehog 
perturbatively using the method employed by Klinkhamer and Manton |T^ for the 
sphaleron of the Weinberg-Salam model. 

In Section 2, we define the topological charge and establish the corresponding 
lower bound on the energy functional. In Section 3 we present the solutions which 
have no electric fields in the first subsection and electrically charged solutions in 
the second subsection. The perturbation analysis of the gauged Skyrmion around 
the (ungauged) hedgehog is carried out in Section 4, and Section 5 is devoted to a 
discussion of our results. 



2 The topological charge and lower bound 

The definition of the topological charge is based on the criterion that it be equal to the 
Baryon number, namely the degree of the map. For the gauged theory however, this 
quantity must be gauge invariant as well. This requirement can be systematically [|TU 
satisfied by arranging the gauge invariant topological charge density to be the sum 
of the usual, gauge variant winding number density 

Qo = ei,ke'''"'dirdj4>'dk4>'4>' , (6) 

plus a total divergence whose surface integral vanishes due to the finite energy condi- 
tions, such that the combined density is gauge invariant. In 3 dimensions, this is given 
explicitly in Refs. |jTO|, |l3l for the SO (3) gauged model, and for the present case of 
interest, namely the SO (2) gauged S^ model, the charge density can be derived from 
that of the SO (3) gauged model by contraction of the gauge group SO (3) down to 
5*0(2). It can also be arrived at directly. To state the definition of the charge, we 
denote the gauge covariant counterpart of (|) by 

Qg = e,,ke'''"'D4'^DjCt>'Dj,ct>^ct>'' , (7) 

so that using the notations (^ and (0) we have the definition of the gauge invariant 
topological charge 

Q = Qo + diD^i, (8) 
= QG + le^,kF.,{e^''(j)''Dk(j)^). (9) 

In (H) the density Qi is the following gauge variant form 

= 3e,,ke^''Ajd,<j)^ 0^. (10) 
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The flux of fij vanislies, as can deduced by anticipating tlie finite energy conditions 
to be stated later. 

Note that the 3-volume integral of go in ® is the degree of the map for the 
ungauged system namely the Baryon number. 

Identifying g, (P), with the naught component j° of the Baryon current, j'^ is 
defined by 

= e^^""" BabcdD^r Dp<f)''D^<f)'<f)^ - ^ e^"''' F,p{eAB<p'' D^^p^). (12) 

The 4-divergence of (|ll]) receives a contribution only from its first term, which 
being locally a total divergence implies that the 3-volume integral of is a conserved 
quantity. Alternatively we consider the 4-divergence of (p!2D , 

d,f = Q6^''P''e^peABD^rDu<P''Dp<p^D,<p'' (13) 

which is analogous to the corresponding quantity in the work of Goldstone and 
Wilczek 0] . This contrasts with the expression for the total divergence of the topolog- 
ical current in the work of d'Hoker and Farhi 0, where a different gauging prescription 
is used leading to that quantity being equal to the local anomaly. 

We now proceed to find a model whose Hamiltonian Tio is bounded from below by 
the topological charge density defined by (P). We will then show that the Hamiltonian 
(|^) is given by Ho plus certain positive definite terms. 

First of all, we reproduce the density qg, (0), in by using the following in- 
equality 

(KsD.r - e,,ke'""'^4D,(l)'Dk(l)''(l)^)'' > (14) 

where the two constants K3 and K2 have the dimensions of length. Expanding the 
square, we get qg on the right hand side of 

KliD^r)^ + 4{D[,rDj]<p'? > 2kskIqg. (15) 

To reproduce the other term in (^, ^EijkFij^e^^ (j)^ dk4>^) , we use the following 
inequality 

(nlFij - ^K^e,^„e^^(f)^Dk(j)y > (16) 

yielding 

4F'j + /«lJ(e^^0^A0^)' > «:^/€4£.,fcF,,(£^^0^Z^fc0^). (17) 

With the special choice for the relative values of the constants Sh^k^ — '^4'^0' 
sum of (0) and (|T5|) yields the following 

4F^, + 4iD^rr + 4{D^^rD,^ct>'f + ^(^^^^^a^^)^ > 2^,4^. (is) 
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The right hand side of (^) is now proportional to the topological charge density 
Q defined by so that the inequality ([T8|) can be interpreted as the topological 
inequality giving the lower bound on the energy density functional if we define the 
latter to be the left hand side of (plSf), namely 



2, A 



4/^0 



(19) 



The Hamiltonian system (|T9D is almost the Hamiltonian of the gauged Skyrme 
model (^ (remember that Aq = 0). It differs from the latter only in its last term. 
Now we can use the identity 



(e 



ABiB 



Di 



kA\2 



(A 



ia\2 



+ (0["A 



and add the positive definite term 



9^ 



i(0[-A0^l)' + (0f"A0^l)' 



(20) 



appearing on the 



right hand side of (0) to Tio in (IH) to end up with the Hamiltonian for the U{1) 
gauged Skyrme model: 



n 



(21) 



which is nothing but the static Hamiltonian (^) in the temporal gauge Aq = 0, and 
where 



^0' 



A, 



Kjr) 



(22) 



By virtue of ([T8|), (pT)) is also bounded from below by namely by a number 

proportional to the topological charge density g. 

We thus see that Hq can be considered as a minimal (f/(l) gauged) model, but 
from now on, we will restrict our attention to the physically more relevant model (pl|) 
and integrate it numerically to find its topologically stable finite energy solitons. 

The soliton solutions to the system (|2l|) can only be found by solving the second- 
order Euler- Lagrange equations, and not some first-order Bogomol'nyi equations since 
saturating the inequalities (|1^) and ([T6|) would not saturate the lower bound on the 
energy density functional Ti. In this context we note that saturating (14) and (|T^) 
does indeed saturate the topological lower bound on the functional Hq by virtue of 
the inequality (0), and should it have turned out that the Bogomol'nyi equations 
arising from the saturation of (|1^) and (|T6D supported non-trivial solutions, then TCq 
would have been a very interesting system to consider. As it turns out however, these 
Bogomol'nyi equations have only trivial solutions in exactly the same way as in the 
case of the (ungauged) Skyrme model |jl|. 

The energy for the static configuration, when the electric field vanishes, is ex- 
pressed as 



^(Ao, Ai, A2) = J d'x [Ao4 + Ai(Ar)' + 



kb\2-l 



(23) 



6 



and performing the dilation x ax, ^ a ^A^, we get 



E{Xo, Ai, A2) = / d'x [^Ff^ + aAi(A0'^)' + -j{Dy,rD,^<p'f]. (24) 



If we choose a = (^)^^^ then we have 



E{\,,K\,) = {\,\,f/'E{^^^) , (25) 

from which we see that we can set Ai = A2 = 1 without any loss of generality. By 
virtue of (^TJ) and (P5|), we can finally state 



E{\o, Xi, A2) > 2{-^y/' / d'x Q. (26) 



A1A2 

H 

Notice that for the usual Skyrme model we have 

Esk{K A2) = jdx^[\Mr? + Hdyird,^(t>''?] 

= (AiA2)i/2^.fc(l,l) (27) 
> 2(AiA2)i/2j^^3^^_ 

We will use (pT]) to compare the numerical solutions of the gauged Skyrme model 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability considerations discussed 
in this Section apply only to the solutions with no electric field, i.e. with Aq = 0. 



3 The soliton and the charged U{1) Skyrmion 

To find the static solutions, we have to look for the largest symmetry group of the 
functional to be subjected to the variational principle, and look for solutions which 
are invariant under that symmetry group. For the solutions in the Aq = gauge 
this is the static Hamiltonian (^), while for the solutions in the 7^ it is the 
Lagrangian (|[). For our choice of gauge group the largest symmetry is the 5*0(2) 
group corresponding to an axial rotation in space-time and a gauge transformation 



on the gauge field. Defining the axial variables r = y + and and z = X3 in terms 
of the coordinates = (xq,, X3), a = 1,2, the most general axially symmetric Ansatz 



|T5| for the fields 0" = (0", 0^) (with a = 1, 2 and A = 3, 4), and, A = (A„, A3), is 



0" = sin / sin g rf , = sin / cos g , = cos / , (28) 
a{r,z)-n ^ C2(r,^) , ci(r,z) h{r, z) 

Aa = eal3 Xf3 H Xa A3 = , Aq = , (29) 
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with = (sinn0, cosn0) in terms of the azimuthal angle and Xa = Xa/r. n in 71°" 
is the vorticity, which for the Nucleons of interest to us here, equals one, n = 1. The 
functions a, b, ci, C2, / and g both depend on r and z. 

Our Ansatz (|29D for the U{1) field consists of decomposing the latter in the most 
general tensor basis possible. We will find out below, when we compute the Euler- 
Lagrange equations, that the functions Ci and C2 vanish identically. Anticipating this, 
we suppress them henceforth. In its final form this Ansatz agrees with that used in 
Ref. |T5[, the latter being arrived at by specialising the Rebbi-Rossi Ansatz for the 
axially symmetric SO {3) field. 

The static Hamiltonian, z.e.the Too component of the energy mometum tensor T^,y, 
is then given by 



H 



+ al + {br - + b 



+ 



Ai 



/2 + /2 + sin^ f{gl + gl) + ^ sin^ / sin^ g 



+2A2 sin^ / 



{frg. - f.grf + s\n'g[^{n + P. + {qI + 9l) /)]] }rdrdz 



(30) 

The boundary conditions for the Skyrmion fields are the same as the boundary 
conditions for the hedgehog ansatz when expressed in the cylindrical coordinates 
where g = tt/2 + arctan(2;/r) and defining R = {z'^ + r'^Y^'^, f = f{R) with /(O) = ir 
and limij__»oo f{R) = 0. From this we can deduce that the function / has a fixed value 
at the origin and at infinity. For smothness along the z axis, each field, that is f, g and 
Aa must satisfy the condition that the partial derivative with respect to of the field at 
r = vanishes. The boundary conditions and the asymptotic behaviours for a and b 
are chosen so that the gauge fields are well defined and Aq looks asymptoticaly like 
a coulomb field (i.e. with an electric charge but no magnetic charge). We also require 
that the total energy be finite. These conditions leads to the following constraints: 



/(0,0) =71 /(r^oo,^^oo) =0 fr{r = 0,z) =0 

g{r = 0,z<0) =0 g{r = 0,z>0) = tc gR\R^oo =0 

a{r = 0,z) =1 a^loo =0 a^joo =0 

ar{r = 0,z) =0 y4o(r 00, 2; ^ 00) =Vo + q/r AQ{r = 0,z) =0 

(31) 

where R = {z'^ + 7-2^1/2 where we have used the notation |^ = a,- etc. Note 
that the field g is undefined at the origin and the resulting discontinuity of g at that 
point is an artefact of the coordinate system used. The asymptotic behaviour of Aq 
at infinity will be discussed in a later Section. To solve the equations numerically, it 
is more conveniant to use the field Aq rather than 6; this is why we have expressed 
the boundary condition in terms of that field. On the other hand, the equations take 
a simpler form when written in term of b, so we shall still use it below. 
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Now the volume integral (with the appropriate normalisation of 127?^) of Qq given 
by @ is the degree of the map, or, the Baryon number. It is straightforward to verify 
that when the Ansatz (^) is substituted in and the volume integral is computed 
subject to the boundary conditions given above, the result will equal the integer n 
defined in ( |28D . Thus, the Baryon number of the field configuration (|28D equals the 
vortex number n. In what follows, we will restrict ourselves to unit Baryon number, 
n = 1, i.e.to the Nucleons. 



Before we proceed to substitute the Anstaz (|28|), (p9D into the field equations, we 
calculate the Baryon current (pA]) for the field configurations (p8|) , (|29|) described by 
the solutions we seek. We express the space-like part of this current ji in the radial 
direction flowing out of the normal to the surface of the cylinder which we denote by 
jr, and in the z direction which we denote by jz- The result is 

6 

jr = -[{fgz - gfz)asm'^ f sing (32) 
r 

1 1 

+ 2 (5"^^ - ^9z) sin / cos fsmg + - {af^ - fa^) cos ^f] 

6 

jz = — [{fgr- gfr)asm'^ fsing (33) 
r 

+ 2(5'^^ ~ ^9r) sin / cos / sin 5( + -{dfr - far) cosg] , 

where we have denoted = f etc. Note that the Baryon current (|32D, ([33| ) are not 
sensitive to the charge of the Nucleon, i. e. that the function b{r, z) does not feature 
in them. 

We now turn to the equations to be solved, namely the Euler-Lagrange equations 
arising from the variational principle applied to the Lagrangian (^), in the static limit. 
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Substituting the ansatz ( PB]) ( P^ into these equations of motion leads to 



^ - a,, - asin2(/) sin^l^^) \^ + '-tW + fl + {dl + dl) sin^(/)] 



Ao 



Ao 



A/ 

+4^ 



-(^r + + sin' g) sin / cos /+ 

+9zifrr9z - fz9rr + fr9rz " frz9r + f r / T 9 z))fiVa. f + {frQ^ - Qrfzf COS / 



+ sin (yf/r^ 



+ Pz - '^{91 + sin^ /) cos / sin g 



+ ((a - h ){frr + fzz - fr/r) + 2/^(00^ - 66^) + 2f^{aa^ - bb^)) sin /sin 51 
+2(0' - b'^){fr9r + fz9z)smf cos g)] = 0. 



^5- +'2{frgr + fz9z) cot f 



-f,2 



+4 



A2 



— sm (7cos(7+ 

fz{fz9rr frr9z ~\~ frz9r fr9rz ~\~ fz9r/'^) 
+ fr{fr9zz - fzz9r + frz9z - fz9rz " fz9z/r) 

+ smg/r\{a' - b'){{gl + gl) sin' /-/,'- /|; 
((a' - 6')(fi'r-r + 9zz - 9r/r) + 2gr{aar - 66r) 
+4(a' - 6') (/r-fifr + /^s-^) cos / sin / sin g^ 



cosg 
+ 2gz 
0. 



-66^;)) sin'/ sin 5f 



(34) 

In the case of the Aq = gauge, equations (|34D coincide with the Euler-Lagrange 
equations derived from the positive definite Hamiltonian density (|^). Moreover in 
that case, those equations also coincide with the Euler-Lagrange equations of the 
reduced two dimensional Hamiltonian obtained by subjecting to axial symmetry 
by substituting the Ansatz (p^)-p9|) into it. This is expected due to the strict impo- 
sition of symmetry. In the Aq ^ gauge, the Euler-Lagrange equations are derived 
from the Lagrangian which is not positive definite. Nonetheless these equations 
coincide with those arising from the reduced two dimensional Lagrangian obtained 
by subjecting the Lagrangian (H) to axial symmetry. (This happens also for the Julia 
Zee dyon 0.) 



3.1 ^0 = 0: U{1) Skyrme soliton 

It is easy to see from ( p4D that there are solutions for which 6 = {i.e. Aq = 0). 
As mentioned before, in that case, equation ( P^ can be obtained by minimising the 
Hamiltonian (pOl). Notice also that setting a = is not compatible with our boundary 
conditions {Ai would not be well defined at the origin). We thus expect our gauged 
solution to carry a non-zero magnetic field. 

To show this we have to solve equations (^4]) numerically for the non-vanishing 
functions /(r, z),g{r, z) and a(r, z). 

We have restricted our numerical integrations to the case where the vortex number 
n appearing in the axially symmetric Ansatz (^) is equal to 1, ie.our soliton carries 
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unit Baryon number. 

Using (|30|), we have found numerically that E{1,1,1) = 247r^l.01 whereas for 
those values of Aq, Ai, A2 the lower bound for the energy given by ( p6D is 247r^0.555. 
In Figure 1, we present the total energy for the gauged Skyrmion as a function of 
Ao, together with the lower bound given by (p6l). Note that the asymptotic value 
of -E'(Ao, 1, 1) is 247r^l.232 as Aq — > cxo. As a comparison, the energy (|27D for the 



ungauged Skyrmion is -^^^(l, 1) = 247r^l.232 with a lower bound set at 247r^. We see 
that -^(Ao = 00, 1, 1) = Esk{l, 1) which means that as Aq 00, the gauge coupling 

1 /2 

I/Aq goes to zero and the gauged Skyrmion becomes in this limit the ungauged 
Skyrmion. 

It is interesting to note that the energy of gauged Skyrmion is smaller than the 
energy of the ungauged Skyrmion, as expected, but that on the other hand, the 
amount by which the energy of the gauged Skyrmion exceeds its topological lower 
bound is larger than the excess of the energy of the ungauged Skyrmion above its 
respective topological lower bound. For example we can clearly see from Figure l.a 
that at Ao = 20, the energy of the gauged Skyrmion 1.22 (in units of 247r^) exceeds 
the lower bound 0.95 by 0.27. This is larger than 0.232, the excess of the ungauged 
Skyrmion energy over its lower bound. For smaller values of Aq Figure l.a shows that 
the excess of the energy of the gauged Skyrmion over its lower bound is even larger, 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, z.e.the term proportional to Aq 
in (^), is decreasing as Aq increases. Notice that we could have used for the Maxwell 
Energy the sum of all the terms involving the gauge field functions a and b in (|30D , 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the energy density of the 
Skyrmion in the r, z plane for Aq = 1. One sees clearly that the effect of the gauged 
field is to make the Skyrmion elongated along the z axis. The magnetic field vectors 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we show the configuration 
of magnetic field using arrows to represent the magnetic field vector at each point on 
the grid. Notice that there is a vortex around the point r = 2, ^ = 0. The magnetic 
field is thus generated by a current flowing on a ring centred around the z axis. 

In terms of the usual physical constants^], we have Ag^ = 4e^, Ai = F^/8 and 
A2^^ = 8a^ where we use a instead of the traditional e for the Skyrme coefficient to 
avoid confusion with the electric charge. 

In our units, c = ^ = 1, we have e = {AnaY^'^ where a = 1/137 is the fine 
structure constant. Choosing = 186MeV, we can find the value for a by requiring 
that the energy of the neutron M„ = 939MeV matches the energy of the Skyrmion: 



M„ = 1^^(4,1,1). 



In Figure l.a, we can read the value of E{Xq, 1, 1) (given in units of 247T'^MeV) with 
Ao = 2o?/e^. We now have to find the value of Ao for which 247r^£'(Ao, 1, 1) = 
2(2Ao)"'^^^eM„/F^. The intersection between the curve :(2Ao)^/^ and the curve 
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Figure 1: a) Energy and topological bound of the gauged Skyrmion in units of 247r^. 
b) Ratio of the electromagnetic and the total energy as a function of Aq . 




Figure 2: a) Energy density for the gauged Skyrmion in the (r, z) plane (Aq = Ai = 
A2 = 1). b) Energy density level curve for the gauged Skyrmion (Aq = Ai = A2 = 1). 
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Figure 3: Magnetic Field of the gauged Skyrmion (Aq = Ai = A2 = 1). 

E{\q, 1, 1) in Figure l.a is located in the region where the energy is virtually equal to 
the asymptotic value E{\o, 1, 1) = 1.232. This means that a ?s 37r2l.232F^/M„ ^ 7.2 
and that Aq ~ 1138. We can thus conclude that the effective impact of the Maxwell 
term we have added to the Skyrme model is relatively small. 

This justifies the procedure used in where the Skyrmion was coupled with 
an external magnetic field of a magnetic monopole. Indeed, as the Maxwell field 
generated by a Skyrmion is very small (for the parameters fitting the actual mass of 
the nucleons) the external field is much larger than the Skyrmion's magnetic field. 

It would be interesting to find the differences between the electromagnetic quan- 
tities obtained from the ungauged model, as in 0, and our f/(l) gauged model. We 
are not able to compute the solutions of the U{1) gauged model for the physical value 
of the parameter Aq as this is too large, but as we now know that the influence of 
the gauge field is very small, we can compute the latter perturbatively around the 
(ungauged) Hedgehog as an induced field. This enables the evalution of the energy 
correction and the induced magnetic moment. This perturbative analysis will be be 
carried out in the Section 4. 

It can also be concluded that if the U{1) gauged Skyrme model were quantised as 
in Q (by quantising the zero modes corresponding to the global gauge transforma- 
tion) but taking into account the electromagnetic field generated classically by the 
Skyrmion, the result would not differ very markedly from what was obtained in . 



,U L, i ,1 , |, 1 , 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a function of Aq 
3.2 Aq ^ : charged U{1) Skyrmion 

We can now look for solutions with a non-zero electric charge by requiring that the 
field b in our ansatz (^) does not vanish. To do this we follow the same procedure 
as Julia and Zee ^ and require that the electric field be asymptotically of the form 
^0 = Vo + q/{r'^ + z'^y^'^ where Vq and q are two constants. In practice, one computes 
solutions for different values of Vq and evaluates q by computing the electric flux. We 
sought only those solutions, for which the electric flux equals 47r times the charge of 
the electron. 

It is important to realise that in this case, equations (|3^) are obtained after min- 
imising the action and thus they do not minimise the Hamiltonian (^0|). 

In our units, the charge of the electron is 0.303. In Figure 4 we show the energy 
as a function of Aq, as well as Vq as a function of Aq, so that q = 0.303. 

One can see that, for a fixed value of Aq, the energy of the charged gauged Skyrmion 
is smaller than the energy of the ungauged Skyrmion when Aq < 7 but it is always 
larger than the energy of the uncharged gauged Skyrmion. If the energies of the 
electrically charged and uncharged gauged Skyrmions were interpreted as the the 
masses of the Proton and the Neutron mp and mA?, then on this purely classical level 
we would have to conclude that {mp — rri]^) > which is not correct. This is expected 
on the basis of its analogy with the dyon . Clearly, to calculate this mass difference 
correctly one would have to perform the collective coordinate quantisation as in Ref. 
P], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is unfortunately very 
difficult to carry out the numerical computations accurately when Aq is very large. 

At this stage, it is worth saying a few words about the numerical methods we 
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have used. To compute the static solutions, we have employed a relaxation method 
using finite differences on a regular grid {dr = dz). This discretisation method is 
similar to the one employed in the numerical computation of the solutions of Skyrme 
models in 2+1 and 3+1 dimensions ||16|, |1^ 0. To compute the electrically charged 
solutions, we have imposed the boundary condition 6(oo, oo) = Vq for different Vq 
and using a dichotomic method, we have determined the values that give a solution 
with the same electric flux as the proton. Most of the simulations where done on 
200 X 400 or 300 x 600 grids. By computing the same solution for various lattice sizes, 
we have empirically obtained the following relation for the expression of the energy 
of a solution : E = Eq + Kdr"^ where Eq is the exact solution, dr = dz the lattice 
spacing and K is a constant which depends on Aq but takes values between 1 and 
0.5. We see thus that the energies we have obtained are accurate to within one or one 
half of a percent. This inaccuracy in the value of the energy is comparable to that 



of many other similar works on 2 dimensional systems ]1^ |T8[, and though it might 



look large, it does not affect any of the conclusions we have drawn. 



4 Perturbation around the Hedgehog 

We have seen from the work of 3.1 above that the energy of the U{1) gauged Skyrmion 
for the physical values of the parameters, namely of the pion decay constant and the 
U{1) coupling, does not differ signiflcantly from that of the ungauged Hedgehog. It 
is therefore justifled, for these values of the parameters, to treat the U{1) fleld as a 
perturbation to the Hedgehog in the same way as Klinhamer and Manton Jl^ treat 
the U{1) fleld as a pertubation to the SU{2) sphaleron. We will then be able to 
compute the magnetic moment of the Neutron, as well as the (small) deviation of its 
mass from that of the Hedgehog. 

The equations for the flelds are derived from the Lagrangian (|l]). The 

equation for will be of the form 

Ao d^F,,^ = . (35) 

The method consists of setting the gauge fleld to zero in the current in (|35D (and 
in the equation for 0") and calculating the resulting induced electromagnetic fleld A^. 
The gauge fleld computed this way can then be interpeted as the U{1) fleld generated 
by the (ungauged) Hedgehog Skyrme fleld. With this perturbative procedure, it is 
possible to calculate the induced static magnetic potential Aj (i = 1, 2, 3), but not the 
static electric potential Aq, which in this scheme vanishes and can only be calculated 
non-perturbatively. The reason simply is that restricting to the use of the static 
Hedgehog, the zeroth component of the current jo at = vanishes, resulting in 
turn in vanishing induced potential Aq according to (|35|). 

As a consequence the electric fleld will be identicaly zero, which implies that we 
can derive the equation from the static Hamiltonian rather than from the Lagrangian. 
Notice also that we could try to compute perturbatively a solution for the electrically 
charges skyrmion by keeping in jo the terms proportional to Aq, instead of setting 
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Aq to 0, and impose the condition that the electric field is assymptoticaly like that 
of the proton. This pertubation method would not make much sense though as one 
would expect the electric field to be quite large close to the skyrmion. 

The relevant energy functional is (^), and the resulting equation arising from the 
variation of the gauge field Ai is 

Ao djFij = ji (36) 



We are concerned here with the case where Ai = {p7[j and the chiral field 0" = 
(0", 0^, 0^) in (P^) describes the Hedgehog, i.e. 

0° = sinF(/2)x" , (f)^ = sin F{R)x\ (f)^ = cos F (R) , (38) 



where now a;" = x"'/R, with R = \fr^+z^. By virtue of equation (R3) the current 



(|37|), given by ( pH]) and Ai = 0, will now induce a (small) U (1) field Aj, with curvature 
Fij = diAj — djAi. 

The shift in the energy of the Hedgehog due to the induced U{1) field Aj is 

= 1 d^x {XoFijFij + 4A, j,) , (39) 



in which jj = ji(0) is the current ( pTf ) for Ai = 0. (Note that all quantities evaluated at 
Ai = are denoted by Roman script, e.g. jj = ji(0), as well as the induced connection 
and curvature Aj and Fj^.) When equation (^) is satified for the induced U{1) field, 

AE = -Ao J d?x YijYij (40) 

= 2 j d'xk^U. (41) 

which is, as expected, a strictly negative quantity. 



The current = (ja, js) in (|39D and (^T]) is given, for the Hedghog field configura- 
tion (11), by 

sin^F /Ai , sm^F\ \ . 

ja = 0. (43) 

We now note that di jj = 0, which means that equation (p6D, for Ai = 0, takes the 
following form in terms of the induced U{1) connection Aj = (A^, 0) 

AoA A„ = -j„ . (44) 

The solution is well known and can be expressed, using the obvious notation ja(x) = 
j(-R) SapXp in terms of (|4^), as 

1 /■ 1 

A,(x) = ---^Sc^p / I ]{R!) xp d-K' , (45) 

47rAo J |x — x'l 
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with ^ ^ 

n -EE 5J^#T mo'A') Yi'HoA) . 

After performing the angular integrations we have 

A„(x) = -I{R) e^pxp , (46) 
with I{^i){R) given by the integral 

^(«' - k (C § J<«'' + r ^<«'' «'^^«') ■ <^^> 

Finally, in the ^ 1 region of interest, the induced U{1) potential is 

A„(x) = -— eaf3 Xfs , with , / = — / ds , (48) 

ii^ oAo "'0 

to be evaluated numerically using the numerically constructed hedgehog profile func- 
tion F{x) (|3|). 

Comparing (|48|) with the usual Maxwell potential of a magnetic dipole ^ 

A(x) = ^ , 

we find that fj, = (0, 0, fj,) is 

/i = 47r/ . (49) 

We can evaluate the magnetic moment (^) and the energy correction (^I]) induced 
by the electromagnetic field by evaluating the integral (^7[) and (^8|) numerically. If 
we take the experimental values = 186MeV and a = 7.2 we obtain 

/i = 0.01393/m = 0.43nm (50) 
AE = -O.lkeV. (51) 

The experimental value for the magnetic moment of the proton and the neutron are 
respectivly fip = 0.0902/m = 2.79nm and /i„ = 0.0617/m = — 1.91nm. If on the 
other hand we take the values of the parameters derived in [0], = 129MeV and 
a = 5.45, we have 

/i = 0.0468/m = 1.449nm (52) 
AE = -0.32keV. (53) 

The magnetic moment of a particle is strictly speaking a quantum property and it 
should be computed by quantising the SU{2) gauge degree of freedom as in |]^. Nev- 
ertheless, we see that if the take the parameters derived in ||^ the classical magnetic 
moment is of the correct order of magnitude. The sign is of course undetermined as 
the classical magnetic moment is a vector. We can thus conclude that our model offers 
a reasonable classical description of Nucleons and affords a method for computing the 
electromagnetic field generated by the Skyrmion, classically. It is quite surprising to 
see that a quantum property like the magnetic moment can be reasonably predicted 
by a purely classical procedure. 
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5 Summary and discussion 



We have shown that the SU(2) Skyrme model gauged with U{1) has two types of 
finite energy static solutions, electrically uncharged and charged respectively. Both 
of these solutions are axially symmetric and carry no magnetic charge but support a 
magnetic field shaped like a torus centred around the axis of symmetry, albeit resulting 
in zero magnetic fiux. The uncharged solutions, like the ungauged Skyrmion, have a 
topological lower bound. The electrically charged solutions are the analogues of the 
Julia- Zee dyons of the Georgi-Glashow model. 

Concerning the stability of the electrically neutral solution, which is expected to be 
stable by virtue of the lower bound on the energy, we have not made any quantitative 
effort to test it. We expect however that the solitons of this gauged Skyrme model are 
stable, or that at least they have stable branches for all values of the parameters in the 
model. This expectation is based on our knowledge of the corresponding situation 
when the Skyrme model is gauged instead with S0{3) ||2^, |2l|], in which case the 
equations arising from the imposition of spherical symmetry were one dimensional, 
and hence technically much more amenable to the numerical integration. In that 
case it was found that in addition to stable branches of solutions, there were also 
some unstable branches bifurcating from the former, the important matter being 
that there were indeed stable branches of solutions, characterised by the (ranges of 
the) parameters of the model. It would be very interesting to carry out the analysis 
corresponding to that of |20, 21], for the considerably more complex case of the axially 
symmetric equations at hand. This however is technically beyond the scope of the 
present work. 

The energies of the gauged uncharged Skyrmions are smaller than the energy of the 

1 /2 

usual ungauged Skyrmion. When the gauge coupling I/Aq goes to 0, the uncharged 
gauged Skyrmion tends to the ungauged Skyrmion. We also note that the energy of 
the electrically charged Skyrmion is higher than the uncharged one, just as the mass 
of the dyon is higher than that of the monopole of the Georgi-Glashow model. 

Perhaps the most interesting physical result of the present work is that when 
parameters in the model are fitted to reproduce physical quantities, it turns out that 
the effect of the Maxwell term in the Skyrme Lagrangian is very small. This is because 
for the physical value of the constant Aq = 1138, the energy of the gauged uncharged 
Skyrmion differs little from that of the ungauged Skyrmion, as seen from Figure lb. 
The gauged Skyrmion field itself is thus nearly radially sjTiimetric (though the gauge 
field is not). 

Having found that the infiuence of the electromagnetic field on the Skyrmion 
is small, we were pointed in the direction of treating the magnetic potential as an 
induced field perturbatively around the (ungauged) Hedgehog. We have been able 
thus, to compute the classical magnetic moment of the (uncharged) Skyrmion of unit 
Baryon charge, namely of the Neutron. The result is that the classical magnetic 
moment of the Skyrmion matches surprisingly well to the experimental values of the 
magnetic moments of the Nucleons. 
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Abstract 

We use a prescription to gauge the SU{2) Skyrme model with a U{1) field, 
characterised by a conserved Baryonic current. This model reverts to the usual 
Skyrme model in the limit of the gauge coupling constant vanishing. We show 
that there exist axially symmetric static solutions with zero magnetic charge, 
which can be electrically either charged or uncharged. The energies of the 
(uncharged) gauged Skyrmions are less than the energy of the (usual) ungauged 
Skyrmion. For physical values of the parameters the impact of the U{1) field 
is very small, so that it can be treated as a perturbation to the (ungauged) 
spherically symmetric Hedgehog. This allows the perturbative calculation of 
the magnetic moment. 
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1 Introduction 



For a long time now, much attention has been paid to the Skyrme [1] model in 3 
dimensions. It is believed to be an effective theory for nucleons in the large N limit 
of QCD at low energies. The classical properties as well as the quantum properties 
of the model are in relatively good agreement with the observed properties of small 
nuclei [2, 3, 4]. 

Gauged Skyrme models have been used in the past. The U{1) gauged model [2, 5] 
was used to study the decay of nucleons in the vicinity of a monopole [5], while the 
SU{2)l gauged model [7] was used to study the decay of nucleons when the Skyrme 
model is coupled to the weak interactions [7]. The Skyrme model has also been used 
to compute the quantum properties of the Skyrmion [3] where the gauge degrees of 
freedom are quantised to compute the low energy eigenstates of a Skyrmion. These 
states were identified as the proton, the neutron and the delta. 

The aim of this work is to show that the Skyrme model can be coupled to a self 
contained electromagnetic field and that this U{1) gauged model has stable classical 
solutions. In addition to these solitons with vanishing magnetic and electric flux, we 
show that this system supports solutions with nonvanishing electric flux which are 
analogous to the dyon solutions of the Georgi-Glashow model, just as the uncharged 
solitons are the analogues of the monoploles [8] of that model. The electrically charged 
lumps have larger energy, or mass, than the uncharged soliton, just like the ,Iulia-Zee 
dyon [9] has larger energy, or is heavier, than the (electrically uncharged) monopole. 
We shall refer to these lumps as charged U{1) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Maxwell gauged Skyrme 
model, the present work is also an example of a soliton in a d-dimensional SO{N) 
gauged 5''^ model with N < d for the case d = 3, N = 2 , extending the results of 
Ref. [10] which were restricted to the N = d cases. (The work of Ref. [10] consists 
of establishing topological lower bounds for the generic case, encompassing earlier 
examples in two [11] and three [12, 13] dimensions respectively.) The gauging pre- 
scription used here by us coincides precisely with that used in Ref. [5] and permits 
the establishing of a topological lower bound which did not feature in Ref. [5] and 
which is carried out here to establish the stability of the soliton. Such lower bounds 
are absent in the other prescription of gauging the Skyrme model as in Refs. [7]. 
(Notice that we name the sigma models after the manifold in which the fields take 
their values rather than using the name of the symmetry group for the model. Thus 
what is sometime called the 0{d+ 1) model in the literature will be refered to as the 
5^ model.) 

The U{1) gauged SU{2) Skyrme model is described by the Lagrangian [5] 

C = ^Tr{D,UD,U^) - ^Jr{[iD,U)U\ iD,U)U^y - \tI^ 

where the C/(l) gauge covariant derivative is 

D^V = d^V ^ieA^\Q,lJ\, (1) 
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5 Summary and discussion 



We have shown that the SU{2) Skyrme model gauged with U{1) has two types of 
finite energy static solutions, electrically uncharged and charged respectively. Both 
of these solutions are axially symmetric and carry no magnetic charge but support a 
magnetic field shaped like a torus centred around the axis of symmetry, albeit resulting 
in zero magnetic flux. The uncharged solutions, like the ungauged Skyrmion, have a 
topological lower bound. The electrically charged solutions are the analogues of the 
Julia-Zee dyons [9] of the Georgi-Glashow model. 

Concerning the stability of the electrically neutral solution, which is expected to be 
stable by virtue of the lower bound on the energy, we have not made any quantitative 
effort to test it. We expect however that the solitons of this gauged Skyrme model are 
stable, or that at least they have stable branches for all values of the parameters in the 
model. This expectation is based on our knowledge of the corresponding situation 
when the Skyrme model is gauged instead with S0{3) [20, 21], in which case the 
equations arising from the imposition of spherical symmetry were one dimensional, 
and hence technically much more amenable to the numerical integration. In that 
case it was found that in addition to stable branches of solutions, there were also 
some unstable branches bifurcating from the former, the important matter being 
that there were indeed stable branches of solutions, characterised by the (ranges of 
the) parameters of the model. It would be very interesting to carry out the analysis 
corresponding to that of [20, 21], for the considerably more complex case of the axially 
symmetric equations at hand. This however is technically beyond the scope of the 
present work. 

The energies of the gauged uncharged Skyrmions are smaller than the energy of the 

1 /2 

usual ungauged Skyrmion. When the gauge coupling l/Ag goes to 0, the uncharged 
gauged Skyrmion tends to the ungauged Skyrmion. We also note that the energy of 
the electrically charged Skyrmion is higher than the uncharged one, just as the mass 
of the dyon is higher than that of the monopole of the Georgi-Glashow model. 

Perhaps the most interesting physical result of the present work is that when 
parameters in the model are fitted to reproduce physical quantities, it turns out that 
the effect of the Maxwell term in the Skyrme Lagrangian is very small. This is because 
for the physical value of the constant Aq = 1138, the energy of the gauged uncharged 
Skyrmion differs little from that of the ungauged Skyrmion, as seen from Figure lb. 
The gauged Skyrmion field itself is thus nearly radially symmetric (though the gauge 
field is not). 

Having found that the influence of the electromagnetic field on the Skyrmion 
is small, we were pointed in the direction of treating the magnetic potential as an 
induced field perturbatively around the (ungauged) Hedgehog. We have been able 
thus, to compute the classical magnetic moment of the (uncharged) Skyrmion of unit 
Baryon charge, namely of the Neutron. The result is that the classical magnetic 
moment of the Skyrmion matches surprisingly well to the experimental values of the 
magnetic moments of the Nucleons. 
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where T^iy = df,Ai, — d^An and where the charge matrix of the quarks is expressed 
as Q = + 0-3). This differs from the covariant derivative of Ref. [5] only in the 
unimportant matter of the sign of i in (1), which we have chosen for consistency of 
the convention used in Ref. [10]. 

In what follows it will be more convenient [10] to parametrise the Skyrme field as 
an valued field = (j)^), a = 1,2, A = 3,4 subject to the constraint |</!)"p = 1. 
The two fields U and (J) are related to each other via the following expression 

U = </)V", = U^ = </)"f" (2) 

where = {ia°',ia^,l) and = {—ia°',—ia^,l), in terms of the Pauli matrices 

The gauge covariant derivative now can be re-expressed as 

L>>" = a>" + A^£"^</.^, D,(t>^ = d^(t>^. (3) 

where A^ = eAf, and F^^ = eT^„. 

The Lagrangian for the U(l) gauged Skyrme model can then be written as 

£ = -XoFl^ + X^D^r? - A2|/^[>"I^.]</'T (4) 

where the square brackets on the indices imply (total) antisymmetrisation and where 
Xq^ = 4e^, Ai = F^/S and A2 ^ = 8a^. The late Greek indices n label the Minkowskian 
coordinates, while the early Greek indices a = 1,2 and the upper case Latin indices 
A = 3, A label the fields = cj)^). 

The static Hamiltonian pertaining to the Lagrangian (4) is 

n = Aoi^2 + All A</i1' + A2|L'[,rD,]</.''|2 

(5) 

+i2x,\dM' + ^^{Aiir p + i6X2[\rmcf>^\' + iiw^m}, 

where the indices i = a, 3 label the space-like coordinates. 

To find the static solutions of the model, one would usualy solve the Euler La- 
grange equations which minimise the Hamiltonian (5), but because of the electric 
potential Aq, one must solve the Euler Lagrange equations derived from the La- 
grangian (4). We then look for static solutions, but, as for the Julia-Zee dyon [9], we 
have to impose the proper asymptotic behaviour for the electric potential to obtain 
static solutions which are electrically charged (in the classical sense, i.e. solutions 
where the flux of the electric held is non zero). 

When the full equations of motion are written down, one flnds as expected that 
there are static solutions for which Aq = 0, i.e. solutions for which the electric field is 
identically zero. For these solutions in the temporal gauge, the equations of motion 
reduce to the equations obtained by minimising the Hamiltonian (5). We study the 
solutions of unit Baryon charge of the U{1) gauged Skyrme model with and whithout 
an electric field, for various values of the U{1) coupling constant (or equivalently the 
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Skyrme coupling). For physical values of these parameters in the model, we find that 
the energy (mass) of the gauged Skyrmion does not differ significantly from that of 
the ungauged charged-1 Skyrmion, namely the familiar hedgehog [1]. This implies 
that for these values of the physical parameters, the U{1) gauged Skyrmion can be 
regarded as a perturbation of the (ungauged) hedgehog, enabling the computation 
of the magnetic moments of the gauged Skyrmion (i.e. the Neutron) and the shift 
of the energy of the gauged Skyrmion away from the energy of the hedgehog [1]. 
perturbatively using the method employed by Klinkhamer and Manton [14] for the 
sphaleron of the Weinberg-Salam model. 

In Section 2, we define the topological charge and establish the corresponding 
lower bound on the energy functional. In Section 3 we present the solutions which 
have no electric fields in the first subsection and electrically charged solutions in 
the second subsection. The perturbation analysis of the gauged Skyrmion around 
the (ungauged) hedgehog is carried out in Section 4, and Section 5 is devoted to a 
discussion of our results. 

2 The topological charge and lower bound 

The definition of the topological charge is based on the criterion that it be equal to the 
Baryon number, namely the degree of the map. For the gauged theory however, this 
quantity must be gauge invariant as well. This requirement can be systematically [10] 
satisfied by arranging the gauge invariant topological charge density to be the sum 
of the usual, gauge variant winding number density 



plus a total divergence whose surface integral vanishes due to the finite energy condi- 
tions, such that the combined density is gauge invariant. In 3 dimensions, this is given 
explicitly in Refs. [10, 13] for the 50(3) gauged model, and for the present case of 
interest, namely the S0{2) gauged model, the charge density can be derived from 
that of the 50(3) gauged model by contraction of the gauge group 50(3) down to 
50(2). It can also be arrived at directly. To state the definition of the charge, we 
denote the gauge covariant counterpart of (6) by 



so that using the notations (6) and (7) we have the definition of the gauge invariant 
topological charge 



^0 = s,jk£'""'%rdj(f>''dk<f>'<f>'' , 



(6) 



(7) 



g = go + difli 



'1 1 



(8) 
(9) 



In (8) the density is the following gauge variant form 



(10) 
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with [19] 

1^ -^1 i=0 m=-t '^'^ ^ ^ ^> 

After performing the angular integrations we have 

A„(x) = -I{R) SafiXfi , (46) 
with I{i){R) given by the integral 

Finally, in the i? 3> 1 region of interest, the induced U{1) potential is 

/ -^1 
A„(x) = -— Xfi , with , I = —J^ s^j{s) ds , (48) 

to be evaluated numerically using the numerically constructed hedgehog profile func- 
tion F{x) (38). 

Comparing (48) with the usual Maxwell potential of a magnetic dipole /i 
we find that /i = (0, 0, //) is 

= 47r7 . (49) 

We can evaluate the magnetic moment (49) and the energy correction (41) induced 
by the electromagnetic field by evaluating the integral (47) and (48) numerically. If 
we take the experimental values = 186MeV and a = 7.2 we obtain 

fi = 0.01393/m = 0.43nm (50) 
AE = -O.lkeV. (51) 

The experimental value for the magnetic moment of the proton and the neutron are 
respectivly fip = 0.0902/m = 2.79nm and = 0.0617/m = —l.dlnm. If on the 
other hand we take the values of the parameters derived in [3], F^^ = 129MeV and 
a = 5.45, we have 

Id = 0.0468/m = l.U9nm (52) 
AE = -0.32keV. (53) 

The magnetic moment of a particle is strictly speaking a quantum property and it 
should be computed by quantising the SU{2) gauge degree of freedom as in [3]. Nev- 
ertheless, we see that if the take the parameters derived in [3] the classical magnetic 
moment is of the correct order of magnitude. The sign is of course undetermined as 
the classical magnetic moment is a vector. We can thus conclude that our model offers 
a reasonable classical description of Nucleons and affords a method for computing the 
electromagnetic field generated by the Skyrmion, classically. It is quite surprising to 
see that a quantum property like the magnetic moment can be reasonably predicted 
by a purely classical procedure. 



17 



^0 to 0, and impose the condition that the electric field is assymptoticaly like that 
of the proton. This pertubation method would not make much sense though as one 
would expect the electric field to be quite large close to the skyrmion. 

The relevant energy functional is (23), and the resulting equation arising from the 
variation of the gauge field Ai is 

Ao djFij = ji (36) 
Ji = (lAi^'^Ar + A2 \{d,\<t>^\'')D,rD^<P^ + 2<t>^D^irD,^<t>^D^ct>^]) . (37) 

We are concerned here with the case where = (37) and the chiral field = 
(</)", f/)^, f/)^) in (37) describes the Hedgehog, i.e. 

= sin F{R) , </)^ = sin F{R) x^ = cos F{R) , (38) 

where now x"" = x"" /R, with R = \/r^ + . By virtue of equation (36) the current 
(37), given by (38) and Ai = 0, will now induce a (small) U{1) field Aj, with curvature 
Fjj = diAj — djAi. 

The shift in the energy of the Hedgehog due to the induced U{1) field Aj is 

AE = J Sx {X^YijYij + Aki ]i) , (39) 

in which jj = ji(0) is the current (37) for Ai = 0. (Note that all quantities evaluated at 
Ai = are denoted by Roman script, e.g. jj = ji{0), as well as the induced connection 
and curvature Aj and Fjj.) When equation (36) is satified for the induced ^7(1) field, 

AE = -Xo J d^x YijYij (40) 

= 2 j d^xAJ,, (41) 

which is, as expected, a strictly negative quantity. 

The current jj = (ja, ja) in (39) and (41) is given, for the Hedghog field configura- 
tion (38), by 

sin^F /Ai , / ,2 sin^FA A ^ 
= -^(y + 2A2(^F'2 + -^j (42) 

j3 = 0. (43) 

We now note that di jj = 0, which means that equation (36), for Ai = 0, takes the 
following form in terms of the induced U{1) connection Aj = (Aq,,0) 

AoA A„ = -u . (44) 

The solution is well known and can be expressed, using the obvious notation ja,(x) = 
i{R) Eapxp in terms of (42), as 

1 /■ 1 

A„(x) = ----£«/3 / jT o!x' , (45) 

47rAo J |x — X I 
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The flux of Q,j vanishes, as can deduced by anticipating the finite energy conditions 
to be stated later. 

Note that the 3-volume integral of Qq in (8) is the degree of the map for the 
ungauged system namely the Baryon number. 

Identifying g, (9), with the naught component of the Baryon current, j'^ is 
defined by 

f = e^''''''eabcddA''dp(l>'da(l>''4>'' ^^e^'''''' dMpe^^^'^^'d,^^^^^ (11) 
= e^'""' EabcdD.r D,^'D„^^'' - ^ e^'^P" F,,{eABct>'' D„ci>^)- (12) 

The 4-divergence of (11) receives a contribution only from its first term, which 
being locally a total divergence implies that the 3-volume integral of j° is a conserved 
quantity. Alternatively we consider the 4-divergence of (12), 

d^f = 6e''-'P-eo^^eABD,(j)''D,<l>^D,<j)^D,(j)'' (13) 

which is analogous to the corresponding quantity in the work of Goldstone and 
Wilczek [6]. This contrasts with the expression for the total divergence of the topolog- 
ical current in the work of d'Hoker and Farhi [7], where a different gauging prescription 
is used leading to that quantity being equal to the local anomaly. 

We now proceed to find a model whose Hamiltonian Hq is bounded from below by 
the topological charge density defined by (9). We will then show that the Hamiltonian 
(5) is given by Tio plus certain positive definite terms. 

First of all, we reproduce the density Qg, (7), in (9) by using the following in- 
equality 

{KsD^r - Sijke'"'"'4Dj(l)''Dk(l>'(l>''Y > (14) 

where the two constants and K2 have the dimensions of length. Expanding the 
square, we get go on the right hand side of 

4{Dirf + 4{D[irDj](f^'f > 2i^s4qg. (15) 

To reproduce the other term in (9), ^SijkFij{s^^ (j)^ dkcj)^) , we use the following 
inequality 

{KlFij - ^K,eijke^''(l)''Dk(l>y > (16) 

yielding 

+ 4^{s^^4>''D,4>^f > KlK,8,j,F,^{8^'^4>''D,4>^)- (17) 

With the special choice for the relative values of the constants Sk^k^ = the 
sum of (13) and (15) yields the following 

+ 4{D.rf + 4{D^^rD,^c|>'f + '^-^{e^^cl>^D,cl>y > 2K,4g. (18) 
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The right hand side of (18) is now proportional to the topological charge density 
Q defined by (9) so that the inequality (18) can be interpreted as the topological 
inequality giving the lower bound on the energy density functional if we define the 
latter to be the left hand side of (18), namely 

^0 = + ^liD^rf + 4{D^^rD,^<t>'f + ^-^{e^^<t>^D4^f. (19) 







The Hamiltonian system (19) is almost the Hamiltonian of the gauged Skyrme 
model (5) (remember that = 0). It differs from the latter only in its last term. 
Now we can use the identity 



(20) 



and add the positive definite term 



A]\2 



appearing on the 



2 " 

right hand side of (20) to 1-Lq in (19) to end up with the Hamiltonian for the U{1) 
gauged Skyrme model: 

n = 4Ff^ + KliD^rf + KliDy^'^D^^ct,')^ > 2^i,4g (21) 

which is nothing but the static Hamiltonian (5) in the temporal gauge Aq = 0, and 
where 

Ai = «:^(l + ^), ^0 = 4, ^2 = 4- (22) 

By virtue of (18), (21) is also bounded from below by 2KsKlg, namely by a number 
proportional to the topological charge density g. 

We thus see that Ho can be considered as a minimal {U{1) gauged) model, but 
from now on, we will restrict our attention to the physically more relevant model (21) 
and integrate it mmierically to find its topologically stable finite energy solitons. 

The soliton solutions to the system (21) can only be found by solving the second- 
order Euler-Lagrange equations, and not some first-order Bogomol'nyi equations since 
saturating the inequalities (14) and (16) would not saturate the lower bound on the 
energy density functional H. In this context we note that saturating (14) and (16) 
does indeed saturate the topological lower bound on the functional T-Lq by virtue of 
the inequality (18), and should it have turned out that the Bogomol'nyi equations 
arising from the saturation of (14) and (16) supported non-trivial solutions, then T-Lq 
would have been a very interesting system to consider. As it turns out however, these 
Bogomol'nyi equations have only trivial solutions in exactly the same way as in the 
case of the (ungauged) Skyrme model [1]. 

The energy for the static configuration, when the electric field vanishes, is ex- 
pressed as 

^(Ao, Ai, A2) = I d'x [Aoi^5 + Ai(A</'")' + \2{Dii(l>''Dj^(l>'f] (23) 
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have used. To compute the static solutions, we have employed a relaxation method 
using finite differences on a regular grid (dr = dz). This discretisation method is 
similar to the one employed in the numerical computation of the solutions of Skyrme 
models in 2+1 and 3+1 dimensions [16, 17, 18]. To compute the electrically charged 
solutions, we have imposed the boundary condition 6(oo, oo) = Vq for different Vq 
and using a dichotomic method, we have determined the values that give a solution 
with the same electric flux as the proton. Most of the simulations where done on 
200 X 400 or 300 x 600 grids. By computing the same solution for various lattice sizes, 
we have empirically obtained the following relation for the expression of the energy 
of a solution : E = Eq + Kdr^ where Eq is the exact solution, dr = dz the lattice 
spacing and K is a constant which depends on Aq but takes values between 1 and 
0.5. We see thus that the energies we have obtained are accurate to within one or one 
half of a percent. This inaccuracy in the value of the energy is comparable to that 
of many other similar works on 2 dimensional systems [17, 18], and though it might 
look large, it does not afl^ect any of the conclusions we have drawn. 

4 Perturbation around the Hedgehog 

We have seen from the work of 3.1 above that the energy of the U{1) gauged Skyrmion 
for the physical values of the parameters, namely of the pion decay constant and the 
U{1) coupling, does not differ significantly from that of the ungauged Hedgehog. It 
is therefore justified, for these values of the parameters, to treat the U{1) field as a 
perturbation to the Hedgehog in the same way as Klinhamer and Manton [14] treat 
the U{1) field as a pertubation to the SU{2) sphaleron. We will then be able to 
compute the magnetic moment of the Neutron, as well as the (small) deviation of its 
mass from that of the Hedgehog. 

The equations for the fields (0°, .4^) are derived from the Lagrangian (1). The 
equation for ^4^ will be of the form 

Ao d^Fi,^ = . (35) 

The method consists of setting the gauge field .4^ to zero in the current in (35) (and 
in the equation for (/)") and calculating the resulting induced electromagnetic field A^. 
The gauge field computed this way can then be interpeted as the U{1) field generated 
by the (ungauged) Hedgehog Skyrme field. With this perturbative procedure, it is 
possible to calculate the induced static magnetic potential A,j (i = 1, 2, 3), but not the 
static electric potential Aq, which in this scheme vanishes and can only be calculated 
non-perturbatively. The reason simply is that restricting to the use of the static 
Hedgehog, the zeroth component of the current jo at = vanishes, resulting in 
turn in vanishing induced potential Aq according to (35). 

As a consequence the electric field will be identicaly zero, which implies that we 
can derive the equation from the static Hamiltonian rather than from the Lagrangian. 
Notice also that we could try to compute perturbatively a solution for the electrically 
charges skyrmion by keeping in jo the terms proportional to Aq, instead of setting 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a function of Aq 

3.2 Aq ^ : charged C/(l) Skyrmion 

We can now look for solutions with a non-zero electric charge by requiring that the 
field b in our ansatz (29) does not vanish. To do this we follow the same procedure 
as Julia and Zee [9] and require that the electric field be asymptotically of the form 

^0 = Vo + q/{r'^ + z'^y/'^ where Vq and q are two constants. In practice, one computes 
solutions for different values of Vq and evaluates q by computing the electric flux. We 
sought only those solutions, for which the electric flux equals in times the charge of 
the electron. 

It is important to realise that in this case, equations (34) are obtained after min- 
imising the action and thus they do not minimise the Hamiltonian (30). 

In our units, the charge of the electron is 0.303. In Figure 4 we show the energy 
as a function of Aq, as well as V'o as a function of Aq, so that q = 0.303. 

One can see that, for a fixed value of Aq, the energy of the charged gauged Skyrmion 
is smaller than the energy of the ungauged Skyrmion when Aq < 7 but it is always 
larger than the energy of the uncharged gauged Skyrmion. If the energies of the 
electrically charged and uncharged gauged Skyrmions were interpreted as the the 
masses of the Proton and the Neutron mp and ruff, then on this purely classical level 
we would have to conclude that [mp — mjsi) > which is not correct. This is expected 
on the basis of its analogy with the dyon [9]. Clearly, to calculate this mass difference 
correctly one would have to perform the collective coordinate quantisation as in Ref. 
[3], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is unfortunately very 
difficult to carry out the numerical computations accurately when Aq is very large. 

At this stage, it is worth saying a few words about the numerical methods we 
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and performing the dilation x ax, ^ a ^A^, we get 

^(Ao, Ai, A2) = / d'x [^Ff^ + a\,{D,rf + ^{D^irDi^<t>'f]. (24) 
If we choose a = {^)^^'^ then we have 

^(Ao,Ai,A2) = (AiA2)^/^^(^,l,l), (25) 

from which we see that we can set Ai = A2 = 1 without any loss of generality. By 
virtue of (21) and (25), we can finally state 




Notice that for the usual Skyrme model we have 

^.;t(Ai,A2) = Sdx^[\Mrf + M{diirdj](l>'f] 

= (AiA2)i/^^.,(l,l) (27) 
> 2(AiA2)i/2j-^^3^^_ 

We will use (27) to compare the numerical solutions of the gauged Skyrme model 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability considerations discussed 
in this Section apply only to the solutions with no electric field, z.e.with Aq = 0. 



3 The soliton and the charged U{1) Skyrmion 

To find the static solutions, we have to look for the largest symmetry group of the 
functional to be subjected to the variational principle, and look for solutions which 
are invariant under that symmetry group. For the solutions in the Aq = gauge 
this is the static Hamiltonian (21), while for the solutions in the / it is the 
Lagrangian (4). For our choice of gauge group the largest symmetry is the 50(2) 
group corresponding to an axial rotation in space-time and a gauge transformation 
on the gauge field. Defining the axial variables r = \Jxi+ X2 and and z = x^'m terms 
of the coordinates Xi = {x^, X3), a = 1,2, the most general axially symmetric Ansatz 
[15] for the fields (f)" = ((/i", c/i^) (with a = 1, 2 and A = 3, 4), and, Ai = (A«, A3), is 

sin / cos g , 4>'^ = cos / , (28) 

Ciir, z) , b(r,z) , , 

A3 = Ao = (29) 



sin / sin g , 



A^ 



a{r, z) — n 



C2{r, z) 



-Xr 
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with n° = {sin n(f), cos n(f)) in terms of the azimuthal angle (f) and Xa = Xa/r. n in n" 
is the vorticity, which for the Nucleons of interest to us here, equals one, n = 1. The 
functions a, b,Ci,C2, f and g both depend on r and z. 

Our Ansatz (29) for the U{1) field consists of decomposing the latter in the most 
general tensor basis possible. We will find out below, when we compute the Euler- 
Lagrange equations, that the functions ci and C2 vanish identically. Anticipating this, 
we suppress them henceforth. In its final form this Ansatz agrees with that used in 
Ref. [15], the latter being arrived at by specialising the Rebbi-Rossi Ansatz for the 
axially symmetric S0{3) field. 

The static Hamiltonian, i e.the Too component of the energy mometum tensor T^^,, 
is then given by 



H 



~ 9 



b\2 
r ) 



/2 + /2 + sin^ + + «!_t6! sin2 / sin2 g 



+2X2 sin^ / 



ifrgz 



9[^{f^ + Pz + + 9l) sin' /)]] }rdrdz. 

(30) 

The boundary conditions for the Skyrmion fields are the same as the boundary 
conditions for the hedgehog ansatz when expressed in the cylindrical coordinates 
where g = 7r/2 + arctan(2;/r) and defining R= {z^ -[- r'^Y^'^, f = f{R) with /(O) = tt 
and lim^_^oo f{R) = 0. From this we can deduce that the function / has a fixed value 
at the origin and at infinity. For smothness along the z axis, each field, that is /, g and 
Aa must satisfy the condition that the partial derivative with respect to of the field at 
r = vanishes. The boundary conditions and the asymptotic behaviours for a and b 
are chosen so that the gauge fields are well defined and Aq looks asymptoticaly like 
a coulomb field (i.e. with an electric charge but no magnetic charge). We also require 
that the total energy be finite. These conditions leads to the following constraints: 



f{r ^ 00, z^ 00) =0 fr{r = 0,z) =0 

g{r = 0,z>0) =7r gulR^oo =0 

Ao{r ^ cx),z ^ cx)) =Vo + q/r Ao{r = 0,z) =0 

(31) 

where i? = (2^ + r^)^/^ and where we have used the notation |^ = etc. Note 
that the field g is undefined at the origin and the resulting discontinuity of g at that 
point is an artefact of the coordinate system used. The asymptotic behaviour of Aq 
at infinity will be discussed in a later Section. To solve the equations numerically, it 
is more conveniant to use the field rather than b; this is why we have expressed 
the boundary condition in terms of that field. On the other hand, the equations take 
a simpler form when written in term of b, so we shall still use it below. 



/(0,0) =7r 

g{r = 0,z<0) =0 

a{r = 0,2;) =1 

ttrir = 0,z) = 
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Figure 3: Magnetic Field of the gauged Skyrmion (Aq = Ai = A2 = 1). 

£'(Ao, 1, 1) in Figure l.a is located in the region where the energy is virtually equal to 
the asymptotic value E{Xo, 1, 1) = 1.232. This means that a 37r2l.232F^/M„ ^ 7.2 
and that Aq ~ 1138. We can thus conclude that the effective impact of the Maxwell 
term we have added to the Skyrme model is relatively small. 

This justifies the procedure used in [5] where the Skyrmion was coupled with 
an external magnetic field of a magnetic monopole. Indeed, as the Maxwell field 
generated by a Skyrmion is very small (for the parameters fitting the actual mass of 
the nucleons) the external field is much larger than the Skyrmion's magnetic field. 

It would be interesting to find the differences between the electromagnetic quan- 
tities obtained from the ungauged model, as in [3], and our U{1) gauged model. We 
are not able to compute the solutions of the U{1) gauged model for the physical value 
of the parameter Aq as this is too large, but as we now know that the influence of 
the gauge field is very small, we can compute the latter perturbatively around the 
(ungauged) Hedgehog as an induced field. This enables the evalution of the energy 
correction and the induced magnetic moment. This perturbative analysis will be be 
carried out in the Section 4. 

It can also be concluded that if the U{1) gauged Skyrme model were quantised as 
in [3] (by quantising the zero modes corresponding to the global gauge transforma- 
tion) but taking into account the electromagnetic field generated classically by the 
Skyrmion, the result would not differ very markedly from what was obtained in [3]. 
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Figure 1: a) Energy and topological bound of the gauged Skyrmion in units of 247r^. 
b) Ratio of the electromagnetic and the total energy as a function of Aq . 




Figure 2: a) Energy density for the gauged Skyrmion in the (r, z) plane (Aq = Ai = 
A2 = 1). b) Energy density level curve for the gauged Skyrmion (Aq = Ai = A2 = 1). 
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Now the volume integral (with the appropriate normalisation of 127r^) of qq given 
by (6) is the degree of the map, or, the Baryon number. It is straightforward to verify 
that when the Ansatz (28) is substituted in qq and the volume integral is computed 
subject to the boundary conditions given above, the result will equal the integer n 
defined in (28). Thus, the Baryon number of the field configuration (28) equals the 
vortex number n. In what follows, we will restrict ourselves to unit Baryon number, 
n = 1, i.e. to the Nucleons. 

Before we proceed to substitute the Anstaz (28), (29) into the field equations, we 
calculate the Baryon current (11) for the field configurations (28), (29) described by 
the solutions we seek. We express the space-like part of this current ji in the radial 
direction flowing out of the normal to the surface of the cylinder which we denote by 
jr, and in the z direction which we denote by j^. The result is 

6 

jr = -[{fg^-gf^)asm'^fsmg (32) 
r 

+ 2 (i'"^ ~ "-dz) sin / cos / sin ^ + - {af^ - fa^) cos g] 
6 

jz = — [{fgr-gfr)asm'^fsmg (33) 
r 

+^{ga'r - agr) sin f cos f sin g + -{afr - fa^) cosg] , 

where we have denoted H = / etc. Note that the Baryon current (32), (33) are not 
sensitive to the charge of the Nucleon, i.e. that the function b{r,z) does not feature 
in them. 

We now turn to the equations to be solved, namely the Euler-Lagrange equations 
arising from the variational principle applied to the Lagrangian (4), in the static limit. 
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Substituting the ansatz (28) (29) into these equations of motion leads to 

arr -f - a.. - asin'if) sin2(^) + ^[/^ + + (^2 ^ ^2^ ,i^2^^^]] ^ q 

A/ - (^,^ + g! + ^ sin^ ^) sin / cos /+ 

+4^ sin / [{gr{fzz9r - fr9zz + fz9rz " /r^^^ " fz/rgz) 

+9z{frr9z - fz9rr + fr9rz " /rzS'r + fr/rgz))sm f + (/^^^ - grfz)"^ COS / 

+ sin ^/r2((a2 - b'^){f^ + /| - 2(^2 ^ sin^ /) cos / sin ^ 

+ ((a2 - 62)(/^^ + /^^ - fr/r) + 2/^(aa^ - bbr) + 2/^(aa^ - 66^)) sin/ sin ^ 
+2(a2 - 62)(/^^^ + /^^^)sin/cos^)] = 0. 

^9 +2ifr9r + fz9z) cot f - sing cos g+ 
+442 

[/^(/zfl'rr frr9z ~l~ frz9r fr9rz ~l~ fz9r/l") 

~^ fr{fr9zz fzz9r ~l~ frz9z fz9rz fz9z/'^) 

+ Sin ^/r2((a2 - 62)((^2 + ^2) ,i^2 ^ _ ^2 _ ^2) .^.^ 

((a^ - 62)(^^^ + g^^ - g^/r) + 25-^(00^ - bbr) + '2gz{aaz - bbj) sin V sin ^ 
+4(a2 - 62)(/^^^ + /^^^)cos/sin/sin^)] = 0. 

(34) 

In the case of the = gauge, equations (34) coincide with the Euler-Lagrange 
equations derived from the positive definite Hamiltonian density (5). Moreover in 
that case, those equations also coincide with the Euler-Lagrange equations of the 
reduced two dimensional Hamiltonian obtained by subjecting (5) to axial symmetry 
by substituting the Ansatz (28)-(29) into it. This is expected due to the strict impo- 
sition of symmetry. In the 7^ gauge, the Euler-Lagrange equations are derived 
from the Lagrangian (4) which is not positive definite. Nonetheless these equations 
coincide with those arising from the reduced two dimensional Lagrangian obtained 
by subjecting the Lagrangian (4) to axial symmetry. (This happens also for the Julia 
Zee dyon [9].) 

3.1 ^0^0: [/(I) Skyrme soliton 

It is easy to see from (34) that there are solutions for which b = {i.e. Aq = 0). 
As mentioned before, in that case, equation (34) can be obtained by minimising the 
Hamiltonian (30). Notice also that setting a = is not compatible with our boundary 
conditions {A^ would not be well defined at the origin). We thus expect our gauged 
solution to carry a non-zero magnetic field. 

To show this we have to solve equations (34) numerically for the non-vanishing 
functions f{r,z),g{r,z) and a{r,z). 

We have restricted our numerical integrations to the case where the vortex number 
n appearing in the axially symmetric Ansatz (28) is equal to 1, z.e.our soliton carries 
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unit Baryon number. 

Using (30), we have found numerically that £'(1,1,1) = 247r^l.01 whereas for 
those values of Ao,Ai,A2 the lower bound for the energy given by (26) is 247r^0.555. 
In Figure 1, we present the total energy for the gauged Skyrmion as a function of 
Aq, together with the lower bound given by (26). Note that the asymptotic value 
of £'(Ao,l,l) is 247r^l.232 as Aq — > oo. As a comparison, the energy (27) for the 
ungauged Skyrmion is Esk{^, 1) = 247r^ 1.232 with a lower bound set at 247r^. We see 
that E{Xo = oo, 1, 1) = Esk{i, 1) which means that as Aq — > oo, the gauge coupling 

1 /2 

1/Aq' goes to zero and the gauged Skyrmion becomes in this limit the ungauged 
Skyrmion. 

It is interesting to note that the energy of gauged Skyrmion is smaller than the 
energy of the ungauged Skyrmion, as expected, but that on the other hand, the 
amount by which the energy of the gauged Skyrmion exceeds its topological lower 
bound is larger than the excess of the energy of the ungauged Skyrmion above its 
respective topological lower bound. For example we can clearly see from Figure l.a 
that at Aq = 20, the energy of the gauged Skyrmion 1.22 (in units of 247r^) exceeds 
the lower boimd 0.95 by 0.27. This is larger than 0.232, the excess of the ungauged 
Skyrmion energy over its lower bound. For smaller values of Aq Figure l.a shows that 
the excess of the energy of the gauged Skyrmion over its lower bound is even larger, 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, i e.the term proportional to Aq 
in (30), is decreasing as Aq increases. Notice that we could have used for the Maxwell 
Energy the sum of all the terms involving the gauge field functions a and b in (30), 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the energy density of the 
Skyrmion in the r,z plane for Aq = 1. One sees clearly that the effect of the gauged 
field is to make the Skyrmion elongated along the z axis. The magnetic field vectors 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we show the configuration 
of magnetic field using arrows to represent the magnetic field vector at each point on 
the grid. Notice that there is a vortex around the point r = 2, z = 0. The magnetic 
field is thus generated by a current flowing on a ring centred around the z axis. 

In terms of the usual physical constants[3], we have Ag^ = 4e^, Ai = and 
= Sa^ where we use a instead of the traditional e for the Skyrme coefficient to 
avoid confusion with the electric charge. 

In our units, c = h = 1, we have e = (47rQ;)^/^ where a = 1/137 is the fine 
structure constant. Choosing F^^ = 186MeV, we can find the value for a by requiring 
that the energy of the neutron M„ = 939MeV matches the energy of the Skyrmion: 

M„ = g^(?;^,i.i). 

In Figure l.a, we can read the value of E{Xo, 1, 1) (given in units of 2A'w'^MeV) with 
Aq = 2a^/e^. We now have to find the value of Aq for which 247r^£'(Ao, 1, 1) = 
2(2Ao)^/^eM„/F^. The intersection between the curve g|^(2Ao)^/^ and the curve 
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Abstract 

We use a prescription to gauge the SU{2) Skyrme model with a U{1) field, 
characterised by a conserved Baryonic current. This model reverts to the usual 
Skyrme model in the limit of the gauge coupling constant vanishing. Wc show 
that there exist axially symmetric static solutions with zero magnetic charge, 
which can be electrically cither charged or uncharged. The energies of the 
(uncharged) gauged Skyrmions are less than the energy of the (usual) ungauged 
Skyrmion. For physical values of the parameters the impact of the U{1) field 
is very small, so that it can be treated as a perturbation to the (ungauged) 
spherically symmetric Hedgehog. This allows the perturbative calculation of 
the magnetic moment. 
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1 Introduction 



For a long time now, much attention has been paid to the Skyrme [1] model in 3 
dimensions. It is believed to be an effective theory for nucleons in the large N limit 
of QCD at low energies. The classical properties as well as the quantum properties 
of the model are in relatively good agreement with the observed properties of small 
nuclei [2, 3, 4]. 

Gauged Skyrme models have been used in the past. The U{1) gauged model [2, 5] 
was used to study the decay of nucleons in the vicinity of a monopole [5] , while the 
SU(2)l gauged model [7] was used to study the decay of nucleons when the Skyrme 
model is coupled to the weak interactions [7]. The Skyrme model has also been used 
to compute the quantum properties of the Skyrmion [3] where the gauge degrees of 
freedom are quantised to compute the low energy eigenstates of a Skyrmion. These 
states were identified as the proton, the neutron and the delta. 

The aim of this work is to show that the Skyrme model can be coupled to a self 
contained electromagnetic field and that this U{1) gauged model has stable classical 
solutions. In addition to these solitons with vanishing magnetic and electric flux, we 
show that this system supports solutions with nonvanishing electric flux which are 
analogous to the dyon solutions of the Georgi-Glashow model, just as the uncharged 
solitons are the analogues of the monoploles [8] of that model. The electrically charged 
lumps have larger energy, or mass, than the uncharged soliton, just like the .luIia-Zcc 
dyon [9] has larger energy, or is heavier, than the (electrically uncharged) monopole. 
Wc shall refer to these lumps as charged U{1) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Maxwell gauged Skyrme 
model, the present work is also an example of a soliton in a li-dimensional SO{N) 
gauged S'' model with N < d ion the case d = 3, N = 2 , extending the results of 
Ref. [10] which were restricted to the N = d cases. (The work of Ref. [10] consists 
of establishing topological lower bounds for the generic case, encompassing earlier 
examples in two [11] and three [12, 13] dimensions respectively.) The gauging pre- 
scription used hero by us coincides precisely with that used in Ref. [5] and permits 
the establishing of a topological lower bound which did not feature in Ref. [5] and 
which is carried out here to establish the stability of the soliton. Such lower bounds 
are absent in the other prescription of gauging the Skyrme model as in Refs. [7]. 
(Notice thai wc name the sigma models after the manifold in which the fields take 
their values rather than using the name of the symmetry group for the model. Thus 
what is sometime called the 0{d+ 1) model in the literature will be refered to as the 
S" model.) 

The U(l) gauged SU{2) Skyrme model is described by the Lagrangian [5] 
C = ^Tr{D,UD,U^) - ^,Tr{[(D,U)UK (D.up^f - 

where the U{1) gauge covariant derivative is 

Di,U = di,U + ieAiAQ,U], (1) 
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Ao to 0, and impose the condition that the electric field is assymptoticaly like that 
of the proton. This pertubation method would not make much sense though as one 
would expect the electric field to be quite large close to the skyrmion. 

The relevant energy functional is (23), and the resulting equation arising from the 
variation of the gauge field At is 

\d,F,,=j, (36) 
ji = -e"^ (iAi(4^A<?" + As \{dj\<p-^)D,ifDjd^ + 2(j^D[i<?-"D;]0-*D;(i6'*]) . (37) 

We are concerned here with the case where Ai = Q (37) and the chiral field (jf = 
{(j)" , (j)^ , 4>^) in (37) describes the Hedgehog, i.e. 

(j," = sin F{R) x" , 4? = sin F{R) i^ c^" = cos F{R) , (38) 



where now x"' = x''/R, with R — ^/r^~+^. By virtue of equation (36) the current 
(37), given by (38) and Ai = 0, will now induce a (small) U(l) field Aj, with curvature 

The shift in the energy of the Hedgehog due to the induced U{1) field Aj is 

AE = y ifx {XoFijFij + 4Ai ji) , (39) 

in which = ji(0) is the current (37) for Ai — 0. (Note that all quantities evaluated at 
Ai = are denoted by Roman script, e.g. j; = ji(0), as well as the induced connection 
and curvature A; and Fy.) When equation (36) is satified for the induced (7(1) field, 

AE = -Ao j d\ FjjFjj (40) 

= 2 j (fxAiii, (41) 

which is, as expected, a strictly negative quantity. 

The current j; = (ja, ja) in (39) and (41) is given, for the Hedghog field configura- 
tion (38), by 

sin"F/Ai , / ,2 sin'FAA , 

j3 = 0. (43) 

We now note that 9; j; = 0, which means that equation (36), for Ai = 0, takes the 
following form in terms of the induced U(l) connection Aj = (Aq, 0) 

AoA A„ = -j„ . (44) 

The solution is well known and can be expressed, using the obvious notation ja(x) = 
j(i?) Sa/iX/i in terms of (42), as 

1 r 1 

Aa(x) = --^£ai} / 1 71 i{R') f^x' , (45) 

47rAo J |x — X I 
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The flux of Oj vanishes, as can deduced by anticipating the finite energy conditions 
to be stated later. 

Note that the 3-volume integral of Qq in (8) is the degree of the map for tht 
ungauged system namely the Baryon number. 

Identifying q, (9), with the naught component j" of the Baryon current, is 
defined by 

f = el'"''' E,,,A4>''d,<t>'d,<t>'<t>' + 7, dM,s^''<l>''d.<^^) (H! 

The 4-divergence of (11) receives a contribution only from its first term, whicf 
being locally a total divergence implies that the 3-volume integral of j° is a conservec 
quantity. Alternatively we consider the 4-divergence of (12), 

d^f = 66''-''"'e,f,6ABDi,<j>''D,<j>CDf,<j>^D,<j>'' [13] 

which is analogous to the corresponding quantity in the work of Goldstone anc 
Wilczek [6]. This contrasts with the expression for the total divergence of the topolog- 
ical current in the work of d'Hoker and Farhi [7], where a different gauging prescriptior 
is used leading to that quantity being equal to the local anomaly. 

We now ijrocoed to find a model whose Hamiltonian "Ho is bounded from below bj 
the topological charge density defined by (9). We will then show that the Hamiltoniar 
(5) is given by Ho plus certain positive definite terms. 

First of all, we reproduce the density Qg, (7), in (9) by using the following in- 
equality 

(KsDir - eijte''''"'4Dj(l>''D;,4,'4,''f > (14; 

where the two constants K3 and K2 have the dimensions of length. Expanding thf 
square, we get Qg on the right hand side of 

i4(Di4,''f + 4{D,^4,"Dj]4/'f > 2K.n,KlQG- (15; 

To reproduce the other term in (9), \£ijkFij{£'^^4'^dk4>^), we use the following 
inequality 

(K^Fi^ - iK4eijfce'**0*Dfc(>'*)' > (16; 

yielding 

<Ff, + Kl^-{e^''^''D,<t>y > 4K,e,,,F„{e^W'D,^^). (If 

With the special choice for the relative values of the constants 3K3K2 = tht 
sum of (13) and (15) yields the following 

4f^, + 4(Dirr + «^(i?[,0"-D,]0')^ + ^(e^«0«A0^)^ > 2«3«^e. (is; 
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The right hand side of (18) is now proportional to the topological charge density 
Q defined by (9) so that the inequality (18) can be interpreted as the topological 
inequality giving the lower bound on the energy density functional if we define the 
latter to be the left hand side of (18), namely 

no = 4Ffi + ^m-f + 4iDi.rDj]^r + ^(e-^^^^A^-^)^. (19) 

The Hamiltonian system (19) is almost the Hamiltonian of the gauged Skyrme 
model (5) (remember that Aq = 0). It differs from the latter only in its last term. 
Now we can use the identity 

(6^''<j>''D,<j>^f = (A<^°)' - [^(<?i'"A/')' + (0'"A<^'*')'] (20) 

and add the positive definite term ^^j^ |j (<?>[" A^''')^ + (^'"A^'^')^! appearing on the 
right hand side of (20) to Ho in (19) to end up with the Hamiltonian for the (7(1) 
gauged Skyrme model; 

H = KlFf^ + /c?(A0")' + 4{D[irDj^<l,''f > 2k,4q (21) 

which is nothing but the static Hamiltonian (5) in the temporal gauge = 0, and 
where 

By virtue of (18), (21) is also bounded from below by 2k3K.I(), namely by a number 
proportional to the topological charge density q. 

We thus see that Tio can be considered as a minimal {U{1) gauged) model, but 
from now on, we will restrict our attention to the physically more relevant model (21) 
and integrate it numerically to find its topologically stable finite energy solitons. 

The soliton solutions to the system (21) can only be found by solving the second- 
order Euler-Lagrange equations, and not some first-order Bogomol'nyi equations since 
saturating the inequalities (14) and (16) would not saturate the lower bound on the 
energy density functional "H. In this context we note that saturating (14) and (16) 
does indeed saturate the topological lower bound on the functional "Ho by virtue of 
the inequality (18), and should it have turned out that the Bogomol'nyi equations 
arising from the saturation of (14) and (16) supported non-trivial solutions, then Ho 
would have been a very interesting system to consider. As it turns out however, these 
Bogomol'nyi equations have only trivial solutions in exactly the same way as in the 
case of the (ungauged) Skyrme model [1]. 

The energy for the static configuration, when the electric field vanishes, is ex- 
pressed as 

EiXo, Ai, As) = / d'x [Ao4 + \i{D,rf + X^iDy,^" D .^tj^'f] (23) 
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have used. To compute the static solutions, we have employed a relaxation method 
using finite differences on a regular grid {dr = dz). This discretisation method is 
similar to the one employed in the numerical computation of the solutions of Skyrme 
models in 2-1-1 and 3-1-1 dimensions [16, 17, 18]. To compute the electrically charged 
solutions, we have imposed the boundary condition 6(oo, oo) — Vq for different V'o 
and using a dichotomic method, we have determined the values that give a solution 
with the same electric flux as the proton. Most of the simulations where done on 
200 X 400 or 300 x 600 grids. By computing the same solution for various lattice sizes, 
we have empirically obtained the following relation for the expression of the energy 
of a solution : E = Eo + Kdr'^ where Eo is the exact solution, dr = dz the lattice 
spacing and K \s a constant which depends on Aq but takes values between 1 and 
0.5. We see thus that the energies we have obtained are accurate to within one or one 
half of a percent. This inaccuracy in the value of the energy is comparable to that 
of many other similar works on 2 dimensional systems [17, 18], and though it might 
look large, it does not affect any of the conclusions we have drawn. 

4 Perturbation around the Hedgehog 

We have seen from the work of 3.1 above that the energy of the U{1) gauged Skyrmion 
for the physical values of the parameters, namely of the pion decay constant and the 
(7(1) coupling, does not differ significantly from that of the ungauged Hedgehog. It 
is therefore justified, for these values of the parameters, to treat the (7(1) field as a 
perturbation to the Hedgehog in the same way as Klinhamer and Manton [14] treat 
the (7(1) field as a pertubation to the SU{2) sphaleron. We will then be able to 
compute the magnetic moment of the Neutron, as well as the (small) deviation of its 
mass from that of the Hedgehog. 

The equations for the fields ((^°,^j,) are derived from the Lagrangian (1). The 
equation for A/^ will be of the form 

Ao d,F^, = . (35) 

The method consists of setting the gauge field A^j to zero in the current jf, in (35) (and 
in the equation for 0°) and calculating the resulting induced electromagnetic field A^^. 
The gauge field computed this way can then be interpeted as the U (1) field generated 
by the (ungauged) Hedgehog Skyrme field. With this perturbative procedure, it is 
possible to calculate the induced static magnetic potential A; {i = 1, 2, 3), but not the 
static electric potential Aq, which in this scheme vanishes and can only be calculated 
non-perturbatively. The reason simply is that restricting to the use of the static 
Hedgehog, the zeroth component of the current jo at A^ = vanishes, resulting in 
turn in vanishing induced potential Aq according to (35). 

As a consequence the electric field will be idcnticaly zero, which implies that wo 
can derive the equation from the static Hamiltonian rather than from the Lagrangian. 
Notice also that we could try to compute perturbatively a solution for the electrically 
charges skyrmion by keeping in jo the terms proportional to Ao, instead of setting 
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Figure 1: a) Energy and topological bound of the gauged Skyrmion in units of 247r^. 
b) Ratio of the electromagnetic and the total energy as a function of Aq . 




Figure 2: a) Energy density for the gauged Skyrmion in the (r, z) plane (Aq = Ai = 
A2 = 1). b) Energy density level curve for the gauged Skyrmion (Aq = Ai = A2 = 1). 
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Now the volume integral (with the appropriate normalisation of 127r^) of go giver 
by (6) is the degree of the map, or, the Baryon number. It is straightforward to verifj 
that when the Ansatz (28) is substituted in Qq and the volume integral is computec 
subject to the boundary conditions given above, the result will equal the integer r 
defined in (28). Thus, the Baryon number of the field configuration (28) equals tli< 
vortex number n. In what follows, wc will restrict ourselves to unit Baryon number 
n = 1, ie.to the Nucleons. 

Before we proceed to substitute the Anstaz (28), (29) into the field equations, wf 
calculate the Baryon current (11) for the field configurations (28), (29) described b) 
the solutions we seek. We express the space-like part of this current j; in the radia 
direction flowing out of the normal to the surface of the cylinder which we denote bj 
jr, and in the z direction which we denote by j^. The result is 

if = -[(/S2 - 3/2)0 sin V sin p (32; 
r 

+ ^{9az - ag^) sin / cos / sin 5 + -(a/^ - fa^) cosg] 

= --[{fgr - gfr)asm^ fsing {33 
r 

+ 2^9(^r - agr) sin /cos/ sin (jr + -{afr - far)cosg] , 

where we have denoted = / etc. Note that the Baryon current (32), (33) are no 
sensitive to the charge of the Nucleon, i.e. that the function b{r, z) does not featurt 
in them. 

We now turn to the equations to be solved, namely the Euler-Lagrange equations 
arising from the variational principle applied to the Lagrangian (4), in the static limit 
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Substituting the ansatz (28) (29) into these equations of motion leads to 
a.. -'f - a.. - asm'if) sin^(fi) [i^^ + ^[/,^ + /| + (g^^ + gl) sin^(/)]] = 
Kr -'f + ^ + b..- bsin'if) sin^(fi) [Ai^ + + /I + (p^2 + gl) si„^(/)]] = 

A/ -(9,^ + 9,^ + ^sin^9)sin/cos/+ 

+4|^ sin/[((7r(/zzSr - /ri7z2 + /zSrz - frzOz - fz/rOz) 

+9z0rr9z - fzOrr + frOrz - frzOr + fr/rgz))sill f + (/rfi^ - flr/z)^ COS / 

+ Sin gy {(a' - 6^) {f^ + fl - 2(gl + gl) sin^ /) cos / sin g 

+ ((a^ - 6^)(/„ + fzz - fr/r) + 2fr{aar - bbr) + 2f^(aa^ - 66^)) sin / sing 
+2(a2 - 6')(/,fi. + /,fc)sin/cosfl)] = 0. 

Afi +2(/r5i, + /25^)cot/ - ^!^sin5Cosg+ 

|^/£(/25rr frrQz ~l~ frzQr frQrz ~l~ fz9r/^^ 

+ fr{frOzz - fzzOr + /rzi7z - fzOrz - fzOz/r) 

+ sin g/r^(^(a'^ - b'^){(gl + gl) siv? f - - /|) cos^i 

((a^ - b'^)(grT + gzz - gr/r) + 2gr(aar - bb^) + 2g^(aa^ - bb ^)) sin^ f sing 
+4(a^ - 62)(/rgr + /2g2) cos/sin/sin s)] = 0. 

(34) 

In the case of the Aq = gauge, equations (34) coincide with the Euler-Lagrange 
equations derived from the positive definite Hamiltonian density (5). Moreover in 
that case, those equations also coincide with the Euler-Lagrange equations of the 
reduced two dimensional Hamiltonian obtained by subjecting (5) to axial symmetry 
by substituting the Ansatz (28)-(29) into it. This is expected due to the strict impo- 
sition of symmetry. In the j4o 7^ gauge, the Euler-Lagrange equations are derived 
from the Lagrangian (4) which is not positive definite. Nonetheless these equations 
coincide with those arising from the reduced two dimensional Lagrangian obtained 
by subjecting the Lagrangian (4) to axial symmetry. (This happens also for the Julia 
Zee dyon [9].) 

3.1 ^0 = 0: [/(I) Skyrme soliton 

It is easy to see from (34) that there are solutions for which 6 = {i.e. Ao = 0). 
As mentioned before, in that case, equation (34) can be obtained by minimising the 
Hamiltonian (30). Notice also that setting a = is not compatible with our boundary 
conditions (Ai would not be well defined at the origin). We thus expect our gauged 
solution to carry a non-zero magnetic field. 

To show this wo have to solve equations (34) numerically for the non-vanishing 
functions f{r, z), g{r, z) and a{r, z). 

We have restricted our numerical integrations to the case where the vortex number 
n appearing in the axially symmetric Ansatz (28) is equal to 1, ie.our soliton carries 



10 



unit Baryon number. 

Using (30), we have found numerically that £(1,1,1) = 247r^l.01 whereas for 
those values of Aq, Ai, A2 the lower bound for the energy given by (26) is 247r^0.555. 
In Figure 1, we present the total energy for the gauged Skyrmion as a function of 
Ao, together with the lower bound given by (26). Note that the asymptotic value 
of i?(Ao,l,l) is 247r^l.232 as Aq — s- 00. As a comparison, the energy (27) for the 
ungauged Skyrmion is £,*(!, 1) = 247r^l.232 with a lower bound set at 247r^. We see 
that E(Xo = 00, 1, 1) = Eski^, 1) which means that as Ao 00, the gauge coupling 
I/aJ^^ goes to zero and the gauged Skyrmion becomes in this limit the ungauged 
Skyrmion. 

It is interesting to note that the energy of gauged Skyrmion is smaller than the 
energy of the ungauged Skyrmion, as expected, but that on the other hand, the 
amount by which the energy of the gauged Skyrmion exceeds its topological lower 
bound is larger llian the excess of the energy of the ungauged Skyrmion above its 
respective topological lower bound. For example we can clearly see from Figure l.a 
that at Ao = 20, the energy of the gauged Skyrmion 1.22 (in units of 247r^) exceeds 
the lower bound 0.95 by 0.27. This is larger than 0.232, the excess of the ungauged 
Skyrmion energy over its lower bound. For smaller values of Ao Figure l.a shows that 
the excess of the energy of the gauged Skyrmion over its lower bound is even larger, 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, ie.the term proportional to Aq 
in (30), is decreasing as Aq increases. Notice that we could have used for the Maxwell 
Energy the sum of all the terms involving the gauge field functions a and b in (30), 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the energy density of the 
Skyrmion in the r, z plane for Aq = 1. One sees clearly that the effect of the gauged 
field is to make the Skyrmion elongated along the z axis. The magnetic field vectors 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we show the configuration 
of magnetic field using arrows to represent the magnetic field vector at each point on 
the grid. Notice that there is a vortex around the point r = 2, 2; = 0. The magnetic 
field is thus generated by a current flowing on a ring centred around the z axis. 

In terms of the usual physical constants[3], we have Aq ' = 4e^, \\ = F^/?: and 
A2 ^ = 8a^ where we use a instead of the traditional e for the Skyrme coefficient to 
avoid confusion with the electric charge. 

In our units, c = fi = 1, wc have e = (47ra)'^^ where a — 1/137 is the fine 
structure constant. Choosing = 186Mey, we can find the value for a by requiring 
that the energy of the neutron M„ = 939Mey matches the energy of the Skyrmion: 

F 2(1^ 
M„ = ^£;(^,1,1). 

In Figure La, we can read the value of E{\, 1, 1) (given in units of 24jr^i'VfeV') with 
Ao = 2(!?jf?. We now have to find the value of Ao for which 247r^ii(Ao, 1, 1) = 
2(2Ao)^^^eM„/F^. The intersection between the curve (2Ao)^^^ and the curve 
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with n° = (sin cosn0) in terms of the azimuthal angle (j) and = Xa/r. n in n° 
is the vorticity, which for the Nucleons of interest to us here, equals one, n = 1. The 
functions a, 6, Ci, C2, / and g both depend on r and z. 

Our Ansatz (29) for the U{1) field consists of decomposing the latter in the most 
general tensor basis possible. We will find out below, when we compute the Euler- 
Lagrango equations, that the functions Ci and C2 vanish identically. Anticipating this, 
we suppress them henceforth. In its final form this Ansatz agrees with that used in 
Ref. [15], the latter being arrived at by specialising the Rebbi-Rossi Ansatz for the 
axially symmetric S0(3) field. 

The static Hamiltonian, ie.the Too component of the energy mometum tensor Tj,„, 
is then given by 

H= jl^^^al + al + {K-^r + bl\ 

+2A2 sin^ / [(As. - hgrf + sin^ 9[^(/,^ + fl + {gl + gl) sin^ /)]] ^drdz. 

(30) 

The boundary conditions for the Skyrmion fields are the same as the boundary 
conditions for the hedgehog ansatz when expressed in the cylindrical coordinates 
where g — 7r/2 + arctan(.z/r) and defining R — {z^ + r^Y^^, f = f{R) w'ith /(O) — tt 
and limg^oo f{R) ~ 0. Prom this we can deduce that the function / has a fixed value 
at the origin and at infinity. For smothness along the z axis, each field, that is /, g and 
Aa must satisfy the condition that the partial derivative with respect to of the field at 
r = vanishes. The boundary conditions and the asymptotic behaviours for a and b 
are chosen so that the gauge fields are well defined and Aq looks asymptoticaly like 
a coulomb field (i.e. with an electric charge but no magnetic charge). We also require 
that the total energy be finite. These conditions leads to the following constraints: 



/(0,0) =7r f(r^oo,z^oo) =0 /r(r = 0,z) =0 

g{r = 0,z<0) =0 g{r = 0, z > 0) = tt teU^oc =0 

a{r = 0,z) =1 Orloo =0 cizloo =0 

ar{r = 0,z) =0 Ao(r ^ (x,z (x) =Vo + q/r Ao(r = 0,z) =0 

(31) 

where R = (z^ + r^Y^^ and where we have used the notation |^ = etc. Note 
that the field g is undefined at the origin and the resulting discontinuity of g at that 
point is an artefact of the coordinate system used. The asymptotic behaviour of Aq 
at infinity will be discussed in a later Section. To solve the equations numerically, it 
is more conveniant to use the field ^0 rather than b; this is why we have expressed 
the boundary condition in terms of that field. On the other hand, the equations take 
a simpler form when written in term of 6, so we shall still use it below. 
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Figure 3: Magnetic Field of the gauged Skyrmion (Ao = Ai = A2 = 1). 



E{\o, 1, 1) in Figure l.a is located in the region where the energy is virtually equal tc 
the asymptotic value E{\o, 1, 1) = 1.232. This means that Ri 37r^l.232F,j/Af„ ^ 7.^ 
and that Ao ~ 1138. We can thus conclude that the cfl'cctivc impact of the Maxwcl 
term we have added to the Skyrmo model is relatively small. 

This justifies the procedure used in [5] where the Skyrmion was coupled witl 
an external magnetic field of a magnetic monopole. Indeed, as the Maxwell fielc 
generated by a Skyrmion is very small (for the parameters fitting the actual mass o: 
the nucleons) the external field is much larger than the Skyrmion's magnetic field. 

It would be interesting to find the differences between the electromagnetic quan- 
tities obtained from the ungauged model, as in [3]. and our U{1) gauged model. We 
are not able to compute the solutions of the U{1) gauged model for the physical valuf 
of the parameter Ao as this is too large, but as we now know that the infiuence o: 
the gauge field is very small, we can compute the latter perturbatively around thf 
(ungauged) Hedgehog as an induced field. This enables the evalution of the energj 
correction and the induced magnetic moment. This perturbative analysis will be bt 
carried out in the Section 4. 

It can also be concluded that if the U{1) gauged Skyrme model were quantised as 
in [3] (by quantising the zero modes corresponding to the global gauge transforma- 
tion) but taking into account the electromagnetic field generated classically by th( 
Skyrmion, the result would not differ very markedly from what was obtained in [3]. 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a function of Aq 
3.2 Ao^O : charged [7(1) Skyrmion 

We can now look for solutions with a non-zero electric charge by requiring that the 
field b in our ansatz (29) does not vanish. To do this we follow the same procedure 
as Julia and Zee [9] and require that the electric field be asymptotically of the form 

^0 = Vo + q/{r'^ + z'^)'-^^ where V'o and q are two constants. In practice, one computes 
solutions for different values of Vo and evaluates q by computing the electric flux. We 
sought only those solutions, for which the electric flux equals Air times the charge of 
the electron. 

It is important to realise that in this case, equations (34) are obtained after min- 
imising the action and thus they do not minimise the Hamiltonian (30). 

In our units, the charge of the electron is 0.303. In Figure 4 we show the energy 
as a function of Aq, as well as Vq as a function of Aq, so that q = 0.303. 

One can see that, for a fixed value of Aq, the energy of the charged gauged Skyrmion 
is smaller than the energy of the ungauged Skyrmion when Aq < 7 but it is always 
larger than the energy of the uncharged gauged Skyrmion. If the energies of the 
electrically charged and uncharged gauged Skyrmions were interpreted as the the 
masses of the Proton and the Neutron mp and mn, then on this purely classical level 
we would have to conclude that {mp — mfi) > which is not correct. This is expected 
on the basis of its analogy with the dyon [9]. Clearly, to calculate this mass difference 
correctly one would have to perform the collective coordinate quantisation as in Ref. 
[3], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is unfortunately very 
difhcult to carry out the numerical computations accurately when Aq is very large. 

At this stage, it is worth saying a few words about the numerical methods we 
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and performing the dilation x — > ax, ^ a ^yl^,, we get 

E(Ao, Ai, A2) = / d'x + aX,(D,rr + ^(C[i0"C^]0')']. (24) 

If we choose a = (f^)^'^ then we have 

E(X„,Xr,X2) = (XrX2y''E(^,l,l), (25) 

from which we see that we can set Ai = A2 = 1 without any loss of generality. By 
virtue of (21) and (25), we can finally state 



EiXo,XuX,)>2i-^Y/' [d'xs. (26) 



A1A2 

Notice that for the usual Skyrme model we have 

Esk{XuX2) = J dx'ix.idir? + x.,(diirdj]^''y] 

= (X,X2Y''E.,(1,1) (27) 

> 2(AiA2)i/Vdx3ft. 

We will use (27) to compare the numerical solutions of the gauged Skyrme model 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability considerations discussed 
in this Section apply only to the solutions with no electric field, i.e. with Aq = 0. 



3 The soliton and the charged U{1) Skyrmion 

To find the static solutions, wc have to look for the largest symmetry group of the 
functional to be subjected to the variational principle, and look for solutions which 
are invariant under that symmetry group. For the solutions in the Ao = gauge 
this is the static Hamiltonian (21), while for the solutions in the ^0 7^ it is the 
Lagrangian (4). For our choice of gauge group the largest symmetry is the SO(2) 
group corresponding to an axial rotation in space-time and a gauge transformation 
on the gauge field. Defining the axial variables r = \l ~c\ + x\ and and z = x-^ in terms 
of the coordinates Xi = [x^jXs), a = 1,2, the most general axially symmetric Ansatz 
[15] for the fields 0" = (0", 0^) (with a = 1, 2 and ^ = 3, 4), and, A = [A^, A3), is 

= sin / sin p n° , (j)^ = sin / cos g , <j>^ = cos / , (28) 
a(r,z)—n Calr^z) , Ci(r,z) , b(r,z) , . 
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Skyrme coupling). For physical values of these parameters in the model, we find that 
the energy (mass) of the gauged Skyrmion does not differ significantly from that of 
the ungauged charged-l Skyrmion, namely the familiar hedgehog [1]. This implies 
that for these values of the physical parameters, the U{1) gauged Skyrmion can be 
regarded as a perturbation of the (ungauged) hedgehog, enabling the computation 
of the magnetic moments of the gauged Skyrmion (i.e. the Neutron) and the shift 
of the energy of the gauged Skyrmion away from the energy of the hedgehog [1], 
perturbatively using the method employed by Klinkhamer and Manton [14] for the 
sphaleron of the Weinberg-Salam model. 

In Section 2, we define the topological charge and establish the corresponding 
lower bound on the energy functional. In Section 3 we present the solutions which 
have no electric fields in the first subsection and electrically charged solutions in 
the second subsection. The perturbation analysis of the gauged Skyrmion around 
the (ungauged) hedgehog is carried out in Section 4, and Section 5 is devoted to a 
discussion of our results. 



2 The topological charge and lower bound 

The definition of the topological charge is based on the criterion that it be equal to the 
Baryon number, namely the degree of the map. For the gauged theory however, this 
quantity must be gauge invariant as well. This requirement can be systematically [10] 
satisfied by arranging the gauge invariant topological charge density to be the sum 
of the usual, gauge variant winding number density 

^0 = e,,te'"'"'dj''d,^''dt^'<t>'' , (6) 

plus a total divergence whose surface integral vanishes due to the finite energy condi- 
tions, such that the combined density is gauge invariant. In 3 dimensions, this is given 
explicitly in Refs. [10, 13] for the 50(3) gauged model, and for the present case of 
interest, namely the SO{2) gauged model, the charge density can be derived from 
that of the SO(3) gauged model by contraction of the gauge group S0{3) down to 
S0{2). It can also be arrived at directly. To state the definition of the charge, we 
denote the gauge covariant counterpart of (6) by 

&G = e,,^""^D,d^'D,dfD,(l/^d^ , (7) 

so that using the notations (6) and (7) we have the definition of the gauge invariant 
topological charge 

Q = ft + SjOi, (8) 
= QG + \eijkFii{e'''',j,''D^,j,''). (9) 

In (8) the density 0; is the following gauge variant form 

Oi = aejj.e^^^ja,^'* 0^. (10) 
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(47; 



with [19] 

= E E ^|fr yii\S'A') Yi'\9,4>) . 
After performing the angular integrations wo have 

Aa(x) = -I{R) SapXp , (46 

with I(t){R) given by the integral 

Finally, in the 1 region of interest, the induced U{1) potential is 

/ - 1 

Aa(x) = eap , with , / = — / s^iis) ds , (48; 
XI oAq Jo 

to be evaluated numerically using the numerically constructed hedgehog profile func- 
tion F(x) (38). 

Comparing (48) with the usual Maxwell potential of a magnetic dipole 
we find that /i = (0, 0, fi) is 

n = Avi . (49; 

We can evaluate the magnetic moment (49) and the energy correction (41) inducec 
by the electromagnetic field by evaluating the integral (47) and (48) numerically. I: 
we take the experimental values F„ = 186MeV and a = 7.2 we obtain 

= 0.01393/m = 0.43nm (50; 
AE = -O.lA^eV. (51; 

The experimental value for the magnetic moment of the proton and the neutron art 
respectivly //p = 0.0902/m = 2.79nm and = 0.0617/m = — 1.91nr«. If on the 
other hand we take the values of the parameters derived in [3], F„ = 129MeV anc 
a = 5.45, we have 

fj. = 0.0468/m = 1.449nm (52; 
AE = -0.32keV. [53] 

The magnetic moment of a particle is strictly speaking a quantum property and it 
should be computed by quantising the SU(2) gauge degree of freedom as in [3]. Nev- 
ertheless, we see that if the take the parameters derived in [3] the classical magnetic 
moment is of the correct order of magnitude. The sign is of course undetermined as 
the classical magnetic moment is a vector. We can thus conclude that our model offers 
a reasonable classical description of Nucleons and affords a method for computing th( 
electromagnetic field generated by the Skyrmion, classically. It is quite surprising tc 
see that a quantum property like the magnetic moment can be reasonably predictec 
by a purely classical procedure. 
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5 Summary and discussion 



We have shown that the SU{2) Skyrme model gauged with U{1) has two types of 
finite energy static solutions, electrically uncharged and charged respectively. Both 
of these solutions are axially symmetric and carry no magnetic charge but support a 
magnetic field shaped like a torus centred around the axis of symmetry, albeit resulting 
in zero magnetic flux. The uncharged solutions, like the ungauged Skyrmion, have a 
topological lower bound. The electrically charged solutions are the analogues of the 
Julia-Zee dyons [9] of the Georgi-Glashow model. 

Concerning the stability of the electrically neutral solution, which is expected to be 
stable by virtue of the lower bound on the energy, we have not made any quantitative 
effort to test it. We expect however that the solitons of this gauged Skyrme model are 
stable, or that at least they have stable branches for all values of the parameters in the 
model. This expectation is based on our knowledge of the corresponding situation 
when the Skyrme model is gauged instead with 50(3) [20, 21], in which case the 
equations arising from the imposition of spherical symmetry were one dimensional, 
and hence technically much more amenable to the numerical integration. In that 
case it was found that in addition to stable branches of solutions, there were also 
some unstable branches bifurcating from the former, the important matter being 
that there were indeed stable branches of solutions, characterised by the (ranges of 
the) parameters of the model. It would be very interesting to carry out the analysis 
corresponding to that of [20, 21], for the considerably more complex case of the axially 
symmetric equations at hand. This however is technically beyond the scope of the 
present work. 

The energies of the gauged uncharged Skyrmions are smaller than the energy of the 
usual ungauged Skyrmion. When the gauge coupling I/aJ'^ goes to 0, the uncharged 
gauged Skyrmion tends to the ungauged Skyrmion. We also note that the energy of 
the electrically charged Skyrmion is higher than the uncharged one, just as the mass 
of the dyon is higher than that of the monopole of the Georgi-Glashow model. 

Perhaps the most interesting physical result of the present work is that when 
parameters in the model are fitted to reproduce physical quantities, it turns out that 
the effect of the Maxwell term in the Skyrme Lagrangian is very small. This is because 
for the physical value of the constant Aq = 1138, the energy of the gauged uncharged 
Skyrmion differs little from that of the ungauged Skyrmion, as seen from Figure lb. 
The gauged Skyrmion field itself is thus nearly radially symmetric (though the gauge 
field is not). 

Having found that the influence of the electromagnetic fleld on the Skyrmion 
is small, we were pointed in the direction of treating the magnetic potential as an 
induced field perturbatively around the (ungauged) Hedgehog. We have been able 
thus, to compute the classical magnetic moment of the (uncharged) Skyrmion of unit 
Baryon charge, namely of the Neutron. The result is that the classical magnetic 
moment of the Skyrmion matches surprisingly well to the experimental values of the 
magnetic moments of the Nucleons. 
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where = d/j^A^ — d^Ai^ and where the charge matrix of the quarks is expressed 
as Q = i(|l + 1T3). This differs from the covariant derivative of Ref [5] only in the 
unimportant matter of the sign of i in (1), which we have chosen for consistency of 
the convention used in Ref. [10]. 

In what follows it will be more convenient [10] to parametrise the Skyrme field as 
an valued field 0" = (0", 4>^), a = l,2,A = 3,i subject to the constraint [(^"l^ = 1. 
The two fields U and (j) are related to each other via the following expression 

U = 4,''t'', u-^ = (yt = 0V (2) 

where t" = (ia",ia^,l) and f" = (—ia°,—ia^,l), in terms of the Pauli matrices 

The gauge covariant derivative now can be re-expressed as 

D,r = diA" + ^^.e""/, Cf.^-* = a^.^-*. (3) 

where = eAi^ and F^^j, = ej-^^. 

The Lagrangian for the U{1) gauged Skyrme model can then be written as 

C = -AoF^l + \i\D^r? - MDy^rD^-^ct-'? (4) 

where the square brackets on the indices imply (total) antisymmetrisation and where 
Aq ' = 4e^, Ai = F^l<i and Aj' = 8g^. The late Greek indices fi label the Minkowskian 
coordinates, while the early Greek indices a = 1, 2 and the upper case Latin indices 
^ = 3, 4 label the fields 0" = (0", (/.■*). 

The static Hamiltonian pertaining to the Lagrangian (4) is 

■H = \oF^j + AilA^I' + X2\DiirDj^,l>''\^ 

(5) 

+i2X,\d,A,p + Al{X,\r? + 16X2[\r?\dJ^? + i|a.(l0-*P)l']}, 

where the indices i = a,3 label the space-like coordinates. 

To find the static solutions of the model, one would usualy solve the Euler La- 
grange equations which minimise the Hamiltonian (5), but because of the electric 
potential Ao, one must solve the Euler Lagrange equations derived from the La- 
grangian (4). We then look for static solutions, but, as for the Julia-Zee dyon [9], we 
have to impose the proper asymptotic behaviour for the electric potential to obtain 
static solutions which are electrically charged (in the classical sense, i.e. solutions 
where the flux of the electric fleld is non zero). 

When the full equations of motion are written down, one flnds as expected that 
there are static solutions for which Ao — 0, i.e. solutions for which the electric fleld is 
identically zero. For these solutions in the temporal gauge, the equations of motion 
reduce to the equations obtained by minimising the Hamiltonian (5). We study the 
solutions of unit Baryon charge of the U{V] gauged Skyrnio model with and whithout 
an electric field, for various values of the C'^(l) coupling constant (or equivalently the 
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contained electromagnetic field and that this U{1) gauged model 
solutions. In addition to these solitons with vanishing magnetic c 
show that this system supports solutions with nonvanishing elei 
analogous to the dyon solutions of the Georgi-Glashow model, ju 
solitons are the analogues of the monoploles [8] of that model. The 
lumps have larger energy, or mass, than the uncharged soliton, ju 
dyon [9] has larger energy, or is heavier, than the (electrically un 
We shall refer to these lumps as charged U(l) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Max\ 
model, the present work is also an example of a soliton in a d-(. 
gauged S'^ model with N < d ioi the case = 3, = 2 , exte: 
Ref. [10] which were restricted to the N = d cases . (The work ^ 
of establishing topological lower bounds for the generic case, e: 
examples in two [11] and three [12, 13] dimensions respectively, 
scription used here by us coincides precisely with that used in I 
the establishing of a topological lower bound which did not fea 
which is carried out here to establish the stability of the soliton. 
are absent in the other prescription of gauging the Skyrme mc 
(Notice that we name the sigma models after the manifold in w 
their values rather than using the name of the symmetry group 1 
what is sometime called the 0(d+ 1) model in the literature will 
S"^ model.) 

The U{1) gauged SU(2) Skyrme model is described by the Li 

where the U(l) gauge covariant derivative is 

D,U = d,U + ieA,lQ,U], 
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where = eA^ and F^^, = eJ^^,,. 

The Lagrangian for the U(l) gauged Skyrme model can then 

^ = -Aoi^^i + M\D,r\'' - x2\D[,rD,](t>'\^ 

where the square brackets on the indices imply (total) antisymme 
Aq ^ = 4e^, Ai = F^/^ and A2 ^ = 8a^. The late Greek indices /i lal 
coordinates, while the early Greek indices a = 1,2 and the uppe 
A = 3, 4 label the fields = c/)^). 

The static Hamiltonian pertaining to the Lagrangian (4) is 

n = Aoi^5 + Xi\D,r? + X2\D[^rDJ](|)'\^ 

+i2Xo\d,Ao\^ + Al{Xr\r\' + 16A2[|r Pl^.^^P + 

where the indices i = a,3 label the space-like coordinates. 

To find the static solutions of the model, one would usualy 
grange equations which minimise the Hamiltonian (5), but bee 
potential Aq, one must solve the Euler Lagrange equations de 
grangian (4). We then look for static solutions, but, as for the Ju 
have to impose the proper asymptotic behaviour for the electric 
static solutions which are electrically charged (in the classical s 
where the fiux of the electric field is non zero). 

When the full equations of motion are written down, one fin^ 
there are static solutions for which Aq = 0, i.e. solutions for whid 
identically zero. For these solutions in the temporal gauge, the ( 
reduce to the equations obtained by minimising the Hamiltoniar 
solutions of unit Baryon charge of the U(l) gauged Skyrme mode! 
an electric field, for various values of the U{1) coupling constant 
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2 The topological charge and lower b( 



The definition of the topological charge is based on the criterion th 
Baryon number, namely the degree of the map. For the gauged t 
quantity must be gauge invariant as well. This requirement can b( 
satisfied by arranging the gauge invariant topological charge dei 
of the usual, gauge variant winding number density 

plus a total divergence whose surface integral vanishes due to the 
tions, such that the combined density is gauge invariant. In 3 dim( 
explicitly in Refs. [10, 13] for the 5*0 (3) gauged model, and foi 
interest, namely the 5*0(2) gauged model, the charge density ^ 
that of the SO (3) gauged model by contraction of the gauge gr( 
SO (2). It can also be arrived at directly. To state the definitio 
denote the gauge covariant counterpart of (6) by 

SO that using the notations (6) and (7) we have the definition of 
topological charge 

Q = Qo + di^i, 

= QG + ^£ijkFij{e^''(l)''Dk(l)^). 

In (8) the density Qi is the following gauge variant form 
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<^/xJ — '^c. c-afjCAB-^fiY J^i^V^ ^pY 

which is analogous to the corresponding quantity in the wort 
Wilczek [6]. This contrasts with the expression for the total diverg 
ical current in the work of d'Hoker and Farhi [7], where a different j 
is used leading to that quantity being equal to the local anomaly 

We now proceed to find a model whose Hamiltonian Tio is bou 
the topological charge density defined by (9). We will then show tl 
(5) is given by I-Lq plus certain positive definite terms. 

First of all, we reproduce the density qg, (7), in (9) by usii 
equality 

{KsD.r - £^Jk£'''''4DJ(|>'Dk(|>'(|>'f > 

where the two constants k,^ and k,2 have the dimensions of leng 
square, we get qg on the right hand side of 

To reproduce the other term in (9), ^£ijkFij{£^^(j)^dk4>^), 
inequality 

yielding 

With the special choice for the relative values of the constant 
sum of (13) and (15) yields the following 
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right hand side of (20) to Tio in (19) to end up with the Hamil 
gauged Skyrme model: 

which is nothing but the static Hamiltonian (5) in the temporal 
where ^ 

\ 9 / -1 \ \ 4 \ 4 

By virtue of (18), (21) is also bounded from below by 2KsKlg, n 
proportional to the topological charge density g. 

We thus see that Tio can be considered as a minimal (U{1) ; 
from now on, we will restrict our attention to the physically more 
and integrate it numerically to find its topologically stable finite 

The soliton solutions to the system (21) can only be found by 
order Euler-Lagrange equations, and not some first-order Bogomo 
saturating the inequalities (14) and (16) would not saturate the 
energy density functional V.. In this context we note that satur 
does indeed saturate the topological lower bound on the functio 
the inequality (18), and should it have turned out that the Bo^ 
arising from the saturation of (14) and (16) supported non-trivia 
would have been a very interesting system to consider. As it turn; 
Bogomol'nyi equations have only trivial solutions in exactly the 
case of the (ungauged) Skyrme model [1]. 

The energy for the static configuration, when the electric fi 
pressed as 

^(Ao, Ai, A2) = I d?x [Aoi^5 + Ai(A^")' + HD[^r-^ 
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Notice that for the usual Skyrme model we have 



(AiA2)i/2^.,(l,l) 



We will use (27) to compare the numerical solutions of the gai 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability cons 
in this Section apply only to the solutions with no electric field, 

3 The soliton and the charged f/(l) 

To find the static solutions, we have to look for the largest sym 
functional to be subjected to the variational principle, and look 
are invariant under that symmetry group. For the solutions in 
this is the static Hamiltonian (21), while for the solutions in t] 
Lagrangian (4). For our choice of gauge group the largest sym: 
group corresponding to an axial rotation in space-time and a gc 



on the gauge field. Defining the axial variables r = Jxl + X2 and 



of the coordinates Xi = (xa, ^3), = 1,2, the most general axialk 



[15] for the fields = {c/)'^, (j)^) (with a = 1, 2 and ^ = 3, 4), and 





sin / sin g n 



a 



sin / cos g 



A 



C2(r, Z) ^ 



A, 



ci(r, z) 



a 



r 



X 

r 



a 



r 
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^ 2 



f,' + + sin' f{gt + g^^) + ±^ sin^ / sin^ g 



+2A2 sin^ / 



The boundary conditions for the Skyrmion fields are the sar 
conditions for the hedgehog ansatz when expressed in the cyli 
where g = 7r/2 + arctan(z/r) and defining R = {z^ -\- r'^Y^'^^ f = 
and lim^^oo f(R) = 0- From this we can deduce that the functior 
at the origin and at infinity. For smothness along the z axis, each f 
Aa must satisfy the condition that the partial derivative with resp 
r = vanishes. The boundary conditions and the asymptotic be] 
are chosen so that the gauge fields A^^ are well defined and Aq look 
a coulomb field (i.e. with an electric charge but no magnetic chari 
that the total energy be finite. These conditions leads to the foil 

/(0,0) = TT /(r^oo,z^oo) =0 

^(r = 0, z < 0) =0 g{r = 0,z > 0) = tt 

a{r = 0,z) =1 CLr\oo = ^ 

ar{r = 0,z) =0 v4o(r — )• 00, z — )• 00) =VQ-\-q/r /. 

where R = (z'^ -\- r^)^/^ and where we have used the notation 
that the field g is undefined at the origin and the resulting discoi 
point is an artefact of the coordinate system used. The asymptc 
at infinity will be discussed in a later Section. To solve the equal 
is more conveniant to use the field ^0 rather than 6; this is why 
the boundary condition in terms of that field. On the other hand 
a simpler form when written in term of b, so we shall still use it 
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+ l^{9az - ag,)smfcosfsmg + -{af, - fa,] 
6 

jz = — [(f Or - gfr)a sin^ f sing 
r 

1 1 

+ -{gar - agr) sin / cos / sin ^ + -{afr - far) 

where we have denoted % = f etc. Note that the Baryon currer 
sensitive to the charge of the Nucleon, i.e. that the function 6(r, 
in them. 

We now turn to the equations to be solved, namely the Euler- 
arising from the variational principle applied to the Lagrangian (4 
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fz{fz9rr frrQz ~l~ frzQr frQrz ~l~ fzQr/'^) 
~^ frifrOzz fzzQr ~l~ frzQz fzQrz fzQz/'^^ 

+ sin ^/r2((a2 - b'^){{g'^ + g'^^) siri^ f - - f^) cos g 
({a^ - l)^){grr + Qzz - Qvlr) + Ig^iaa^ - bbr) + 2g^{aa, 



+4(ft2 - 6^) {frgr + fzOz) COS / sin / sin g^ 



0. 



In the case of the = gauge, equations (34) coincide with 
equations derived from the positive definite Hamiltonian densit 
that case, those equations also coincide with the Euler-Lagran, 
reduced two dimensional Hamiltonian obtained by subjecting (5^ 
by substituting the Ansatz (28)-(29) into it. This is expected dm 
sition of symmetry. In the Aq gauge, the Euler-Lagrange ec 
from the Lagrangian (4) which is not positive definite. Nonethe 
coincide with those arising from the reduced two dimensional L 
by subjecting the Lagrangian (4) to axial symmetry. (This happe 
Zee dyon [9].) 



3.1 Ao = : U{1) Skyrme soliton 

It is easy to see from (34) that there are solutions for which b 
As mentioned before, in that case, equation (34) can be obtainec 
Hamiltonian (30). Notice also that setting a = is not compatible 
conditions (Ai would not be well defined at the origin). We thus 
solution to carry a non-zero magnetic field. 

To show this we have to solve equations (34) numerically fo 
functions f{r,z),g(r,z) and a{r,z). 

We have restricted our numerical integrations to the case when 
n appearing in the axially symmetric Ansatz (28) is equal to 1, «. 
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that at Ao = 20, the energy of the gauged Skyrmion 1.22 (in un 
the lower bound 0.95 by 0.27. This is larger than 0.232, the exci 
Skyrmion energy over its lower bound. For smaller values of Aq Fi 
the excess of the energy of the gauged Skyrmion over its lower b 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, i.e.the terr 
in (30), is decreasing as Aq increases. Notice that we could have u 
Energy the sum of all the terms involving the gauge field functi( 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the ei 
Skyrmion in the r, z plane for Aq = 1. One sees clearly that the 
field is to make the Skyrmion elongated along the z axis. The m; 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we she 
of magnetic field using arrows to represent the magnetic field vec 
the grid. Notice that there is a vortex around the point r = 2,z 
field is thus generated by a current flowing on a ring centred aro 

In terms of the usual physical constants[3], we have Ag ^ = ^ 
A2 ^ = Sa^ where we use a instead of the traditional e for the S] 
avoid confusion with the electric charge. 

In our units, c = ^ = 1, we have e = {AiraY^'^ where a - 
structure constant. Choosing F^^ = 186MeV, we can find the vali 
that the energy of the neutron = 939MeV matches the ener^ 

F 2r7^ 

In Figure l.a, we can read the value of E(Xq^ 1,1) (given in units 
Aq = 2a^/e^. We now have to find the value of Aq for whicl 
2(2Ao)^/^eM„/F^. The intersection between the curve -^gri^X 
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5 10 15 20 

^0 



20 40 61 

^0 



Figure 1: a) Energy and topological bound of the gauged Skyrmi 
b) Ratio of the electromagnetic and the total energy as a functio 




Figure 2: a) Energy density for the gauged Skyrmion in the (r , 2 
A2 = 1). b) Energy density level curve for the gauged Skyrmion ( 
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Figure 3: Magnetic Field of the gauged Skyrmion (Aq = } 

E{Xq, 1, 1) in Figure l.a is located in the region where the energy 
the asymptotic value E{Xq, 1, 1) = 1.232. This means that a ^ Sti 
and that Aq ~ 1138. We can thus conclude that the effective imj 
term we have added to the Skyrme model is relatively small. 

This justifies the procedure used in [5] where the Skyrmioi 
an external magnetic field of a magnetic monopole. Indeed, a; 
generated by a Skyrmion is very small (for the parameters fitting 
the nucleons) the external field is much larger than the Skyrmioi 

It would be interesting to find the differences between the ele 
titles obtained from the ungauged model, as in [3], and our U(l) 
are not able to compute the solutions of the U{1) gauged model fc 
of the parameter Aq as this is too large, but as we now know t] 
the gauge field is very small, we can compute the latter pertur 
(ungauged) Hedgehog as an induced field. This enables the eval 
correction and the induced magnetic moment. This perturbative 
carried out in the Section 4. 

It can also be concluded that if the f/(l) gauged Skyrme mod( 
in [3] (by quantising the zero modes corresponding to the globa 
tion) but taking into account the electromagnetic field generate 
Skyrmion, the result would not differ very markedly from what m 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a fi 

3.2 ^0 7^ • charged U{1) Skyrmion 

We can now look for solutions with a non-zero electric charge b;" 
field b in our ansatz (29) does not vanish. To do this we follow 
as Julia and Zee [9] and require that the electric field be asympt 
^0 = + q/{r'^ + z'^y^'^ where Vq and q are two constants. In pra 
solutions for different values of Vq and evaluates q by computing t 
sought only those solutions, for which the electric fiux equals Att 
the electron. 

It is important to realise that in this case, equations (34) are 
imising the action and thus they do not minimise the Hamiltonic 

In our units, the charge of the electron is 0.303. In Figure 4 
as a function of Aq, as well as Vq as a function of Aq, so that q = 

One can see that, for a fixed value of Aq, the energy of the charg 
is smaller than the energy of the ungauged Skyrmion when Aq < 
larger than the energy of the uncharged gauged Skyrmion. If 
electrically charged and uncharged gauged Skyrmions were int( 
masses of the Proton and the Neutron nip and m^v, then on this \ 
we would have to conclude that (mp — rrifq) > which is not corre 
on the basis of its analogy with the dyon [9]. Clearly, to calculate 
correctly one would have to perform the collective coordinate qm 
[3], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is 
difficult to carry out the numerical computations accurately whe 

At this stage, it is worth saying a few words about the nur 
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4 Perturbation around the Hedgehog 



We have seen from the work of 3.1 above that the energy of the U ( 
for the physical values of the parameters, namely of the pion dec; 
U(l) coupling, does not differ significantly from that of the ung^ 
is therefore justified, for these values of the parameters, to treat 
perturbation to the Hedgehog in the same way as Klinhamer am 
the f/(l) field as a pertubation to the SU{2) sphaleron. We v 
compute the magnetic moment of the Neutron, as well as the (sn 
mass from that of the Hedgehog. 

The equations for the fields ((/)^, A^) are derived from the 
equation for will be of the form 

The method consists of setting the gauge field A^ to zero in the cu: 
in the equation for 0^) and calculating the resulting induced elect] 
The gauge field computed this way can then be interpeted as the i 
by the (ungauged) Hedgehog Skyrme field. With this perturbal 
possible to calculate the induced static magnetic potential (i = 
static electric potential Aq, which in this scheme vanishes and cai 
non-perturbatively. The reason simply is that restricting to tl 
Hedgehog, the zeroth component of the current jo at ^4^ = Vc 
turn in vanishing induced potential Aq according to (35). 

As a consequence the electric field will be identicaly zero, wt 
can derive the equation from the static Hamiltonian rather than fi 
Notice also that we could try to compute perturbatively a solutioi 
charges skyrmion by keeping in jo the terms proportional to A(^ 
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The shift in the energy of the Hedgehog due to the induced L 

= I d'x (AoF,,F,, + 4A J,) , 

in which = ji{0) is the current (37) for Ai = 0. (Note that all quj 
Ai = are denoted by Roman script, e.g. = jz(0), as well as the 
and curvature and F^j.) When equation (36) is satified for the 

= -Ao / d^x FijFij 



2 j d^x hi] 



which is, as expected, a strictly negative quantity. 

The current = (ja, ja) in (39) and (41) is given, for the Hed^ 
tion (38), by 

sin^F/Ai / ,2 ^\^F\\ 

j3 = 0. 

We now note that 5^ = 0, which means that equation (36), f 
following form in terms of the induced f/(l) connection = (A^ 

AqA Ac^ = —\oL ■ 

The solution is well known and can be expressed, using the obvio 
i{R) Eapxp in terms of (42), as 

1 /■ 1 

Aa(x) = -— --£a/3 / ] J-\ K^') d:^' , 

47rAo J |x — X I 
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tion F{x) (38). 

Comparing (48) with the usual Maxwell potential of a magne 

*w - ■ 

we find that fi = (0, 0, /i) is 

/i = 47r/ . 

We can evaluate the magnetic moment (49) and the energy cor: 
by the electromagnetic field by evaluating the integral (47) and 
we take the experimental values F^^ = ISQMeV and a = 7.2 we o 

/i = 0.01393/m = 0.43nm 
= -O.lkeV. 

The experimental value for the magnetic moment of the proton < 
respectivly fip = 0.0902/m = 2.79nm and fin = 0.0617/m = - 
other hand we take the values of the parameters derived in [3], 
a = 5.45, we have 

/i = 0.0468/m = 1.449nm 
= -0.32keV. 

The magnetic moment of a particle is strictly speaking a quan 
should be computed by quantising the SU (2) gauge degree of fret 
ertheless, we see that if the take the parameters derived in [3] th 
moment is of the correct order of magnitude. The sign is of coui 
the classical magnetic moment is a vector. We can thus conclude t 
a reasonable classical description of Nucleons and affords a metho 
electromagnetic field generated by the Skyrmion, classically. It ie 
see that a quantum property like the magnetic moment can be n 
by a purely classical procedure. 
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equations arising from the imposition of spherical symmetry we 
and hence technically much more amenable to the numerical ii 
case it was found that in addition to stable branches of solutic 
some unstable branches bifurcating from the former, the impc 
that there were indeed stable branches of solutions, characterise 
the) parameters of the model. It would be very interesting to ca 
corresponding to that of [20, 21], for the considerably more compk 
symmetric equations at hand. This however is technically beyo: 
present work. 

The energies of the gauged uncharged Skyrmions are smaller th 

1 /2 

usual ungauged Skyrmion. When the gauge coupling I/Aq goes 
gauged Skyrmion tends to the ungauged Skyrmion. We also not( 
the electrically charged Skyrmion is higher than the uncharged o 
of the dyon is higher than that of the monopole of the Georgi-Gl 
Perhaps the most interesting physical result of the present 
parameters in the model are fitted to reproduce physical quantiti 
the effect of the Maxwell term in the Skyrme Lagrangian is very sn 
for the physical value of the constant Aq = 1138, the energy of th 
Skyrmion differs little from that of the ungauged Skyrmion, as S( 
The gauged Skyrmion field itself is thus nearly radially symmetric 
field is not). 

Having found that the influence of the electromagnetic flel 
is small, we were pointed in the direction of treating the magn^ 
induced field perturbatively around the (ungauged) Hedgehog, 
thus, to compute the classical magnetic moment of the (uncharge 
Baryon charge, namely of the Neutron. The result is that the 
moment of the Skyrmion matches surprisingly well to the experii 
magnetic moments of the Nucleons. 
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1 Introduction 



For a long time now, much attention has been paid to the Skyrme [1] model in c 
dimensions. It is believed to be an effective theory for nucleons in the large A'' limii 

of QCD at low energies. The classical properties as well as the quantum properties 
of the model are in relatively good agreement with the observed properties of smal 
nuclei [2. 3. 1]. 

Gauged Skyrme models have been used in the past. The U{1) gauged model [2, 5 
was used to study the decay of nucleons in the vicinity of a monopole [5], while thf 
SU[2)l gauged model [7] was used to study the decay of nucleons when the Skyrmf 
model is coupled to the weak interactions [7]. The Skyrme model has also been usee 
to compute the quantum properties of the Skyrmion [3] where the gauge degrees o: 
freedom are quantised to compute the low energy eigenstates of a Skyrmion. These 
states were identified as the proton, the neutron and the delta. 

The aim of this work is to show that the Skyrme model can be coupled to a sel: 
contained electromagnetic field and that this f ^(1) gauged model has stable classica' 
solutions. In addition to these solitons with vanishing magnetic and electric flux, we 
show that this system supports solutions with nonvanishing electric flux which are 
analogous to the dyon solutions of the Georgi-Glashow model, just as the unchargec 
solitons are the analogues of the monoplolos [8] of that model. The electrically chargec 
lumps have larger energy, or mass, than the uncharged soliton, just like the Julia-Zee 
dyon [9] has larger energy, or is heavier, than the (electrically uncharged) monopole 
We shall refer to these lumps as charged U{V) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Maxwell gauged Skyrme 
model, the present work is also an example of a soliton in a cl-dimensional SO{N] 
gauged S"* model with N < d for the case d = 3, A' = 2 , extending the results o: 
Ref. [10] which were restricted to the N = d cases. (The work of Ref. [10] consists 
of establishing topological lower bounds for the generic case, encompassing earliei 
examples in two [11] and three [12, 13] dimensions respectively.) The gauging pre- 
scription used here by us coincides precisely with that used in Ref. [5] and permits 
the establishing of a topological lower bound which did not feature in Ref. [5] anc 
which is carried out here to establish the stability of the soliton. Such lower bounds 
are absent in the other prescription of gauging the Skyrme model as in Refs. [7] 
(Notice that we name the sigma models after the manifold in which the fields take 
their values rather than using the name of the symmetry group for the model. Thus 
what is sometime called the 0{d+l) model in the literature will be refered to as the 
S"* model.) 

The U{1) gauged SU{2) Skyrme model is described by the Lagrangian [5] 

where the (7(1) gauge covariant derivative is 

D^U = d^U + ieA^[Q,U], (i; 
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with [19] 

After performing the angular integrations we have 

A„(x) = -I(R) e^nxn , (46) 
with (R) given by the integral 

Finally, in the 1 region of interest, the induced U(l) potential is 
/ - 1 

= EapXp , with , I =^ ds , (48) 

H oAq Jo 

to be evaluated numerically using the numerically constructed hedgehog profile func- 
tion F(x) (38). 

Comparing (48) with the usual Maxwell potential of a magnetic dipole 

A(x) - 

^ ' ~ 4^i?3 ' 

we find that fj, = (0, 0, /i) is 

li = 4ni . (49) 

We can evaluate the magnetic moment (49) and the energy correction (41) induced 
by the electromagnetic field by evaluating the integral (47) and (48) numerically. If 
we take the experimental values F„ = 186MeV and a = 7.2 we obtain 

II = 0.01393/m = 0.43nm (50) 
AE = -O.lkeV. (51) 

The experimental value for the magnetic moment of the proton and the neutron are 

respectivly /tj, = 0.0902/7n = 2.79nni and /i„ = 0.0617/m = — 1.91nm. If on the 
other hand we take the values of the parameters derived in [3], F„ = 129MeV and 
a = 5.45, we have 

fj. = 0.0468/m = 1.449nm (52) 
AE = -0.32keV. (53) 

The magnetic moment of a particle is strictly speaking a quantum property and it 
should be computed by quantising the SU(2) gauge degree of freedom as in [3]. Nev- 
ertheless, we see that if the take the parameters derived in [3] the classical magnetic 
moment is of the correct order of magnitude. The sign is of course undetermined as 
the classical magnetic moment is a vector. We can thus conclude that our model offers 
a reasonable classical description of Nucleons and affords a method for computing the 
electromagnetic field generated by the Skyrmion, classically. It is quite surprising to 
see that a quantum property like the magnetic moment can be reasonably predicted 
by a purely classical procedure. 
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Skyrme coupling). For physical values of these parameters in the model, we find that 
the energy (mass) of the gauged Skyrmion does not differ significantly from that o: 
the ungauged charged-l Skyrmion, namely the familiar hedgehog [1]. This implie; 
that for these values of the physical parameters, the U{1) gauged Skyrmion can be 
regarded as a perturbation of the (ungauged) hedgehog, enabling the computatior 
of the magnetic moments of the gauged Skyrmion (i.e. the Neutron) and the shift 
of the energy of the gauged Skyrmion away from the energy of the hedgehog [1] 
perturbatively using the method employed by Klinkhamer and Manton [14] for thf 
sphaleron of the Weinberg-Salam model. 

In Section 2, we define the topological charge and establish the corresponding 
lower bound on the energy functional. In Section 3 we present the solutions whicf 
have no electric fields in the first subsection and electrically charged solutions ir 
the second subsection. The perturbation analysis of the gauged Skyrmion arounc 
the (ungauged) hedgehog is carried out in Section 4, and Section 5 is devoted to £ 
discussion of our results. 



2 The topological charge and lower bound 

The definition of the topological charge is based on the criterion that it be equal to thf 
Baryon number, namely the degree of the map. For the gauged theory however, this 
quantity must be gauge invariant as well. This requirement can be systematically [10 
satisfied by arranging the gauge invariant topological charge density to be the surr 
of the usual, gauge variant winding number density 

^0 = ei,ke'""%4>''di>t>''d,4>''4>'' , (6; 

plus a total divergence whose surface integral vanishes due to the finite energy condi- 
tions, such that the combined density is gauge invariant. In 3 dimensions, this is giver 
explicitly in Refs. [10, 13] for the 50(3) gauged model, and for the present case o: 
interest, namely the SO{2) gauged model, the charge density can be derived fron 
that of the 50(3) gauged model by contraction of the gauge group SO (3) down tc 
SO (2). It can also be arrived at directly. To state the definition of the charge, w( 
denote the gauge covariant counterpart of (6) by 

1>G = e,,,e''''"'D,d>"D,dfDk(lf (l>'' , (7; 

so that using the notations (6) and (7) we have the definition of the gauge invariant 
topological charge 

Q = Qo + dMi, (a; 

= QG + \sijuFij[e'''',l,''D„j,''). (9; 

In (8) the density Oj is the following gauge variant form 

Oi = 3eiije'**ylja,0'* 0^. (lO; 
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have used. To compute the static solutions, we have employed a relaxation method 
using finite difl'erences on a regular grid {dr = dz). This discretisation method is 
similar to the one employed in the numerical computation of the solutions of Skyrme 
models in 2+1 and 3+1 dimensions [16, 17, 18]. To compute the electrically charged 
solutions, we have imposed the boundary condition i(oo, oo) = V'o for different Vo 
and using a dichotomic method, wo have determined the values that give a solution 
with the same electric flux as the proton. Most of the simulations whore done on 
200 X 400 or 300 x 600 grids. By computing the same solution for various lattice sizes, 
we have empirically obtained the following relation for the expression of the energy 
of a solution : E = Eo + Kdr^ where £"o is the exact solution, dr = dz the lattice 
spacing and if is a constant which depends on Aq but takes values between 1 and 
0.5. We see thus that the energies we have obtained are accurate to within one or one 
half of a percent. This inaccuracy in the value of the energy is comparable to that 
of many other similar works on 2 dimensional systems [17, 18], and though it might 
look large, it does not affect any of the conclusions we have drawn. 

4 Perturbation around the Hedgehog 

We have seen from the work of 3.1 above that the energy of the U (1) gauged Skyrmion 
for the physical values of the parameters, namely of the pion decay constant and the 
f7(l) coupling, does not differ significantly from that of the ungauged Hedgehog. It 
is therefore justified, for these values of the parameters, to treat the 6'(1) field as a 
perturbation to the Hedgehog in the same way as Klinhamer and Manton [14] treat 
the U{1) field as a pertubation to the SU{2) sphaleron. We will then be able to 
compute the magnetic moment of the Neutron, as well as the (small) deviation of its 
mass from that of the Hedgehog. 

The equations for the fields (<t>'',Ai^) are derived from the Lagrangian (1). The 
equation for will be of the form 

Ao d,F^, = . (35) 

The method consists of setting the gauge field A^^ to zero in the current j,, in (35) (and 
in the equation for (j)") and calculating the resulting induced electromagnetic field A,,. 
The gauge field computed this way can then be interpeted as the U{1) field generated 
by the (ungauged) Hedgehog Skyrme field. With this perturbative procedure, it is 
possible to calculate the induced static magnetic potential Aj {i = 1, 2, 3), but not the 
static electric potential Aq, which in this scheme vanishes and can only be calculated 
non-perturbatively. The reason simply is that restricting to the use of the static 
Hedgehog, the zeroth component of the current jo at A^, = vanishes, resulting in 
turn in vanishing induced potential Aq according to (35). 

As a consequence the electric field will be identicaly zero, which implies that we 
can derive the equation from the static Hamiltonian rather than from the Lagrangian. 
Notice also that we could try to compute perturbatively a solution for the electrically 
charges skyrmion by keeping in jo the terms proportional to Aq, instead of setting 
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The right hand side of (18) is now proportional to the topological charge densitj 
Q defined by (9) so that the inequality (18) can be interpreted as the topologica 
inequality giving the lower bound on the energy density functional if we define tht 
latter to be the left hand side of (18), namely 

% = + ^m")' + 4iDiirD,-,4>''r + ^-f^(e^''rD,^y. (19; 

The Hamiltonian system (19) is almost the Hamiltonian of the gauged Skyrmf 
model (5) (remember that Aq = 0). It differs from the latter only in its last term 
Now we can use the identity 

(£'*^,>-^A</'-^)' = (A*^")" - [^((/-["A*^"')" + ((^'"A-^-^')"] (2o; 

and add the positive definite term |^5(0'"A<f''^')^ + (<?'"£'i<S'*')^j appearing on tht 
right hand side of (20) to 'Ho in (19) to end up with the Hamiltonian for the U{1] 
gauged Skyrme model: 

H = + i^iDi^f + 4{DiirDj]<l>''f > 2k^kIq (21] 

which is nothing but the static Hamiltonian (5) in the temporal gauge Ao = 0, anc 
where 

Ai = K^(l + ^), Ao = 4 \2 = 4. (22 

By virtue of (18), (21) is also bounded from below by 2Ki>;3^Q, namely by a numbei 
proportional to the topological charge density q. 

We thus see that tio can be considered as a minimal (U(l) gauged) model, bul 
from now on, we will restrict our attention to the physically more relevant model (21] 
and integrate it numerically to find its topologically stable finite energy solitons. 

The soliton solutions to the system (21) can only be found by solving the second- 
order Euler-Lagrange equations, and not some first-order Bogomol'nyi equations since 
saturating the inequalities (14) and (16) would not saturate the lower bound on the 
energy density functional "H. In this context we note that saturating (14) and (16] 
does indeed saturate the topological lower bound on the functional "Ho by virtue o: 
the inequality (18), and should it have turned out that the Bogomol'nyi equations 
arising from the saturation of (14) and (16) supported non-trivial solutions, then 'H^ 
would have been a very interesting system to consider. As it turns out however, these 
Bogomol'nyi equations have only trivial solutions in exactly the same way as in the 
case of the (ungauged) Skyrme model [1]. 

The energy for the static configuration, when the electric field vanishes, is ex- 
pressed as 

E{\o, Ai, A2) = j d^x [\oFfj + \(Di<lff + \^(Dyi<lfDj^<l>''f\ (23] 
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Figure 3; Magnetic Field of the gauged Skyrmion (Aq = Ai = A2 = 1). 



E{\o, 1, 1) in Figure l.a is located in the region where the energy is virtually equal to 
the asymptotic value £;(Ao, 1, 1) = 1.232. This means that a Ri 37t^1.232F,/M„ ^ 7.2 
and that Aq ~ 1138. Wo can thus conclude that the effective impact of the Maxwell 
term we have added to the Skyrmo model is relatively small. 

This justifies the procedure used in [5] where the Skyrmion was coupled with 
an external magnetic field of a magnetic monopole. Indeed, as the Maxwell field 
generated by a Skyrmion is very small (for the parameters fitting the actual mass of 
the nucleons) the external field is much larger than the Skyrmion's magnetic field. 

It would be interesting to find the differences between the electromagnetic quan- 
tities obtained from the ungauged model, as in [3], and our U{1) gauged model. We 
are not able to compute the solutions of the U{1) gauged model for the physical value 
of the parameter Aq as this is too large, but as we now know that the influence of 
the gauge field is very small, we can compute the latter perturbatively around the 
(ungauged) Hedgehog as an induced field. This enables the evalution of the energy 
correction and the induced magnetic moment. This perturbative analysis will be be 
carried out in the Section 4. 

It can also be concluded that if the U{1) gauged Skyrme model were quantised as 
in [3] (by quantising the zero modes corresponding to the global gauge transforma- 
tion) but taking into account the electromagnetic field generated classically by the 
Skyrmion, the result would not differ very markedly from what was obtained in [3]. 



13 



with n" = {sin ruj), cos n(j)) in terms of the azimuthal angle (l> and = Xa/r. n in n" 
is the vorticity, which for the Nucleons of interest to us here, equals one, n = I. Thf 
functions a, 6, Ci, C2, / and g both depend on r and z. 

Our Ansatz (29) for the U{1) field consists of decomposing the latter in the mosi 
general tensor basis possible. We will find out below, when we compute the Euler- 
Lagrange equations, that the functions Ci and Ca vanish identically. Anticipating this 
we suppress them henceforth. In its final form this Ansatz agrees with that used ir 
Ref. [15], the latter being arrived at by specialising the Rebbi-Rossi Ansatz for thf 
axially symmetric SO (3) field. 

The static Hamiltonian, i. e.the Too component of the energy mometum tensor Tj,„ 
is then given by 

H= f{^[a'r + al + ibr-lr + bl^ 

+2\2 sin^ / [Urg. - Szgrf + sin^ 9[^(/,^ + ft + {gl + gf) sin^ /)]] }rdrd^. 

(30; 

The boundary conditions for the Skyrmion fields are the same as the boundar} 
conditions for the hedgehog ansatz when expressed in the cylindrical coordinates 
where i; = 7r/2 + arctan(2/r) and defining R = (z^ + r'^fl^, J = J{R) with /(O) = 7, 
and lima^oo f{R) ~ 0. Prom this we can deduce that the function / has a fixed valu< 
at the origin and at infinity. For smothness along the z axis, each field, that is /, g anc 
Aa must satisfy the condition that the partial derivative with respect to of the field a1 
r = Q vanishes. The boundary conditions and the asymptotic behaviours for a and / 
are chosen so that the gauge fields are well defined and looks asymptoticaly likt 
a coulomb field (i.e. with an electric charge but no magnetic charge). We also requirt 
that the total energy be finite. These conditions leads to the following constraints; 



/(r- -> 00, z ^ 00) =0 /r(r = 0,z) =0 

g{r = f),z >0) = JT flRlffi^oo = 

'^r 1 00 '^sloo ^ 

Aa{r ^ 00, z ^ 00) =Vo + q/r Aa{r = Q,z) =0 

(3i; 

where R = [z"^ + r'^yl'^ and where we have used the notation |^ = etc. Notf 
that the field g is undefined at the origin and the resulting discontinuity of g at thai 
point is an artefact of the coordinate system used. The asymptotic behaviour of A^ 
at infinity will be discussed in a later Section. To solve the equations numerically, ii 
is more conveniant to use the field ^0 rather than b; this is why we have expressec 
the boundary condition in terms of that field. On the other hand, the equations take 
a simpler form when written in term of 6, so we shall still use it below. 



/(0,0) =7r 

S(r = 0,z<0) =0 

a(r = Q,z) =1 

arir = Q,z) =0 
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unit Baryon number. 

Using (30), we have found numerically that i?(l,l,l) = 247r^l.01 whereas for 

those values of Aq, Ai, A2 the lower bound for the energy given by (26) is 24jr^0.555. 
In Figure 1, we present the total energy for the gauged Skyrmion as a function of 
Ao, together with the lower bound given by (26). Note that the asymptotic value 
of i?(Ao,l,l) is 24-^1.232 as Ao — > 00. As a comparison, the energy (27) for the 
ungaugcd Skyrmion is Egi^(l, 1) = 247r^l.232 with a lower bound sot at 247r^. We see 
that E{Xo = 00, 1, 1) = E,k{l, 1) which means that as Ao — > 00, the gauge coupling 
1/Ay^ goes to zero and the gauged Skyrmion becomes in this limit the ungauged 
Skyrmion. 

It is interesting to note that the energy of gauged Skyrmion is smaller than the 
energy of the ungauged Skyrmion, as expected, but that on the other hand, the 
amount by which the energy of the gauged Skyrmion exceeds its topological lower 
bound is larger than the excess of the energy of the ungauged Skyrmion above its 
respective topological lower bound. For example we can clearly see from Figure l.a 
that at Ao = 20, the energy of the gauged Skyrmion 1.22 (in units of 247r^) exceeds 
the lower bound 0.95 by 0.27. This is larger than 0.232, the excess of the ungauged 
Skyrmion energy over its lower bound. For smaller values of Ao Figure l.a shows that 
the excess of the energy of the gauged Skyrmion over its lower bound is even larger, 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, i.e. the term proportional to Ao 
in (30), is decreasing as Ao increases. Notice that we could have used for the Maxwell 
Energy the sum of all the terms involving the gauge field functions a and b in (30), 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the energy density of the 
Skyrmion in the r, z plane for Ao = 1. One sees clearly that the effect of the gauged 
field is to make the Skyrmion elongated along the z axis. The magnetic field vectors 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we show the configuration 
of magnetic field using arrows to represent the magnetic field vector at each point on 
the grid. Notice that there is a vortex around the point r = 2, z = 0. The magnetic 
field is thus generated by a current flowing on a ring centred around the z axis. 

In terms of the usual physical constants[3], we have Ao^ = 4e^, Ai = and 
A2 ^ = 8a^ where we use a instead of the traditional e for the Skyrme coefficient to 
avoid confusion with the electric charge. 

In our units, c = ft = 1, we have e = (4jra')'^^ where a = 1/137 is the fine 
structure constant. Choosing = ISdMeV, we can find the value for o by requiring 
that the energy of the neutron M„ = Q39MeV matches the energy of the Skyrmion: 

M„ = |^E(^,1,1). 

In Figure l.a, we can read the value of E{Xo, 1, 1) (given in units of 24:Tt^ MeV) with 
Ao — 2o^/e^. We now have to find the value of Ao for which 24jr^ii(Ao, 1, 1) = 
2(2Ao)^^^eM„/F„. The intersection between the curve ^^'^ (2Ao)^^^ and the curve 
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Substituting the ansatz (28) (29) into these equations of motion leads to 
arr - - ^i°'(9) + '-t^n + n + (9i + 9l) sin^(/)]] = 

Kr -'f + ^+ b.. - 6sin^(/) Sin^(s) [Al + ^[/2 + /2 + (^2 + ^2) ,i„2(^)]j ^ 

A/ -(a? + 9K^sin'fl)™/cos/+ 

+4ff sin / [{gr{f,:,gr - frfjzz + fzQrz " frzQz " fz/r(Jz) 

+gz(frTgz - fzQrr + fr9rz " frzQr + fr/rgz))sm f + (frQz " OrfzY COS / 

+ Sin ((a^ - b')if? + - 2(5,^ + gl) sin^ /) cos / sin g 

+ ((a^ - 6^)(/rr + fzz - fr/r) + 2/r(aGr - bhr) + 2f^(aa^ - 66,)) sin/sins 
+2(a2 - b^)(f,g, + /,<?,) sin /cos 9)] = 0. 

^g +'^{frgr + fzgz) eotf - sing cos g+ 

+4^7 [fzifzgrr - frrgz + /r^Sr " /rSrz + fzgr/r) 
+ ir(Jrgzz - fzzgr + /rzSz " /zSrz " /zSz/r) 

+ smg/r'{ia' - b')i{g^ + gl) sin^ f - - /|) cosfl 

((a^ - 6^)(g„ + fc - gr/r) + 2gr{aar - bbr) + 2g^{aa^ - 66, ) ) sin V sin g 
+4(a^ - 6^)(/,p, + /^5,) cos /sin /sins)] = 0- 

(31 

In the case of the ^0 = gauge, equations (34) coincide with the Euler-Lagrangx 
equations derived from the positive definite Hamiltonian density (5). Moreover ir 
that case, those equations also coincide with the Euler-Lagrange equations of thf 
reduced two dimensional Hamiltonian obtained by subjecting (5) to axial symmetr} 
by substituting the Ansatz (28)-(29) into it. This is expected due to the strict impo- 
sition of symmetry. In the ^0 7^ gauge, the Euler-Lagrange equations are derivec 
from the Lagrangian (4) which is not positive definite. Nonetheless these equations 
coincide with those arising from the reduced two dimensional Lagrangian obtainec 
by subjecting the Lagrangian (4) to axial symmetry. (This happens also for the Julii 
Zee dyon [9].) 

3.1 = : [/(I) Skyrme soliton 

It is easy to see from (34) that there are solutions for which 6 = (i.e. = 0) 
As mentioned before, in that case, equation (34) can be obtained by minimising thf 
Hamiltonian (30). Notice also that setting a = is not compatible with our boundar) 
conditions (Ai would not be well defined at the origin). We thus expect our gaugec 
solution to carry a non-zero magnetic field. 

To show this we have to solve equations (34) numerically for the non-vanishing 
functions j(r,z),g(r,z) and a(r,z). 

We have restricted our numerical integrations to the case where the vortex numbei 
n appearing in the axially symmetric Ansatz (28) is equal to 1, i.e. our soliton carries 
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Now the volume integral (with the appropriate normalisation of 127r^) of Qo given 
by (6) is the degree of the map, or, the Baryon number. It is straightforward to verify 
that when the Ansatz (28) is substituted in go and the volume integral is computed 
subject to the boundary conditions given above, the result will equal the integer n 
defined in (28). Thus, the Baryon number of the field configuration (28) equals the 
vortex number n. In what follows, wc will restrict ourselves to unit Baryon number, 
n = 1, i.e.to the Nucleons. 

Before we proceed to substitute the Anstaz (28), (29) into the field equations, we 
calculate the Baryon current (11) for the field configurations (28), (29) described by 
the solutions we seek. We express the space-like part of this current jj in the radial 
direction flowing out of the normal to the surface of the cylinder which we denote by 
jr, and in the z direction which we denote by j^. The result is 

jr = -[(/S2 -s/JasinVsinp (32) 
r 

+ 2(5«z - ag,) sin / cos / sin 5 + -(d/^ - /aj cosg] 
6 

jz = — [(/(/r - p/r)a sin V sing (33) 
r 

+ 2(3"^ -d(7r)sin/cos/sin(jr+ -(d/r -/ar)cossr] , 

where we have denoted = / etc. Note that the Baryon current (32), (33) are not 
sensitive to the charge of the Nucleon, i.e. that the function b{r,z) does not feature 
in them. 

We now turn to the equations to be solved, namely the Euler-Lagrange equations 
arising from the variational principle applied to the Lagrangian (4), in the static limit. 
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Figure 1: a) Energy and topological bound of the gauged Skyrmion in units of 247r^ 
b) Ratio of the electromagnetic and the total energy as a function of Ao . 




Figure 2: a) Energy density for the gauged Skyrmion in the {r, z) plane (Ao = Ai 
A2 = 1). b) Energy density level curve for the gauged Skyrmion (Aq = Ai = A2 = 1) 
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and performing the dilation x — > ax, Af^^ ^ a ^^j,, we get 

E{Xo, Ai, A^) = / d'x + a\,(Dirr + (24) 

If we choose a = (^)^^^ then we have 

i?(Ao,Ai,A2) = (AiA2)^/^E(^,l,l), (25) 

from which we see that we can set Ai = A2 = 1 without any loss of generality. By 
virtue of (21) and (25), we can finally state 




Notice that for the usual Skyrme model we have 

E,,(Ai,A2) = J dx'lXMr? + Hdlirdj-^tj,")'] 

= (AiA2)i/^i?.,(l,l) (27) 

> 2(Xi\,y/^Jdx'i>o- 

We will use (27) to compare the numerical solutions of the gauged Skyrme model 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability considerations discussed 
in this Section apply only to the solutions with no electric field, ie.with Ao = 0. 

3 The soliton and the charged U{1) Skyrmion 

To find the static solutions, we have to look for the largest symmetry group of the 
functional to be subjected to the variational principle, and look for solutions which 
are invariant under that symmetry group. For the solutions in the Aq = gauge 
this is the static Hamiltonian (21), while for the solutions in the Aq ^ it is the 
Lagrangian (4). For our choice of gauge group the largest symmetry is the S0{2) 
group corresponding to an axial rotation in space-time and a gauge transformation 
on the gauge field. Defining the axial variables r — y'xf + and and z — X3 in terms 
of the coordinates Xi = (zq, Z3), a = 1,2, the most general axially symmetric Ansatz 
[15] for the fields tj)" = {(j)", (j)^) (with a = 1, 2 and yl = 3, 4), and, A = (A^, A3), is 

0° = sin / sin g n" , <t>^ = sin / cos g , 0'* = cos / , (28) 

a(r, z)-n , , C2{r,z) ^ ci{r,z) b{r, z) 
Ac = £a0 H Xa As = , Ao = , (29) 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a function of Ao 
3.2 Aq^O : charged U{1) Skyrmion 

We can now look for solutions with a non-zero electric charge by requiring that tht 
field b in our ansatz (29) does not vanish. To do this we follow the same procedure 
as Julia and Zee [9] and require that the electric field be asymptotically of the forir 
^0 = Vo + q/{T^ + z^)'-/^ where Vq and g are two constants. In practice, one compute! 
solutions for different values of Vq and evaluates q by computing the electric flux. W( 
sought only those solutions, for which the electric flux equals 47r times the charge o; 
the electron. 

It is important to realise that in this case, equations (34) are obtained after min- 
imising the action and thus they do not minimise the Hamiltonian (30). 

In our units, the charge of the electron is 0.303. In Figure 4 we show the energj 
as a function of Ao, as well as I'o as a function of Ao, so that q = 0.303. 

One can see that, for a fixed value of Aq, the energy of the charged gauged Skyrmior 
is smaller than the energy of the ungauged Skyrmion when Ao < 7 but it is always 
larger than the energy of the uncharged gauged Skyrmion. If the energies of tht 
electrically charged and uncharged gauged Skyrmions were interpreted as the tht 
masses of the Proton and the Neutron mp and mjv, then on this purely classical leve 
we would have to conclude that (mp — mN) > which is not correct. This is expectec 
on the basis of its analogy with the dyon [9] . Clearly, to calculate this mass difference 
correctly one would have to perform the collective coordinate quantisation as in Ref 
[3], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is unfortunately verj 
difficult to carry out the numerical computations accurately when Ao is very large. 

At this stage, it is worth saying a few words about the numerical methods we 
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The flux of Oj vanishes, as can deduced by anticipating the flnite energy conditions 
to be stated later. 

Note that the 3-volume integral of go in (8) is the degree of the map for the 
ungauged system namely the Baryon number. 

Identifying g, (9), with the naught component of the Baryon current, is 
defined by 

/ = e^'''''e,,,Ard,,4>%^'4>' + 3e''-'''^ d^iA,e^"(l>''d„0'') (H) 
= e,,,aD.r D^'D.^'?''' - ^ e"""" F,,(eAB<p'' D,^-^). (12) 

The 4-divergence of (11) receives a contribution only from its first term, which 
being locally a total divergence implies that the 3-volume integral of j" is a conserved 
quantity. Alternatively we consider the 4-divergence of (12), 

di^j" = &£''"'"'£^D£ABD^rD,<l>fD,<j>''D,<l>'' (13) 

which is analogous to the corresponding quantity in the work of Goldstone and 
Wilczek [6]. This contrasts with the expression for the total divergence of the topolog- 
ical current in the work of d'Hoker and Far hi [7], where a different gauging prescription 
is used leading to that quantity being equal to the local anomaly. 

We now proceed to find a model whose Hamiltonian "Ho is bounded from below by 
the topological charge density defined by (9). We will then show that the Hamiltonian 
(5) is given by Ho plus certain positive definite terms. 

First of all, we reproduce the density Qq, (7), in (9) by using the following in- 
equality 

(KaA^" - e„k£'''""'i4Di<t>''D,,4,'<t>^f > (14) 

where the two constants K3 and K2 have the dimensions of length. Expanding the 
square, we get Qg on the right hand side of 

njiDir? + 4iD[.rDj]<l>''Y > 2k3kIqg. (15) 

To reproduce the other term in (9), \eijkFij{e^^ (j)^ dk<t>^), we use the following 
inequality 

(nlFij - ^K4£ijkS^'',j)''Dk<j)^Y > (16) 

yielding 

+ 4\{e^"'P''D,<j}^)' > 4^^4e„tF„{e^''<j}''D,<t,^). (17) 

With the special choice for the relative values of the constants Sk^k^ = the 
sum of (13) and (15) yields the following 

4F^, + 4(D,rr + 4{Dv^-D,]^''f + ^(.^«0^A0^)^ > 2^,4e- (18) 
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^0 to 0, and impose the condition that the electric field is assymptoticaly like thai 
of the proton. This pertubation method would not make much sense though as out 
would expect the electric field to be quite large close to the skyrmion. 

The relevant energy functional is (23), and the resulting equation arising from th( 
variation of the gauge field Ai is 

Ao djFij = ji (36~ 

We are concerned here with the case where Ai = (37) and the chiral field 0" = 
, (j)^ , 4>*) in (37) describes the Hedgehog, i.e. 

0° = sin F{R) x" , <t>^ = sin F(R) , c^" = cos F(R) , (38] 



where now x" = x" /R, with R — ^/r^~+^. By virtue of equation (36) the current 
(37), given by (38) and Ai = 0, will now induce a (small) U(l) field A;, with curvature 
^ij — diA.j 5jAj. 

The shift in the energy of the Hedgehog due to the induced U{1) field A; is 

AE = I ifx (XoFijFij + 4Ai ji) , (39; 

in which = ji(0) is the current (37) for Ai = 0. (Note that all quantities evaluated ai 
Ai = are denoted by Roman script, e.g. j; = ji{0), as well as the induced connectior 
and curvature A; and Fy.) When equation (36) is satified for the induced U{1) field 

AE = -Ao j d^x FjjFjj (40; 

= 2 1 (fxAiii, (4i; 

which is, as expected, a strictly negative quantity. 

The current j; = (ja, ja) in (39) and (41) is given, for the Hedghog field configura- 
tion (38), by 

sin'F/Ai , / ,2 sin'F\\ , 
}a = ^^(y + 2^2(^ (42. 

j.3 = 0. (43 

We now note that 9; j; = 0, which means that equation (36), for Ai = 0, takes th( 
following form in terms of the induced U{1) connection Aj = (A„, 0) 

AoA A„ = -j„ . (44; 

The solution is well known and can be expressed, using the obvious notation ja(x) = 
j(i?) Sa/iX/i in terms of (42), as 

1 f 1 

Aa(x) = -— — Ea/J / , j(i?') dx' , (45; 

47rAo J |x — X 
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where = d^Ay — d^Ai^ and where the charge matrix of the quarks is expressed 
as Q = + as). This differs from the covariant derivative of Ref. [5] only in the 
unimportant matter of the sign of i in (1), which we have chosen for consistency of 
the convention used in Ref. [10]. 

In what follows it will be more convenient [10] to parametrise the Skyrme field as 
an valued field (j)" = (ip", a = 1, 2, = 3, 4 subject to the constraint [(^"[^ = 1. 
The two fields U and (j) are related to each other via the following expression 

U = 0V°, U-^ = U^ = (^"f" (2) 

where r" = (ia",ia^,l) and f° = (— ia", — ia^, 1), in terms of the Pauli matrices 
(a^a^a='). 

The gauge covariant derivative now can be re-expressed as 

Z)„r = d.^j," + A.e^^^j,", D.ct,^ = d,<p^. (3) 

where = eA^ and Fp„ = ej^.^. 

The Lagrangian for the U{1) gauged Skyrme model can then be written as 

■C = -\^Fl + Ai|D^,^"|^ - A2|D[^<^"D„]<^T (4) 

where the square brackets on the indices imply (total) antisymmetrisation and where 
Aq ' = 4e^, Ai = F^/<i and Aj ' = <io?. The late Greek indices /i label the Minkowskian 
coordinates, while the early Greek indices a = 1, 2 and the upper case Latin indices 
^ = 3, 4 label the fields (j)" = (j)^). 

The static Hamiltonian pertaining to the Lagrangian (4) is 

H = Aoi^5 + AilA'^l' + A2|D[i0"Dj]0T 

(5) 

+z2Ao|s,Ap + Aliwr? + i6A2[|ri'|a.0-^r + iia.d^-'l')!']}, 

where the indices i = a, 3 label the space-like coordinates. 

To find the static solutions of the model, one would usualy solve the Euler La- 
grange equations which minimise the Hamiltonian (5), but because of the electric 
potential Aq, one must solve the Euler Lagrange equations derived from the La- 
grangian (4). We then look for static solutions, but, as for the Julia-Zee dyon [9], we 
have to impose the proper asymptotic behaviour for the electric potential to obtain 
static solutions which are electrically charged (in the classical sense, i.e. solutions 
where the fiux of the electric field is non zero). 

When the full equations of motion are written down, one finds as expected that 
there are static solutions for which = 0, i.e. solutions for which the electric field is 
identically zero. For these solutions in the temporal gauge, the equations of motion 
reduce to the equations obtained by minimising the Hamiltonian (5). Wo study the 
solutions of unit Baryon charge of the U (1) gauged Skyrme model with and whithout 
an electric field, for various values of the (7(1) coupling constant (or equivalently the 
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5 Summary and discussion 



We have shown that the SU(2) Skyrme model gauged with U{1) has two types o: 
finite energy static solutions, electrically uncharged and charged respectively. Bott 
of these solutions are axially symmetric and carry no magnetic charge but support j 
magnetic field shaped like a torus centred around the axis of symmetry, albeit resulting 
in zero magnetic fiux. The uncharged solutions, like the ungauged Skyrmion, have e 
topological lower bound. The electrically charged solutions are the analogues of th( 
Julia-Zee dyons [9] of the Georgi-Glashow model. 

Concerning the stability of the electrically neutral solution, which is expected to b( 
stable by virtue of the lower bound on the energy, we have not made any quantitativf 
effort to test it. Wo expect however that the solitons of this gauged Skyrme model ar( 
stable, or that at least they have stable branches for all values of the parameters in the 
model. This expectation is based on our knowledge of the corresponding situatior 
when the Skyrme model is gauged instead with SO(3) [20, 21], in which case th( 
equations arising from the imposition of spherical symmetry were one dimensional 
and hence technically much more amenable to the numerical integration. In thai 
case it was found that in addition to stable branches of solutions, there were alsc 
some unstable branches bifurcating from the former, the important matter beinf 
that there were indeed stable branches of solutions, characterised by the (ranges o: 
the) parameters of the model. It would be very interesting to carry out the analysii 
corresponding to that of [20, 21], for the considerably more complex case of the axiallj 
symmetric equations at hand. This however is technically beyond the scope of thf 
present work. 

The energies of the gauged uncharged Skyrmions are smaller than the energy of the 
usual ungauged Skyrmion. When the gauge coupling l/A^^ goes to 0, the unchargec 
gauged Skyrmion tends to the ungauged Skyrmion. We also note that the energy o; 
the electrically charged Skyrmion is higher than the uncharged one, just as the mass 
of the dyon is higher than that of the monopole of the Georgi-Glashow model. 

Perhaps the most interesting physical result of the present work is that whoi 
parameters in the model are fitted to reproduce physical quantities, it turns out thai 
the effect of the Maxwell term in the Skyrme Lagrangian is very small. This is because 
for the physical value of the constant Aq = 1138, the energy of the gauged unchargec 
Skyrmion differs little from that of the ungauged Skyrmion, as seen from Figure lb 
The gauged Skyrmion field itself is thus nearly radially symmetric (though the gauge 
field is not). 

Having found that the influence of the electromagnetic held on the Skyrmior 
is small, we were pointed in the direction of treating the magnetic potential as ar 
induced field perturbatively around the (ungauged) Hedgehog. We have been able 
thus, to compute the classical magnetic moment of the (uncharged) Skyrmion of uni 
Baryon charge, namely of the Neutron. The result is that the classical magnetie 
moment of the Skyrmion matches surprisingly well to the experimental values of the 
magnetic moments of the Nucleons. 



18 



Static solutions in the U{1) gauged Skyrme model 



B.M.A.G. Piette* 

Department of Mathematical Sciences 
University of Durham 
Durham DHl 3LE, UK 

and 

D. H. Tchrakiant 

Department of Mathematical Physics 
National University of Ireland, Maynooth 
Maynooth, Ireland 
and 

School of Theoretical Physics 
Dublin Institute for Advanced Studies 
10 Burlington Road 
Dublin 4, Ireland 

Abstract 

We use a prescription to gauge the SU{2) Skyrme model with a U{1) field, 
characterised by a conserved Baryonic current. This model reverts to the usual 
Skyrme model in the limit of the gauge coupling constant vanishing. We show 
that there exist axially symmetric static solutions with zero magnetic charge, 
which can be electrically cither charged or uncharged. The energies of the 
(uncharged) gauged Skyrmions arc less than the energy of the (usual) ungaugcd 
Skyrmion. For physical values of the parameters the impact of the U{1) field 
is very small, so that it can be treated as a perturbation to the (ungauged) 
spherically symmetric Hedgehog. This allows the perturbative calculation of 
the magnetic moment. 



*e-mail address: B.M.A.G. Piette@durham.ac.uk 
^e-mail address: tigran@thphys.may.ie 



1 



[20] Y. Brihaye and D. H. Tchrakian, Nonlinearity 11 891 (1998) 

[21] Y. Brihaye, B. Kleihaus and D. H. Tchrakian, J. Math. Phys 40 1136 (1999) 



20 



[20] Y. Brihaye and D. H. Tchrakian, Nonlinearity 11 891 (1998) 

[21] Y. Brihaye, B. Kleihaus and D. H. Tchrakian, J. Math. Phys 40 1136 (1999) 



20 



Static solutions in the U{1) gauged Skyrme model 



B.M.A.G. Piette* 

Department of Mathematical Sciences 
University of Durham 
Durham DHl 3LE, UK 

and 

D. H. Tchrakiant 

Department of Mathematical Physics 
National University of Ireland, Maynooth 
Maynooth, Ireland 
and 

School of Theoretical Physics 
Dublin Institute for Advanced Studies 
10 Burlington Road 
Dublin 4, Ireland 



Abstract 

We use a prescription to gauge the SU{2) Skyrme model with a U{1) field, 
characterised by a conserved Baryonic current. This model reverts to the usual 
Skyrme model in the limit of the gauge coupling constant vanishing. Wc show 
that there exist axially symmetric static solutions with zero magnetic charge, 
which can be electrically cither charged or uncharged. The energies of the 
(uncharged) gauged Skyrmions are less than the energy of the (usual) ungauged 
Skyrmion. For physical values of the parameters the impact of the U{1) field 
is very small, so that it can be treated as a perturbation to the (ungauged) 
spherically symmetric Hedgehog. This allows the perturbative calculation of 
the magnetic moment. 



*e-mail address: B.M.A.G. Piette@durham.ac.uk 
^e-mail address: tigran@thphys.may.ie 



1 



1 Introduction 



For a long time now, much attention has been paid to the Skyrme [1] model in 3 
dimensions. It is believed to be an effective theory for nucleons in the large N limit 
of QCD at low energies. The classical properties as well as the quantum properties 
of the model are in relatively good agreement with the observed properties of small 
nuclei [2, 3, 4]. 

Gauged Skyrme models have been used in the past. The U{1) gauged model [2, 5] 
was used to study the decay of nucleons in the vicinity of a monopole [5] , while the 
SU(2)l gauged model [7] was used to study the decay of nucleons when the Skyrme 
model is coupled to the weak interactions [7]. The Skyrme model has also been used 
to compute the quantum properties of the Skyrmion [3] where the gauge degrees of 
freedom are quantised to compute the low energy eigenstates of a Skyrmion. These 
states were identified as the proton, the neutron and the delta. 

The aim of this work is to show that the Skyrme model can be coupled to a self 
contained electromagnetic field and that this U{1) gauged model has stable classical 
solutions. In addition to these solitons with vanishing magnetic and electric flux, we 
show that this system supports solutions with nonvanishing electric flux which are 
analogous to the dyon solutions of the Georgi-Glashow model, just as the uncharged 
solitons are the analogues of the monoploles [8] of that model. The electrically charged 
lumps have larger energy, or mass, than the uncharged soliton, just like the .luIia-Zcc 
dyon [9] has larger energy, or is heavier, than the (electrically uncharged) monopole. 
Wc shall refer to these lumps as charged U{1) Skyrmions. 

In addition to its intrinsic interest as a soliton in the Maxwell gauged Skyrme 
model, the present work is also an example of a soliton in a li-dimensional SO{N) 
gauged S'' model with N < d ion the case d = 3, N = 2 , extending the results of 
Ref. [10] which were restricted to the N = d cases. (The work of Ref. [10] consists 
of establishing topological lower bounds for the generic case, encompassing earlier 
examples in two [11] and three [12, 13] dimensions respectively.) The gauging pre- 
scription used hero by us coincides precisely with that used in Ref. [5] and permits 
the establishing of a topological lower bound which did not feature in Ref. [5] and 
which is carried out here to establish the stability of the soliton. Such lower bounds 
are absent in the other prescription of gauging the Skyrme model as in Refs. [7]. 
(Notice thai wc name the sigma models after the manifold in which the fields take 
their values rather than using the name of the symmetry group for the model. Thus 
what is sometime called the 0{d+ 1) model in the literature will be refered to as the 
S" model.) 

The U(l) gauged SU{2) Skyrme model is described by the Lagrangian [5] 
C = ^Tr{D,UD,U^) - ^,Tr{[(D,U)UK (D.up^f - 

where the U{1) gauge covariant derivative is 

Di,U = di,U + ieAiAQ,U], (1) 
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Ao to 0, and impose the condition that the electric field is assymptoticaly like that 
of the proton. This pertubation method would not make much sense though as one 
would expect the electric field to be quite large close to the skyrmion. 

The relevant energy functional is (23), and the resulting equation arising from the 
variation of the gauge field At is 

\d,F,,=j, (36) 
ji = -e"^ (iAi(4^A<?" + As \{dj\<p-^)D,ifDjd^ + 2(j^D[i<?-"D;]0-*D;(i6'*]) . (37) 

We are concerned here with the case where Ai = Q (37) and the chiral field (jf = 
{(j)" , (j)^ , 4>^) in (37) describes the Hedgehog, i.e. 

(j," = sin F{R) x" , 4? = sin F{R) i^ c^" = cos F{R) , (38) 



where now x"' = x''/R, with R — ^/r^~+^. By virtue of equation (36) the current 
(37), given by (38) and Ai = 0, will now induce a (small) U(l) field Aj, with curvature 

The shift in the energy of the Hedgehog due to the induced U{1) field Aj is 

AE = y ifx {XoFijFij + 4Ai ji) , (39) 

in which = ji(0) is the current (37) for Ai — 0. (Note that all quantities evaluated at 
Ai = are denoted by Roman script, e.g. j; = ji(0), as well as the induced connection 
and curvature A; and Fy.) When equation (36) is satified for the induced (7(1) field, 

AE = -Ao j d\ FjjFjj (40) 

= 2 j (fxAiii, (41) 

which is, as expected, a strictly negative quantity. 

The current j; = (ja, ja) in (39) and (41) is given, for the Hedghog field configura- 
tion (38), by 

sin"F/Ai , / ,2 sin'FAA , 

j3 = 0. (43) 

We now note that 9; j; = 0, which means that equation (36), for Ai = 0, takes the 
following form in terms of the induced U(l) connection Aj = (Aq, 0) 

AoA A„ = -j„ . (44) 

The solution is well known and can be expressed, using the obvious notation ja(x) = 
j(i?) Sa/iX/i in terms of (42), as 

1 r 1 

Aa(x) = --^£ai} / 1 71 i{R') f^x' , (45) 

47rAo J |x — X I 
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The flux of Oj vanishes, as can deduced by anticipating the finite energy conditions 
to be stated later. 

Note that the 3-volume integral of Qq in (8) is the degree of the map for tht 
ungauged system namely the Baryon number. 

Identifying q, (9), with the naught component j" of the Baryon current, is 
defined by 

f = el'"''' E,,,A4>''d,<t>'d,<t>'<t>' + 7, dM,s^''<l>''d.<^^) (H! 

The 4-divergence of (11) receives a contribution only from its first term, whicf 
being locally a total divergence implies that the 3-volume integral of j° is a conservec 
quantity. Alternatively we consider the 4-divergence of (12), 

d^f = 66''-''"'e,f,6ABDi,<j>''D,<j>CDf,<j>^D,<j>'' [13] 

which is analogous to the corresponding quantity in the work of Goldstone anc 
Wilczek [6]. This contrasts with the expression for the total divergence of the topolog- 
ical current in the work of d'Hoker and Farhi [7], where a different gauging prescriptior 
is used leading to that quantity being equal to the local anomaly. 

We now ijrocoed to find a model whose Hamiltonian "Ho is bounded from below bj 
the topological charge density defined by (9). We will then show that the Hamiltoniar 
(5) is given by Ho plus certain positive definite terms. 

First of all, we reproduce the density Qg, (7), in (9) by using the following in- 
equality 

(KsDir - eijte''''"'4Dj(l>''D;,4,'4,''f > (14; 

where the two constants K3 and K2 have the dimensions of length. Expanding thf 
square, we get Qg on the right hand side of 

i4(Di4,''f + 4{D,^4,"Dj]4/'f > 2K.n,KlQG- (15; 

To reproduce the other term in (9), \£ijkFij{£'^^4'^dk4>^), we use the following 
inequality 

(K^Fi^ - iK4eijfce'**0*Dfc(>'*)' > (16; 

yielding 

<Ff, + Kl^-{e^''^''D,<t>y > 4K,e,,,F„{e^W'D,^^). (If 

With the special choice for the relative values of the constants 3K3K2 = tht 
sum of (13) and (15) yields the following 

4f^, + 4(Dirr + «^(i?[,0"-D,]0')^ + ^(e^«0«A0^)^ > 2«3«^e. (is; 
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The right hand side of (18) is now proportional to the topological charge density 
Q defined by (9) so that the inequality (18) can be interpreted as the topological 
inequality giving the lower bound on the energy density functional if we define the 
latter to be the left hand side of (18), namely 

no = 4Ffi + ^m-f + 4iDi.rDj]^r + ^(e-^^^^A^-^)^. (19) 

The Hamiltonian system (19) is almost the Hamiltonian of the gauged Skyrme 
model (5) (remember that Aq = 0). It differs from the latter only in its last term. 
Now we can use the identity 

(6^''<j>''D,<j>^f = (A<^°)' - [^(<?i'"A/')' + (0'"A<^'*')'] (20) 

and add the positive definite term ^^j^ |j (<?>[" A^''')^ + (^'"A^'^')^! appearing on the 
right hand side of (20) to Ho in (19) to end up with the Hamiltonian for the (7(1) 
gauged Skyrme model; 

H = KlFf^ + /c?(A0")' + 4{D[irDj^<l,''f > 2k,4q (21) 

which is nothing but the static Hamiltonian (5) in the temporal gauge = 0, and 
where 

By virtue of (18), (21) is also bounded from below by 2k3K.I(), namely by a number 
proportional to the topological charge density q. 

We thus see that Tio can be considered as a minimal {U{1) gauged) model, but 
from now on, we will restrict our attention to the physically more relevant model (21) 
and integrate it numerically to find its topologically stable finite energy solitons. 

The soliton solutions to the system (21) can only be found by solving the second- 
order Euler-Lagrange equations, and not some first-order Bogomol'nyi equations since 
saturating the inequalities (14) and (16) would not saturate the lower bound on the 
energy density functional "H. In this context we note that saturating (14) and (16) 
does indeed saturate the topological lower bound on the functional "Ho by virtue of 
the inequality (18), and should it have turned out that the Bogomol'nyi equations 
arising from the saturation of (14) and (16) supported non-trivial solutions, then Ho 
would have been a very interesting system to consider. As it turns out however, these 
Bogomol'nyi equations have only trivial solutions in exactly the same way as in the 
case of the (ungauged) Skyrme model [1]. 

The energy for the static configuration, when the electric field vanishes, is ex- 
pressed as 

EiXo, Ai, As) = / d'x [Ao4 + \i{D,rf + X^iDy,^" D .^tj^'f] (23) 
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have used. To compute the static solutions, we have employed a relaxation method 
using finite differences on a regular grid {dr = dz). This discretisation method is 
similar to the one employed in the numerical computation of the solutions of Skyrme 
models in 2-1-1 and 3-1-1 dimensions [16, 17, 18]. To compute the electrically charged 
solutions, we have imposed the boundary condition 6(oo, oo) — Vq for different V'o 
and using a dichotomic method, we have determined the values that give a solution 
with the same electric flux as the proton. Most of the simulations where done on 
200 X 400 or 300 x 600 grids. By computing the same solution for various lattice sizes, 
we have empirically obtained the following relation for the expression of the energy 
of a solution : E = Eo + Kdr'^ where Eo is the exact solution, dr = dz the lattice 
spacing and K \s a constant which depends on Aq but takes values between 1 and 
0.5. We see thus that the energies we have obtained are accurate to within one or one 
half of a percent. This inaccuracy in the value of the energy is comparable to that 
of many other similar works on 2 dimensional systems [17, 18], and though it might 
look large, it does not affect any of the conclusions we have drawn. 

4 Perturbation around the Hedgehog 

We have seen from the work of 3.1 above that the energy of the U{1) gauged Skyrmion 
for the physical values of the parameters, namely of the pion decay constant and the 
(7(1) coupling, does not differ significantly from that of the ungauged Hedgehog. It 
is therefore justified, for these values of the parameters, to treat the (7(1) field as a 
perturbation to the Hedgehog in the same way as Klinhamer and Manton [14] treat 
the (7(1) field as a pertubation to the SU{2) sphaleron. We will then be able to 
compute the magnetic moment of the Neutron, as well as the (small) deviation of its 
mass from that of the Hedgehog. 

The equations for the fields ((^°,^j,) are derived from the Lagrangian (1). The 
equation for A/^ will be of the form 

Ao d,F^, = . (35) 

The method consists of setting the gauge field A^j to zero in the current jf, in (35) (and 
in the equation for 0°) and calculating the resulting induced electromagnetic field A^^. 
The gauge field computed this way can then be interpeted as the U (1) field generated 
by the (ungauged) Hedgehog Skyrme field. With this perturbative procedure, it is 
possible to calculate the induced static magnetic potential A; {i = 1, 2, 3), but not the 
static electric potential Aq, which in this scheme vanishes and can only be calculated 
non-perturbatively. The reason simply is that restricting to the use of the static 
Hedgehog, the zeroth component of the current jo at A^ = vanishes, resulting in 
turn in vanishing induced potential Aq according to (35). 

As a consequence the electric field will be idcnticaly zero, which implies that wo 
can derive the equation from the static Hamiltonian rather than from the Lagrangian. 
Notice also that we could try to compute perturbatively a solution for the electrically 
charges skyrmion by keeping in jo the terms proportional to Ao, instead of setting 
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Figure 1: a) Energy and topological bound of the gauged Skyrmion in units of 247r^. 
b) Ratio of the electromagnetic and the total energy as a function of Aq . 




Figure 2: a) Energy density for the gauged Skyrmion in the (r, z) plane (Aq = Ai = 
A2 = 1). b) Energy density level curve for the gauged Skyrmion (Aq = Ai = A2 = 1). 
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Now the volume integral (with the appropriate normalisation of 127r^) of go giver 
by (6) is the degree of the map, or, the Baryon number. It is straightforward to verifj 
that when the Ansatz (28) is substituted in Qq and the volume integral is computec 
subject to the boundary conditions given above, the result will equal the integer r 
defined in (28). Thus, the Baryon number of the field configuration (28) equals tli< 
vortex number n. In what follows, wc will restrict ourselves to unit Baryon number 
n = 1, ie.to the Nucleons. 

Before we proceed to substitute the Anstaz (28), (29) into the field equations, wf 
calculate the Baryon current (11) for the field configurations (28), (29) described b) 
the solutions we seek. We express the space-like part of this current j; in the radia 
direction flowing out of the normal to the surface of the cylinder which we denote bj 
jr, and in the z direction which we denote by j^. The result is 

if = -[(/S2 - 3/2)0 sin V sin p (32; 
r 

+ ^{9az - ag^) sin / cos / sin 5 + -(a/^ - fa^) cosg] 

= --[{fgr - gfr)asm^ fsing {33 
r 

+ 2^9(^r - agr) sin /cos/ sin (jr + -{afr - far)cosg] , 

where we have denoted = / etc. Note that the Baryon current (32), (33) are no 
sensitive to the charge of the Nucleon, i.e. that the function b{r, z) does not featurt 
in them. 

We now turn to the equations to be solved, namely the Euler-Lagrange equations 
arising from the variational principle applied to the Lagrangian (4), in the static limit 
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Substituting the ansatz (28) (29) into these equations of motion leads to 
a.. -'f - a.. - asm'if) sin^(fi) [i^^ + ^[/,^ + /| + (g^^ + gl) sin^(/)]] = 
Kr -'f + ^ + b..- bsin'if) sin^(fi) [Ai^ + + /I + (p^2 + gl) si„^(/)]] = 

A/ -(9,^ + 9,^ + ^sin^9)sin/cos/+ 

+4|^ sin/[((7r(/zzSr - /ri7z2 + /zSrz - frzOz - fz/rOz) 

+9z0rr9z - fzOrr + frOrz - frzOr + fr/rgz))sill f + (/rfi^ - flr/z)^ COS / 

+ Sin gy {(a' - 6^) {f^ + fl - 2(gl + gl) sin^ /) cos / sin g 

+ ((a^ - 6^)(/„ + fzz - fr/r) + 2fr{aar - bbr) + 2f^(aa^ - 66^)) sin / sing 
+2(a2 - 6')(/,fi. + /,fc)sin/cosfl)] = 0. 

Afi +2(/r5i, + /25^)cot/ - ^!^sin5Cosg+ 

|^/£(/25rr frrQz ~l~ frzQr frQrz ~l~ fz9r/^^ 

+ fr{frOzz - fzzOr + /rzi7z - fzOrz - fzOz/r) 

+ sin g/r^(^(a'^ - b'^){(gl + gl) siv? f - - /|) cos^i 

((a^ - b'^)(grT + gzz - gr/r) + 2gr(aar - bb^) + 2g^(aa^ - bb ^)) sin^ f sing 
+4(a^ - 62)(/rgr + /2g2) cos/sin/sin s)] = 0. 

(34) 

In the case of the Aq = gauge, equations (34) coincide with the Euler-Lagrange 
equations derived from the positive definite Hamiltonian density (5). Moreover in 
that case, those equations also coincide with the Euler-Lagrange equations of the 
reduced two dimensional Hamiltonian obtained by subjecting (5) to axial symmetry 
by substituting the Ansatz (28)-(29) into it. This is expected due to the strict impo- 
sition of symmetry. In the j4o 7^ gauge, the Euler-Lagrange equations are derived 
from the Lagrangian (4) which is not positive definite. Nonetheless these equations 
coincide with those arising from the reduced two dimensional Lagrangian obtained 
by subjecting the Lagrangian (4) to axial symmetry. (This happens also for the Julia 
Zee dyon [9].) 

3.1 ^0 = 0: [/(I) Skyrme soliton 

It is easy to see from (34) that there are solutions for which 6 = {i.e. Ao = 0). 
As mentioned before, in that case, equation (34) can be obtained by minimising the 
Hamiltonian (30). Notice also that setting a = is not compatible with our boundary 
conditions (Ai would not be well defined at the origin). We thus expect our gauged 
solution to carry a non-zero magnetic field. 

To show this wo have to solve equations (34) numerically for the non-vanishing 
functions f{r, z), g{r, z) and a{r, z). 

We have restricted our numerical integrations to the case where the vortex number 
n appearing in the axially symmetric Ansatz (28) is equal to 1, ie.our soliton carries 
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unit Baryon number. 

Using (30), we have found numerically that £(1,1,1) = 247r^l.01 whereas for 
those values of Aq, Ai, A2 the lower bound for the energy given by (26) is 247r^0.555. 
In Figure 1, we present the total energy for the gauged Skyrmion as a function of 
Ao, together with the lower bound given by (26). Note that the asymptotic value 
of i?(Ao,l,l) is 247r^l.232 as Aq — s- 00. As a comparison, the energy (27) for the 
ungauged Skyrmion is £,*(!, 1) = 247r^l.232 with a lower bound set at 247r^. We see 
that E(Xo = 00, 1, 1) = Eski^, 1) which means that as Ao 00, the gauge coupling 
I/aJ^^ goes to zero and the gauged Skyrmion becomes in this limit the ungauged 
Skyrmion. 

It is interesting to note that the energy of gauged Skyrmion is smaller than the 
energy of the ungauged Skyrmion, as expected, but that on the other hand, the 
amount by which the energy of the gauged Skyrmion exceeds its topological lower 
bound is larger llian the excess of the energy of the ungauged Skyrmion above its 
respective topological lower bound. For example we can clearly see from Figure l.a 
that at Ao = 20, the energy of the gauged Skyrmion 1.22 (in units of 247r^) exceeds 
the lower bound 0.95 by 0.27. This is larger than 0.232, the excess of the ungauged 
Skyrmion energy over its lower bound. For smaller values of Ao Figure l.a shows that 
the excess of the energy of the gauged Skyrmion over its lower bound is even larger, 
hence this is a general feature. 

In Figure l.b, we also see that the Maxwell Energy, ie.the term proportional to Aq 
in (30), is decreasing as Aq increases. Notice that we could have used for the Maxwell 
Energy the sum of all the terms involving the gauge field functions a and b in (30), 
but this would lead to a figure similar to Figure l.b. 

In Figure 2, we show the profile and the level curve for the energy density of the 
Skyrmion in the r, z plane for Aq = 1. One sees clearly that the effect of the gauged 
field is to make the Skyrmion elongated along the z axis. The magnetic field vectors 
of the Skyrmion are parallel to the r, z plane. In Figure 3, we show the configuration 
of magnetic field using arrows to represent the magnetic field vector at each point on 
the grid. Notice that there is a vortex around the point r = 2, 2; = 0. The magnetic 
field is thus generated by a current flowing on a ring centred around the z axis. 

In terms of the usual physical constants[3], we have Aq ' = 4e^, \\ = F^/?: and 
A2 ^ = 8a^ where we use a instead of the traditional e for the Skyrme coefficient to 
avoid confusion with the electric charge. 

In our units, c = fi = 1, wc have e = (47ra)'^^ where a — 1/137 is the fine 
structure constant. Choosing = 186Mey, we can find the value for a by requiring 
that the energy of the neutron M„ = 939Mey matches the energy of the Skyrmion: 

F 2(1^ 
M„ = ^£;(^,1,1). 

In Figure La, we can read the value of E{\, 1, 1) (given in units of 24jr^i'VfeV') with 
Ao = 2(!?jf?. We now have to find the value of Ao for which 247r^ii(Ao, 1, 1) = 
2(2Ao)^^^eM„/F^. The intersection between the curve (2Ao)^^^ and the curve 
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with n° = (sin cosn0) in terms of the azimuthal angle (j) and = Xa/r. n in n° 
is the vorticity, which for the Nucleons of interest to us here, equals one, n = 1. The 
functions a, 6, Ci, C2, / and g both depend on r and z. 

Our Ansatz (29) for the U{1) field consists of decomposing the latter in the most 
general tensor basis possible. We will find out below, when we compute the Euler- 
Lagrango equations, that the functions Ci and C2 vanish identically. Anticipating this, 
we suppress them henceforth. In its final form this Ansatz agrees with that used in 
Ref. [15], the latter being arrived at by specialising the Rebbi-Rossi Ansatz for the 
axially symmetric S0(3) field. 

The static Hamiltonian, ie.the Too component of the energy mometum tensor Tj,„, 
is then given by 

H= jl^^^al + al + {K-^r + bl\ 

+2A2 sin^ / [(As. - hgrf + sin^ 9[^(/,^ + fl + {gl + gl) sin^ /)]] ^drdz. 

(30) 

The boundary conditions for the Skyrmion fields are the same as the boundary 
conditions for the hedgehog ansatz when expressed in the cylindrical coordinates 
where g — 7r/2 + arctan(.z/r) and defining R — {z^ + r^Y^^, f = f{R) w'ith /(O) — tt 
and limg^oo f{R) ~ 0. Prom this we can deduce that the function / has a fixed value 
at the origin and at infinity. For smothness along the z axis, each field, that is /, g and 
Aa must satisfy the condition that the partial derivative with respect to of the field at 
r = vanishes. The boundary conditions and the asymptotic behaviours for a and b 
are chosen so that the gauge fields are well defined and Aq looks asymptoticaly like 
a coulomb field (i.e. with an electric charge but no magnetic charge). We also require 
that the total energy be finite. These conditions leads to the following constraints: 



/(0,0) =7r f(r^oo,z^oo) =0 /r(r = 0,z) =0 

g{r = 0,z<0) =0 g{r = 0, z > 0) = tt teU^oc =0 

a{r = 0,z) =1 Orloo =0 cizloo =0 

ar{r = 0,z) =0 Ao(r ^ (x,z (x) =Vo + q/r Ao(r = 0,z) =0 

(31) 

where R = (z^ + r^Y^^ and where we have used the notation |^ = etc. Note 
that the field g is undefined at the origin and the resulting discontinuity of g at that 
point is an artefact of the coordinate system used. The asymptotic behaviour of Aq 
at infinity will be discussed in a later Section. To solve the equations numerically, it 
is more conveniant to use the field ^0 rather than b; this is why we have expressed 
the boundary condition in terms of that field. On the other hand, the equations take 
a simpler form when written in term of 6, so we shall still use it below. 
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Figure 3: Magnetic Field of the gauged Skyrmion (Ao = Ai = A2 = 1). 



E{\o, 1, 1) in Figure l.a is located in the region where the energy is virtually equal tc 
the asymptotic value E{\o, 1, 1) = 1.232. This means that Ri 37r^l.232F,j/Af„ ^ 7.^ 
and that Ao ~ 1138. We can thus conclude that the cfl'cctivc impact of the Maxwcl 
term we have added to the Skyrmo model is relatively small. 

This justifies the procedure used in [5] where the Skyrmion was coupled witl 
an external magnetic field of a magnetic monopole. Indeed, as the Maxwell fielc 
generated by a Skyrmion is very small (for the parameters fitting the actual mass o: 
the nucleons) the external field is much larger than the Skyrmion's magnetic field. 

It would be interesting to find the differences between the electromagnetic quan- 
tities obtained from the ungauged model, as in [3]. and our U{1) gauged model. We 
are not able to compute the solutions of the U{1) gauged model for the physical valuf 
of the parameter Ao as this is too large, but as we now know that the infiuence o: 
the gauge field is very small, we can compute the latter perturbatively around thf 
(ungauged) Hedgehog as an induced field. This enables the evalution of the energj 
correction and the induced magnetic moment. This perturbative analysis will be bt 
carried out in the Section 4. 

It can also be concluded that if the U{1) gauged Skyrme model were quantised as 
in [3] (by quantising the zero modes corresponding to the global gauge transforma- 
tion) but taking into account the electromagnetic field generated classically by th( 
Skyrmion, the result would not differ very markedly from what was obtained in [3]. 
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Figure 4: a) Energy of the charged Skyrmion b) Vq as a function of Aq 
3.2 Ao^O : charged [7(1) Skyrmion 

We can now look for solutions with a non-zero electric charge by requiring that the 
field b in our ansatz (29) does not vanish. To do this we follow the same procedure 
as Julia and Zee [9] and require that the electric field be asymptotically of the form 

^0 = Vo + q/{r'^ + z'^)'-^^ where V'o and q are two constants. In practice, one computes 
solutions for different values of Vo and evaluates q by computing the electric flux. We 
sought only those solutions, for which the electric flux equals Air times the charge of 
the electron. 

It is important to realise that in this case, equations (34) are obtained after min- 
imising the action and thus they do not minimise the Hamiltonian (30). 

In our units, the charge of the electron is 0.303. In Figure 4 we show the energy 
as a function of Aq, as well as Vq as a function of Aq, so that q = 0.303. 

One can see that, for a fixed value of Aq, the energy of the charged gauged Skyrmion 
is smaller than the energy of the ungauged Skyrmion when Aq < 7 but it is always 
larger than the energy of the uncharged gauged Skyrmion. If the energies of the 
electrically charged and uncharged gauged Skyrmions were interpreted as the the 
masses of the Proton and the Neutron mp and mn, then on this purely classical level 
we would have to conclude that {mp — mfi) > which is not correct. This is expected 
on the basis of its analogy with the dyon [9]. Clearly, to calculate this mass difference 
correctly one would have to perform the collective coordinate quantisation as in Ref. 
[3], which we do not do here. 

The energy of the charged Skyrmion increases with Aq. It is unfortunately very 
difhcult to carry out the numerical computations accurately when Aq is very large. 

At this stage, it is worth saying a few words about the numerical methods we 



14 



and performing the dilation x — > ax, ^ a ^yl^,, we get 

E(Ao, Ai, A2) = / d'x + aX,(D,rr + ^(C[i0"C^]0')']. (24) 

If we choose a = (f^)^'^ then we have 

E(X„,Xr,X2) = (XrX2y''E(^,l,l), (25) 

from which we see that we can set Ai = A2 = 1 without any loss of generality. By 
virtue of (21) and (25), we can finally state 



EiXo,XuX,)>2i-^Y/' [d'xs. (26) 



A1A2 

Notice that for the usual Skyrme model we have 

Esk{XuX2) = J dx'ix.idir? + x.,(diirdj]^''y] 

= (X,X2Y''E.,(1,1) (27) 

> 2(AiA2)i/Vdx3ft. 

We will use (27) to compare the numerical solutions of the gauged Skyrme model 
with the solutions of the (ungauged) Skyrme model. 

We would like to point out that the topological stability considerations discussed 
in this Section apply only to the solutions with no electric field, i.e. with Aq = 0. 



3 The soliton and the charged U{1) Skyrmion 

To find the static solutions, wc have to look for the largest symmetry group of the 
functional to be subjected to the variational principle, and look for solutions which 
are invariant under that symmetry group. For the solutions in the Ao = gauge 
this is the static Hamiltonian (21), while for the solutions in the ^0 7^ it is the 
Lagrangian (4). For our choice of gauge group the largest symmetry is the SO(2) 
group corresponding to an axial rotation in space-time and a gauge transformation 
on the gauge field. Defining the axial variables r = \l ~c\ + x\ and and z = x-^ in terms 
of the coordinates Xi = [x^jXs), a = 1,2, the most general axially symmetric Ansatz 
[15] for the fields 0" = (0", 0^) (with a = 1, 2 and ^ = 3, 4), and, A = [A^, A3), is 

= sin / sin p n° , (j)^ = sin / cos g , <j>^ = cos / , (28) 
a(r,z)—n Calr^z) , Ci(r,z) , b(r,z) , . 



7 



Skyrme coupling). For physical values of these parameters in the model, we find that 
the energy (mass) of the gauged Skyrmion does not differ significantly from that of 
the ungauged charged-l Skyrmion, namely the familiar hedgehog [1]. This implies 
that for these values of the physical parameters, the U{1) gauged Skyrmion can be 
regarded as a perturbation of the (ungauged) hedgehog, enabling the computation 
of the magnetic moments of the gauged Skyrmion (i.e. the Neutron) and the shift 
of the energy of the gauged Skyrmion away from the energy of the hedgehog [1], 
perturbatively using the method employed by Klinkhamer and Manton [14] for the 
sphaleron of the Weinberg-Salam model. 

In Section 2, we define the topological charge and establish the corresponding 
lower bound on the energy functional. In Section 3 we present the solutions which 
have no electric fields in the first subsection and electrically charged solutions in 
the second subsection. The perturbation analysis of the gauged Skyrmion around 
the (ungauged) hedgehog is carried out in Section 4, and Section 5 is devoted to a 
discussion of our results. 



2 The topological charge and lower bound 

The definition of the topological charge is based on the criterion that it be equal to the 
Baryon number, namely the degree of the map. For the gauged theory however, this 
quantity must be gauge invariant as well. This requirement can be systematically [10] 
satisfied by arranging the gauge invariant topological charge density to be the sum 
of the usual, gauge variant winding number density 

^0 = e,,te'"'"'dj''d,^''dt^'<t>'' , (6) 

plus a total divergence whose surface integral vanishes due to the finite energy condi- 
tions, such that the combined density is gauge invariant. In 3 dimensions, this is given 
explicitly in Refs. [10, 13] for the 50(3) gauged model, and for the present case of 
interest, namely the SO{2) gauged model, the charge density can be derived from 
that of the SO(3) gauged model by contraction of the gauge group S0{3) down to 
S0{2). It can also be arrived at directly. To state the definition of the charge, we 
denote the gauge covariant counterpart of (6) by 

&G = e,,^""^D,d^'D,dfD,(l/^d^ , (7) 

so that using the notations (6) and (7) we have the definition of the gauge invariant 
topological charge 

Q = ft + SjOi, (8) 
= QG + \eijkFii{e'''',j,''D^,j,''). (9) 

In (8) the density 0; is the following gauge variant form 

Oi = aejj.e^^^ja,^'* 0^. (10) 
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(47; 



with [19] 

= E E ^|fr yii\S'A') Yi'\9,4>) . 
After performing the angular integrations wo have 

Aa(x) = -I{R) SapXp , (46 

with I(t){R) given by the integral 

Finally, in the 1 region of interest, the induced U{1) potential is 

/ - 1 

Aa(x) = eap , with , / = — / s^iis) ds , (48; 
XI oAq Jo 

to be evaluated numerically using the numerically constructed hedgehog profile func- 
tion F(x) (38). 

Comparing (48) with the usual Maxwell potential of a magnetic dipole 
we find that /i = (0, 0, fi) is 

n = Avi . (49; 

We can evaluate the magnetic moment (49) and the energy correction (41) inducec 
by the electromagnetic field by evaluating the integral (47) and (48) numerically. I: 
we take the experimental values F„ = 186MeV and a = 7.2 we obtain 

= 0.01393/m = 0.43nm (50; 
AE = -O.lA^eV. (51; 

The experimental value for the magnetic moment of the proton and the neutron art 
respectivly //p = 0.0902/m = 2.79nm and = 0.0617/m = — 1.91nr«. If on the 
other hand we take the values of the parameters derived in [3], F„ = 129MeV anc 
a = 5.45, we have 

fj. = 0.0468/m = 1.449nm (52; 
AE = -0.32keV. [53] 

The magnetic moment of a particle is strictly speaking a quantum property and it 
should be computed by quantising the SU(2) gauge degree of freedom as in [3]. Nev- 
ertheless, we see that if the take the parameters derived in [3] the classical magnetic 
moment is of the correct order of magnitude. The sign is of course undetermined as 
the classical magnetic moment is a vector. We can thus conclude that our model offers 
a reasonable classical description of Nucleons and affords a method for computing th( 
electromagnetic field generated by the Skyrmion, classically. It is quite surprising tc 
see that a quantum property like the magnetic moment can be reasonably predictec 
by a purely classical procedure. 
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5 Summary and discussion 



We have shown that the SU{2) Skyrme model gauged with U{1) has two types of 
finite energy static solutions, electrically uncharged and charged respectively. Both 
of these solutions are axially symmetric and carry no magnetic charge but support a 
magnetic field shaped like a torus centred around the axis of symmetry, albeit resulting 
in zero magnetic flux. The uncharged solutions, like the ungauged Skyrmion, have a 
topological lower bound. The electrically charged solutions are the analogues of the 
Julia-Zee dyons [9] of the Georgi-Glashow model. 

Concerning the stability of the electrically neutral solution, which is expected to be 
stable by virtue of the lower bound on the energy, we have not made any quantitative 
effort to test it. We expect however that the solitons of this gauged Skyrme model are 
stable, or that at least they have stable branches for all values of the parameters in the 
model. This expectation is based on our knowledge of the corresponding situation 
when the Skyrme model is gauged instead with 50(3) [20, 21], in which case the 
equations arising from the imposition of spherical symmetry were one dimensional, 
and hence technically much more amenable to the numerical integration. In that 
case it was found that in addition to stable branches of solutions, there were also 
some unstable branches bifurcating from the former, the important matter being 
that there were indeed stable branches of solutions, characterised by the (ranges of 
the) parameters of the model. It would be very interesting to carry out the analysis 
corresponding to that of [20, 21], for the considerably more complex case of the axially 
symmetric equations at hand. This however is technically beyond the scope of the 
present work. 

The energies of the gauged uncharged Skyrmions are smaller than the energy of the 
usual ungauged Skyrmion. When the gauge coupling I/aJ'^ goes to 0, the uncharged 
gauged Skyrmion tends to the ungauged Skyrmion. We also note that the energy of 
the electrically charged Skyrmion is higher than the uncharged one, just as the mass 
of the dyon is higher than that of the monopole of the Georgi-Glashow model. 

Perhaps the most interesting physical result of the present work is that when 
parameters in the model are fitted to reproduce physical quantities, it turns out that 
the effect of the Maxwell term in the Skyrme Lagrangian is very small. This is because 
for the physical value of the constant Aq = 1138, the energy of the gauged uncharged 
Skyrmion differs little from that of the ungauged Skyrmion, as seen from Figure lb. 
The gauged Skyrmion field itself is thus nearly radially symmetric (though the gauge 
field is not). 

Having found that the influence of the electromagnetic fleld on the Skyrmion 
is small, we were pointed in the direction of treating the magnetic potential as an 
induced field perturbatively around the (ungauged) Hedgehog. We have been able 
thus, to compute the classical magnetic moment of the (uncharged) Skyrmion of unit 
Baryon charge, namely of the Neutron. The result is that the classical magnetic 
moment of the Skyrmion matches surprisingly well to the experimental values of the 
magnetic moments of the Nucleons. 
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where = d/j^A^ — d^Ai^ and where the charge matrix of the quarks is expressed 
as Q = i(|l + 1T3). This differs from the covariant derivative of Ref [5] only in the 
unimportant matter of the sign of i in (1), which we have chosen for consistency of 
the convention used in Ref. [10]. 

In what follows it will be more convenient [10] to parametrise the Skyrme field as 
an valued field 0" = (0", 4>^), a = l,2,A = 3,i subject to the constraint [(^"l^ = 1. 
The two fields U and (j) are related to each other via the following expression 

U = 4,''t'', u-^ = (yt = 0V (2) 

where t" = (ia",ia^,l) and f" = (—ia°,—ia^,l), in terms of the Pauli matrices 

The gauge covariant derivative now can be re-expressed as 

D,r = diA" + ^^.e""/, Cf.^-* = a^.^-*. (3) 

where = eAi^ and F^^j, = ej-^^. 

The Lagrangian for the U{1) gauged Skyrme model can then be written as 

C = -AoF^l + \i\D^r? - MDy^rD^-^ct-'? (4) 

where the square brackets on the indices imply (total) antisymmetrisation and where 
Aq ' = 4e^, Ai = F^l<i and Aj' = 8g^. The late Greek indices fi label the Minkowskian 
coordinates, while the early Greek indices a = 1, 2 and the upper case Latin indices 
^ = 3, 4 label the fields 0" = (0", (/.■*). 

The static Hamiltonian pertaining to the Lagrangian (4) is 

■H = \oF^j + AilA^I' + X2\DiirDj^,l>''\^ 

(5) 

+i2X,\d,A,p + Al{X,\r? + 16X2[\r?\dJ^? + i|a.(l0-*P)l']}, 

where the indices i = a,3 label the space-like coordinates. 

To find the static solutions of the model, one would usualy solve the Euler La- 
grange equations which minimise the Hamiltonian (5), but because of the electric 
potential Ao, one must solve the Euler Lagrange equations derived from the La- 
grangian (4). We then look for static solutions, but, as for the Julia-Zee dyon [9], we 
have to impose the proper asymptotic behaviour for the electric potential to obtain 
static solutions which are electrically charged (in the classical sense, i.e. solutions 
where the flux of the electric fleld is non zero). 

When the full equations of motion are written down, one flnds as expected that 
there are static solutions for which Ao — 0, i.e. solutions for which the electric fleld is 
identically zero. For these solutions in the temporal gauge, the equations of motion 
reduce to the equations obtained by minimising the Hamiltonian (5). We study the 
solutions of unit Baryon charge of the U{V] gauged Skyrnio model with and whithout 
an electric field, for various values of the C'^(l) coupling constant (or equivalently the 
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